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Abstract. We show the existence of global-in-time weak solutions to a general class of coupled 
bead-spring chain models that arise from the kinetic theory of dilute solutions of nonhomoge- 
neous polymeric liquids with noninteracting polymer chains, with finitely extensible nonlinear 
elastic (FENE) spring potentials. The class of models under consideration involves the unsteady 
incompressible Navier-Stokes equations with variable density and density-dependent dynamic 
viscosity in a bounded domain in R"^, o! = 2 or 3, for the density, the velocity and the pres- 
sure of the fluid, with an elastic extra-stress tensor appearing on the right-hand side in the 
momentum equation. The extra-stress tensor stems from the random movement of the poly- 
mer chains and is defined by the Kramers expression through the associated probability density 
function that satisfies a Fokker— Planck-type parabolic equation, a crucial feature of which is 
the presence of a centre-of-mass diffusion term and a nonlinear density-dependent drag coeffi- 
cient. We require no structural assumptions on the drag term in the Fokker— Planck equation; 
in particular, the drag term need not be corotational. With initial density po S [PmimPmax] 
for the continuity equation, where p^in > 0; a square-integrable and divergence-free initial ve- 
locity datum uq for the Navier-Stokes equation; and a nonnegative initial probability density 
function jpo for the Fokker-Planck equation, which has finite relative entropy with respect to 
the Maxwellian M associated with the spring potential in the model, we prove, via a limiting 
procedure on certain regularization parameters, the existence of a global-in-time weak solution 
t i—> {p{t),u(t),ip(t)) to the coupled Navier— Stokes-Fokker-Planck system, satisfying the initial 
condition {p{0),u{0),ip{0)) = (po, JiO: V'o), such that t M> p{t) G [pmin , Pmax] , t i-> u{t) belongs 
to the classical Leray space and t i— > tpit) has bounded relative entropy with respect to M and 
t I— >■ •il){t)/M has integrable Fisher information (w.r.t. the Gibbs measure dv := M{q)dq dx) 
over any time interval [0,T], T > 0. 

Keywords: Kinetic polymer models, FENE chain, Navier-Stokes-Fokkcr-Planck system, vari- 
able density, nonhomogeneous dilute polymer 



1. Introduction 

This paper establishes the existence of global-in-time weak solutions to a large class of bead- 
spring chain models with finitely extensible nonlinear elastic (FENE) type spring potentials, — 
a system of nonlinear partial differential equations that arises from the kinetic theory of dilute 
polymer solutions. The solvent is an incompressible, viscous, isothermal Newtonian fluid with 
variable density and viscosity confined to a bounded open Lipschitz domain 51 C ., d = 2 or 3, 
with boundary dO,. For the sake of simplicity of presentation, we shall suppose that has a 'solid 
boundary' 917; the velocity field u will then satisfy the no-slip boundary condition u = on dVL. 
The polymer chains, which are suspended in the solvent, are assumed not to interact with each 
other. The equations of continuity, balance of linear momentum and incompressibility then have 
the form of the incompressible Navier-Stokes equations with variable density and viscosity (cf. 
Antontsev, Kazhikhov & Monakhov 2|j Feireisl & Novotny [24 , Lions [39] or Simon [52]) in which 
the elastic extra-stress tensor r (i.e., the polymeric part of the Cauchy stress tensor) appears as a 
source term: 
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Given T e M>o, find p : {x,t) e n x [0, T] i-^ p{x, i) e M, w : (x, i) e O x [0, T] i-^ t) £ R"^ 
and p : (x, G X (0, T] H- p(a;, i) G M such that 

(1.1a) inr!x(0,r], 

(1.1b) p{x,o) = po{x) Wxen, 

d{pu) 

(1.1c) ^^+V, •(pii®w)-V, •(M(p)I)(ii)) + V,p = p/ + V, -r inl7x(0,r], 

(l.ld) V:r-w = inf7x(0,r], 

(l.le) u = on (9r2 X (0,T], 

(l.lf) (pw)(a;,0) = ipoUo)ix) Wx e Q. 

It is assumed that each of the equations above has been written in its nondimensional form; p 
denotes a nondimensional solvent density, u a nondimensional solvent velocity, defined as the 
velocity field scaled by the characteristic flow speed Uq. Here P{v) := \ {V x v + (Yx u)'^) ^^'^ 
rate of strain tensor, with {Yx'!i){x,t) G W^^*^ and {Yxv}ij = The scaled dynamic viscosity 
of the solvent, /i(-) G M>oi is density-dependent; in addition, p is the nondimensional pressure and 
/ is the nondimensional density of body forces. 

In a bead-spring chain model, consisting of iiT + 1 beads coupled with K elastic springs to 
represent a polymer chain, the extra-stress tensor r is defined by the Kramers expression as 
a weighted average of ip, the probability density functi on o f the (random) conformation vector 



q := {qj , . . . , qj^)'^ G R^'^ of the chain (see equation (1.81 below), with qi representing the d- 
component conformation/orientation vector of the ith spring. The Kolmogorov equation satisfied 
by "0 is a second-order parabolic equation, the Fokker-Planck equation, whose transport coefficients 
depend on the velocity field u, and the viscous drag coefficient appearing in the Fokker-Planck 
equation is a linear function of the dynamic viscosity /i (Stokes drag is assumed), which, in turn, is 
a nonlinear function of the density p. The domain D of admissible conformation vectors D C M.^'^ 
is a if-fold Cartesian product I?i x • • • x Dk of balanced convex open sets Di C M'', i — 1, . . . ,K; 
the term balanced means that qi G if, and only if, —qi G Di. Hence, in particular, G Di, 

i = 1, . . . , K . Typically Di is the whole of M'^ or a bounded open c?-dimensional ball centred at the 
origin G M'' for each i = 1, . . . ,K. When K ~ 1, the model is referred to as the dumbbell model. 
Let Oi C [0, oo) denote the image of Di under the mapping qi G Di h- > ^jqip, and consider 

the spring potential G C^Cj; K>o) n W^^^ {0^;M.>a), i = l,...,K. Clearly, G Oi. We shall 
suppose that Ui{Q) — and that Ui is unbounded on Oi for each i — 1,...,K. The elastic 
spring -force F, : D,C M^' of the ith spring in the chain is defined by 

(1.2) F,iq,):^U',C^\qJ^)q,, i = l,...,K. 

The partial Maxwellian Mi, associated with the spring potential Ui, is defined by 

The (total) Maxwellian in the model is then 

K 



K 

(1.3) M(g) := J]M,(9,) "iq -.^ {qj , . . . ,ql)^ £ D -.^ X 

~ i=i ~ ~ ~ ~ i^i 

Observe that, for i = 1, . . . , if, 

(1.4) Af(g)y,jAf(g)]-i - -[M(g)]-iy,,M(g) = y,,c/,(i|g,|2) = u[(\\qj-')q„ 

and, by definition, 

M{q)dq = 1. 
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Example 1.1. In the Hookean dumbbell model K = 1, and the spring force is defined by F{q) = q, 
with q E D = M'^, corresponding to U{s) = s, s E O = [0, oo). More generally, in a Hookean 
bead-spring chain model, K > 1, Fi{qi) = qi, corresponding to Ui{s) = s, i = 1, . . . , K , and D is 
the Cartesian product of K copies of M.'^. The associated Maxwellian is 

M{qJ = M,{q^)---MK{qK)^^e , 

with |gp := |qip + • • • + |<ZifP and Z := Zi ■ ■ ■ Z^ = (27r)^''/^. Hookean dumbbell and Hookean 
bead-spring chain models are physically unrealistic as they admit arbitrarily large extensions, o 

A more realistic class of models assumes that the springs in the bead-spring chain have finite 

extension: the domain D is then taken to be a Cartesian product of K bounded open balls Di cM."^', 

centred at the origin S K'^, i = 1, . . . , K, with K > 1. The spring potentials Ui : s E [0, bi/2) i— > 

Ui{s) E [0,(X)), with 6i > 0, i = 1, . . . , iiT, are in that case nonlinear and unbounded functions, 

and the associated bead-spring chain model is referred to as a FENE (finitely extensible nonlinear 

elastic) model; in the case of = 1, the corresponding model is called a FENE dumbbell model. 

Here we shall be concerned with finitely extensible nonlinear bead-spring chain models, with 
11 1 
D 5(0, ) X • • • X 5(0, &^), where bi > Q, i = I, . . . , K , and 5(0, b^ ) is a bounded open baU 

in of radius b^ , centred at S M''. We shall adopt the following structural hypotheses on the 
spring potentials Ui and the associated partial Maxwellians Mi, i = 1, . . . ,K. 

We shall suppose that for i = 1, . . . ,K there exist constants > 0, j = 1, 2, 3, 4, and 7i > 1 
such that the spring potential Ui satisfies 

(1.5a) Qi [dist(g„ dDi)P^ < M,{qi) < ca [dist(g„ dDi)^' Vg, E A, 

(1.5b) c,3 < [dist(g„ dDi)] t/;(i|g.|2) < c,i Vg, E A- 



Since [t^i(| |9iP)]^ = (— logMi(qi) -|- Const. )'^, it follows from (1.5a b) that (if 7^ > 1, as has been 



^2 

assumed here 



(1.6) / l + [C/,(ik,|2)]2 + [C/;(i|g,|2)]2 M,(q,0 dg, < 00, i = l,...,K. 



2 ' 



Example 1.2. In the FENE (finitely extensible nonlinear elastic) dumbbell model, introduced by 
Warner [55], K = 1 and the spring force is given by F{q) = (1 — \q\'^ /b)^^ q, q E D = 5(0,65), 

corresponding to U{s) = - | log (l - ^), s € O = [0, |), 6 > 2. More generally, in a FENE bead 
spring chain, one considers K + 1 beads linearly coupled with K springs, each with a FENE spring 
potential. Direct calculations show that the partial Maxwellians Mi a nd th e elastic potentials U, 
i = 1,. . . ,K, of the FENE bead spring chain satisfy the conditions (1.5a b) with and 7; 
provided that hi > 2, i = 1, . . . ,K . Thus, ( 1.6 1 also holds and bi > 2, i = 1, . . . ,K . 

It is interesting to note that in the (equivalent) stochastic version of the FENE dumbbell model 
{K = 1) a solution to the system of stochastic differential equations associated with the Fokker- 
Planck equation exists and has trajectorial uniqueness if, and only if, 7 = | > 1; (cf. Jourdain, 
Lelievre & Le Bris [5^ for details). Thus, in the general class of FENE-type bead-spring chain 
models considered here, the assumption 7^ > 1, i = 1, . . . , ii', is the weakest reasonable requirement 
on the decay-rate of Mi in (1.5a| as <l\st{qi,dDi) — > 0. 

Another example of a finitely extensible nonlinear elastic potential is Cohen's Fade approximant 
to the inverse Langevin function (CPAIL) (cf. [2]), with 



(1.7) 



U^{S) 



3 3 



1 - — 



and 



F^{q^) 



36. 



1 



19.1 



where bi is again a positive parameter. Again, direct calculations show that t he M axwellian M and 
the elastic potential U of the CPAIL dumbell model satisfy the conditions ( 1.5a b) with 7^ y , 
provided that 5^ > 3, i = 1, . . . , i^. Thus, (|1.6| also holds and bi> i, i = 1, . . . , K. 
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We note in passing that both of these force laws are approximations to the inverse Langevin 
force law 



F^{q^ 




where the Langevin function L is defined by L{t) :— coth(i) 
monotonia increasing on [0,oo) and tends to 1 as its argument tends to oo, it foUows that the 
function \qi\ € [0, %/5i ) '"^ L^^{\Qi\/Vbi ) G [0,cx)) is strictly monotonic increasing, with a vertical 
asymptote at \qi\ — \fbi. o 



1/t on [0, cx)). As L is strictly 



The governing equations of the general nonhomogeneous bead-spring chain models with centre- 
of-mass diffusion considered in this paper are ( |l.lc| -d), where the extra-stress tensor r is defined 

by the Kramers expression: 



(1.8) 



t{x, t) — k 



^^1^(x,g,^)g^g?C/;(^k^P) dq-KQ{x,t)l\ , 



with the density of polymer chains located at x at time t given by 



(1.9) 



i>{x,q,t) Aq, 



D 



(not to be confused with the solvent density p(x,t)). The probability density function i/j is a 
solution of the Fokker-Planck (forward Kolmogorov) equation 

K 



d->jj 

1h 



+ {U ■ V^)V + Y«, • (^(^) <l^ V') 



(1.10) 



Cip) 



1 

4A 



K K 

j=i 



\ap)M 



in f7 X £) X (0,T], 



with a{v) = Va; V and a density-dependent scaled drag coefficient (■(•) G M>o. A concise derivation 
of the Fokker-Planck equation ( |1.10 1 is presented below. 

The nondimensional constant fc > featuring in (1.8) is a constant multiple of the product of 



the polymer number density (the number of polymer molecules per unit volume), the Boltzmann 
constant fcs, and the absolute temperature T. In (1.10 1, e > is the centre-of-mass diffusion 
coefficient defined as e := (fo/^o)^/(4(if-)- 1)A) with Lq a characteristic length-scale of the solvent 
flow, ^0 '■— v^fcsT/H signifying the characteristic microscopic length-scale and A :— (Co/4H)(C/o/io), 
where Co > is a characteristic drag coefficient and H > is a spring-constant. The nondimensional 
parameter A S M>o, called the Weissenberg number (and usually denoted by Wi), characterizes 
the elastic relaxation property of the fluid, and A — (^ij)fj^i is the symmetric positive deflnite 
Rouse matrix, or connectivity matrix; for example, A — tridiag [—1, 2, —1] in the case of a linear 
chain; see, Nitta [45]. Concerning these scalings and notational conventions, we remark that the 



factor 1/(4A) in equation ( |1.10 ) above appears as a factor 1/(2 A) in the Fokker-Planck equation 
in our earlier papers [HI liot 1111 11^ . 



Definition 1.1. The collection of equations and structural hypotheses (l.la-f) -( [l.lO l will be re- 
ferred to throughout the paper as model (P), or as the general nonhomogeneous FENE-type bead- 
spring chain model with centre-of-mass diffusion. 



We continue with a brief derivation of the Fokker-Planck equation (1.10) in the general non- 



homogeneous FENE-type bead-spring chain model (P) with centre-of-mass diffusion. Henceforth, 



up to and including equation (1.19) below where our nondimensionalization is introduced, the 
symbols p, w, ?/;, g, /i(p), C(p)j Hi, F^, x, q, t will refer to the dimensional forms of these functions 
and variables, unless otherwise stated. For the sake of simplicity of the exposition, we shall confine 
ourselves to the case of linear chains, when the Rouse matrix A — tridiag [—1, 2, —1]; we empha- 
size however that the results that we prove later on in the paper are completely independent of 
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Figure 1. Rouse chain with K springs and K + \ beads. 



the actual structure of A: we shah only require that A is symmetric and positive definite. Thus, 
each polymer molecule is modelled by a linear arrangement oi K -\- \ beads with mass ni that 
are joined by K massless elastic springs. For each K G N>i, the state of a bead-spring chain in 
such an arrangement is fully determined by the position vectors of \ts K + 1 beads, denoted by 
for £ = 1, . . . , K + 1. Equivalently the state of a bead-spring chain can be described by the 
center of mass of the bead system, Tc := Sfci^ Zi^ together with the K connector vectors 
Qi :— Tj+i — ri, i — 1, . . . , K , which are to be understood as column vectors in E'' (cf. Fig. jlj. 

In what follows we shall write r := {rj, . . . , Zk+i)'^ 9 ■= ilJ j ■ ■ ■ : lli)'^ ■ Each spring exerts 
an elastic conservative force on the beads it connects along the corresponding connector vector 
and has a magnitude that depends isotropically on it. We model the spring force by a spring force 
function Fi, which has the form Fi[qi) = HC/^(i|gip) g^, where H > is a spring constant, and 
Ui is a spring potential, i — 1, . . . ,K . It will be assumed that \qi\ < qf^'^^, where q^^^ denotes 
the maximum admissible length of the ith spring. We let Di := {qi E M'^ : \qi\ < qf^^^}. It is 
immediate that Fi{qi) ~ —F{-~qi) for all qi e -Di, i — 1, . . . ,K. 

In the absence of external forces and neglecting inertial effects the Langevin equation for the 
£-th bead in this model is 

K 

(1.11) Co [drj - u{rj, ■) dt^ ^ Bidt + Y,GHF,{q,)dt, e=l,...,K + l. 

i=l ~ ~ 

Here Co is a characteristic drag coefficient, which we shall assume for the moment to be a fixed 
positive constant; Bi denotes a Brownian force acting on the ith bead; and the (K + l) x K matrix 
G with entries Gij, i = I, . . . ^ K + 1, j = 1, . . . , K, called the graph incidence matrix, which (in 
our case here of a linear bead-spring chain with A' + 1 beads and K springs) is defined by 

/ 1 \ 

^11 1+1 if spring j starts from bead z, 

G := ■.. G IR'^^+^^^^; i.e., Gij .= < —1 if spring j teminates in bead i, 

_l I I otherwise. 

The Rouse matrix A G M.^^^ is related to the matrix G through the equality A = G^G. The 
form of the matrix G reflects the fact that the £th bead is pulled by the (£ — l)st spring in the 
—qi^i = r£-i — re direction and by the £th spring in the qe = Zt+i ~ direction with proper 
provision for the beads at the ends of the chain, which are only pulled by a single spring each. Note 
that qi ^ — Gii rt,, i = 1,. . . ,K. For £ — 1, . . . ,K + 1, the Brownian force is defined 
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by a d-component vectorial Wiener process Wg via Bgdi = y2fcj3Y(o dW^, where fc^ 1.38 x 
-[^Q-23 ni^kgs^^K"^ is Boltzmann's constant and T is the absolute temperature measured in Kelvin, 
K. The coefficient V^fcs T Co is due to the Einstein-Smoluchowski relation, which determines the 
diffusion coefficient in Brownian motion. Therefore, ( |1.11| can be written in the following form 



dt 

















= < 




+ Co"'G 




> 






. w(?:a'+i(*),0 . 




_ EK{ZK+lit)-ZK{t)) _ 


> 



(1.12) 




0^1 w 



By defining 



m ■■= 



lit) 



Wit) 



Wlit) 

W2it) 
WK+lit) 




kiXit)) 



uiziit),t) 

Uit2it),t) 



Fiiz2it) 
^2(2:3(0 



ziit)) 

t2it)) 



uitK+lit).t) J [ EKiZK+lit)-ZKit)) 

equation (1.121 can be rewritten as the following stochastic differential equation: 

(1-13) dX{t) ^b{X{t))dt + a{X{t))dW{t), X{Q)^X, 

for the {K + l)(i-component vectorial random variable X{t), t e [0,r]. 

By recalling a classical result from stochastic analysis stated in Theorem 1 1 . 1| below, we can now 
write down the associated forward Kolmogorov equation (Fokker-Planck equation), a parabolic 
partial differential equation governing the evolution of the probability density function of the 
stochastic process 1 1-> X{t) (see, for example. Corollary 5.2.10 in |3S]). 

Theorem 1.1. Suppose that the (K + l)d-component vectorial random variable X{t) has a prob- 
ability density function {z,t) i— > il)[z,t) of class C'^'^iW'-^'^^^'^ x [0,cx))) (i.e., twice continuously 
differentiable with respect to z & and once with respect to t), and let X{Q) = X be 

a square-integrable random variable with probability density function tpQ S C^(M^^^^-''^). Also, 
suppose that b and a in (1.13) are globally Lipschitz continuous, and c{z) — a{z) a{z)'^ . Then, 

(K+l)d 



{K+l)d ... , 

j=i ■> j,j=i 
in X [0,00) where ^'(5,0) = 7/'o(g) for z e 



dzidz.j 



(ctji^), 



The Hookean spring force satisfies the global Lipschitz continuity assumption in Theorem |1.1[ 
whereas FENE-type spring forces do not. Indeed, FENE-type forces are only locally Lipschitz on 
D and are not defined on all of M'*. For the purposes of the present informal derivation of equation 
(1.10) we shall ignore this technicality, and will proceed as if Theorem 1 1 . 1 1 held in the FENE case 
also. Thus, on applying Theorem |1.1| to (1.131, we arrive at the following forward Kolmogorov 
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equation: 



d'4' 



K+l 



dt 



(1.15) 



K+l 



A' 



Co 



1,A'+1 



where ip^'^+^ denotes the probability density function with respect to the {r^, . 
nates. By performing the linear change of variables 

r (ri, . . . , tk+i) G K^^+i)'* ^ (r„ q) := (r„ gi, . . . , g^) G K^^+i)'^ 
and letting Qi ■= i 9)' obtain 



-A+ 



^) coordi- 



dtp 

Ik 



a: 



K 



1 \ 

£=1 / 



i=l 




A'+l A 



A A 



(1.16) 

\ /\ -I- 1 [ n / 

i=l i = l 

where A = G^G is the Rouse matrix. We shall use the approximations 

A+l A+l 
> r II . ) ft^ ii.(r „ . t\ . 

K 



^ A+l A+l 

e=i i=i 



for 



I,. . . ,K, in the second and the third term on the left-hand side of (1.16), respectively, 



and suppose that the associated approximation errors are negligible; we note in particular that 
if u is linear in its spatial variable, then the above approximations are exact. Otherwise, the 
approximation errors are nonzero and their magnitudes depend on the extent to which x i— >■ u{x,-) 
deviates from a linear function on the characteristic microscopic length scale ig. Hence, 



'dt 



K I A , 

v.. • (y(rc,-)^) + E^9- • {'^ru{rc,■)q^)^-Y,A,,F,{q,)'^ 



(1.17) 



More generally, since the (Stokes) drag coefficient depends linearly on the dynamic viscosity, 
which, in turn, has been assumed to depend nonlinearly on the (variable) density, we shall replace 
Co in (1.17) with CoC(p(Sji))i where C(p) is a smooth nondimensional function of the (unknown) 
nondimensionalized density p. After renaming the centre of mass Vc into x, we thus have that 



9-0 
'dt 



A / ^ K ,\ 



(1.18) 



knT 1 



Co K+l 



A, 



Co 



K A 



EE^»^-^^ 

i=i j=i 



which, in the special case of (^{p{x, t)) = 1 collapses to (1.17). Concerning alternative models with 
variable, configuration-dependent drag, CoC('?)i where C is a certain fixed function of g, we refer 
to the comments at the end of this section. 
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The quantity g{x,t) = J^ip{x,q,t) dq, signifying the density of polymer chains at the point 
X £ n and time t G [0, T], has units of length to the power minus d. We shall assume temporarily, 
for the purposes of nondimensionalization, that there are no initial polymer molecule concentration 
gradients; i.e., that ^(^O) is constant throughout the spatial domain f2; we denote this constant, 
called the polymer density number, by Up. It then follows on integrating (1.18) over D, using 
the divergence theorem and dropping boundary integrals (cf. (1.24a b) below) that g is constant 
through Q and throughout time. We emphasize, that our analysis in subsequent sections does not 
require that g is constant in space or in time: we shall only demand that g{-,0) ~ V'o(': ?) dq G 
L°°(S7), in fact (cf. (3.3)). In any case, with Up thus defined, we introduce nondimensionalized 



(hatted) variables in terms of their nonhatted counterparts as follows: 

(1.19) x::^Lqx, qi'.^ioq., u -.^ Uqu, t -.^ {Lo/Uo)% and 4' := {up/ig'^) ip, 



where £o := \/fcs T/H is the characteristic length scale of a spring and Lq and Uq are the character- 
istic macroscopic length and velocity, respectively. Upon nondimensionalization using the hatted 



variables, the Fokker-Planck equation (1.18) becomes 



K 



1 



K 



1 

4A {K 



1) 



3 = 

K K 




where A := {(o/AH) (Uo/Lq) is the characteristic relaxation time of a spring (the Weissenberg 
number, usually denoted by Wi: the ratio of the microscopic to macroscopic time scales), 

u^{s)^e^^u,iils), F^{q,) = L^;(i|g.n?. = (H4)-'F,(4?.), 



and the spatial, configurational and temporal variables x, q and t now belong to the rescaled 
domains 

n:=n/Lo, D:=D/eo, and [0,f] with f := C/oT/io, 
respectively. Discarding the hats and writing e for the centre-of-mass diffusion coefficient; i.e., 

2 

e 



4:\{K+1) 



we have 



~dt 



K 

4=1 



{^xu)qi^- 



1 ^ 
4A ^ 



A^,FM3) 




K K 



4 = 1 j = l 



Finally, on recalling the definition (1.3) of the Maxwellian M and noting the identities (1.4), the 
two terms involving the entries Aij of the Rouse matrix A in the last equation can be merged into 



a single term, resulting in the Fokker-Planck equation (1.10 1. From this nondimensionalization 



procedure one also obtains the important nondimensional parameter 



which measures how the maximal admissible extension qf^^^ of the ith spring in the chain compares 
with the characteristic microscopic length-scale ^o- Having defined bi we can express the non- 
dimensionalized configuration domain for the ith spring as Di — {qi £ : \qi\ < \/bi]. 



A noteworthy feature of equation (1.10) in the model (P) compared to classical Fokker-Planck 
equations for bead-spring models in the literature is the presence of the x-dissipative centre-of- 
mass diffusion term eA^,?/' on the right-hand side of the Fokker-Planck equation (1.10). We 
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refer to Barrett & Siili [7] for the derivatfon of (1.10 1 in the case oi K = 1 and constant p; 
see also the article by Schieber |49| concerning generalized dumbbell models with centre-of-mass 
diffusion, and the recent paper of Degond & Liu [IHl for a careful justification of the presence of 
the centre-of-mass diffusion term through asymptotic analysis. In standard derivations of bead- 
spring models the centre-of-mass diffusion term is routinely omitted on the grounds that it is 
several orders of magnitude smaller than the other terms in the equation. Indeed, when the 
characteristic macroscopic length-scale Lq « 1, (for example, Lq = diam(ri)), Bhave, Armstrong 
& Brown [T3] estimate the ratio £q/Lq to be in the range of about 10~^ to 10"''. However, the 
omission of the term e A^tp from (1.10) in the case of a heterogeneous solvent velocity u{x,t) 
is a mathematically counterproductive model reduction. When e Axij' is absent, (1.10) becomes 
a degenerate parabolic equation exhibiting hyperbolic behaviour with respect to {x,t). Since 
the study of weak solutions to the coupled problem requires one to work with velocity fields u 
that have very limited Sobolev regularity (typically u S L'^{0,T; L'^{n)) n L'^{0,T; Hl{n))), one 
is then forced into the technically unpleasant framework of hyperbolically degenerate parabolic 
equations with rough transport coefficients (cf. Ambrosio |1J, DiPerna & Lions [20], Mucha 
The resulting difficulties are further exacerbated by the fact that, when D is bounded, a typical 
spring force F{q) for a finitely extensible model (such as FENE) explodes as q approaches dD; 
see Example 1 1 . 2 1 above . For these reasons, as in our earlier papers in this field (cf. [7] [8l iQj ITO]). 
we shall retain the centre-of-mass diffusion term in (1.10 1. 

We continue with a brief literature survey. Unless otherwise stated, the centre-of-mass diffusion 
term is absent from the model considered in the cited reference (i.e., e is set to 0) and the density p 
of the solvent is assumed to be constant; also, in all references cited K = 1, i.e., a simple dumbbell 
model is considered rather than a bead-spring chain model. 

An early contribution to the existence and uniqueness of local-in-time solutions to a family of 
dumbbell type polymeric flow models is due to Renardy |3H]- While the class of potentials F{q) 
considered by Renardy [35^ (cf. hypotheses (F) and (F') on pp. 314-315) does include the case of 
Hookean dumbbells, it excludes the practically relevant case of the FENE dumbbell model (see 
Example 1.2 above). More recently, E, Li & Zhang ,22, and Li, Zhang & Zhang [37j have revisited 
the question of local existence of solutions for dumbbell models. A further development in this 
direction is the work of Zhang & Zhang [S5], where the local existence of regular solutions to 
FENE-type dumbbell models has been shown. All of these papers require high regularity of the 
initial data. Constantin |15j considered the Navier-Stokes equations coupled to nonlinear Fokker- 
Planck equations describing the evolution of the probability distribution of the particles interacting 
with the fluid. Otto fc Tzavaras [47] investigated the Doi model (which is similar to a Hookean 
model (cf. Example 1 1 . 1 1 above) , except that D = S^) for suspensions of rod- like molecules in the 
dilute regime. Jourdain, Lelievre & Le Bris [53] studied the existence of solutions to the FENE 
dumbbell model in the case of a simple Couette flow. By using tools from the theory of stochastic 
differential equations, they showed the existence of a unique local-in-time solution to the FENE 
dumbbell model for d = 2 when the velocity field u is unidirectional and of the particular form 

U{X1,X2) = (1X1(0:2), 0)'^. 

In the case of Hookean dumbbells (K — 1), and assuming e = and constant solvent density 
p, the coupled microscopic-macroscopic model described above yields, formally, taking the second 
moment of g 1— >■ '0('?jS;^)i the fully macroscopic, Oldroyd-B model of viscoelastic flow. Lions & 
Masmoudi [ID] have shown the existence of global-in-time weak solutions to the Oldroyd-B model 
in a simplified corotational setting (i.e., with a{u) = YxU replaced by -^(Vxy: — {Yxu)"^)) by 
exploiting the propagation in time of the compactness of the solution (i.e., the property that if one 
takes a sequence of weak solutions that converges weakly and such that the corresponding sequence 
of initial data converges strongly, then the weak limit is also a solution) and the DiPerna-Lions 
[20| theory of renormalized solutions to linear hyperbolic equations with nonsmooth transport 
coefficients. It is not known if an identical global existence result for the Oldroyd-B model also 
holds in the absence of the crucial assumption that the drag term is corotational. With £ > and 
constant solvent density p, the coupled microscopic-macroscopic model above yields, taking the 
appropriate moments in the case of Hookean dumbbells, a dissipative version of the Oldroyd-B 
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model. In this sense, the Hookean dumbbell model has a macroscopic closure: it is the Oldroyd-B 
model when e = 0, and a dissipative version of Oldroyd-B when e > (cf. Barrett & Siili [7j)- 
Barrett & Boyaval [1] have proved a global existence result for this dissipative Oldroyd-B model in 
two space dimensions. In contrast, the FENE model is not known to have an exact closure at the 
macroscopic level, though Du, Yu & Liu [21] and Yu, Du & Liu [57] have recently considered the 
analysis of approximate closures of the FENE dumbbell model. Lions & Masmoudi [H] proved 
the global existence of weak solutions for the corotational FENE dumbbell model, once again 
corresponding to the case of e = 0, constant solvent density p, and K = 1, and the Doi model, also 
called the rod model; see also the work of Masmoudi [15] . Recently, Masmoudi [H] has extended 
this analysis to the noncorotational case. 

Previously, El-Kareh & Leal '23' had proposed a steady macroscopic model, with added dissi- 
pation in the equation satisfied by the conformation tensor, defined as 



P{x) / qq^U'il\q\^),Pix,q)dq, 



ID~" 

in order to account for Brownian motion across streamlines; the model can be thought of as an 
approximate macroscopic closure of a FENE-type micro-macro model with centre-of-mass diffu- 
sion. 

Barrett, Schwab & Siili 6 showed the existence of global weak solutions to the coupled 



microscopic- macroscopic model (l.lc-d), (1.10) with e = 0, K = 1, constant solvent-density 



p, an x-mollified velocity gradient in the Fokker-Planck equation and an x-moUified probability 
density function ip in the Kramers expression, admitting a large class of potentials U (including 
the Hookean dumbbell model and general FENE-type dumbbell models) ; in addition to these mol- 
lifications, u in the x-convective term (u • Yx)''!' in th^ Fokker-Planck equation was also mollified. 
Unlike Lions & Masmoudi [10], the arguments in Barrett, Schwab & Siili [6 did not require that 
the drag term Vg • {(7{u) q tp) in the Fokker-Planck equation was corotational in the FENE case. 

In Barrett & Siili [^ , we derived the coupled Navier-Stokes-Fokker-Planck model with centre- 
of-mass diffusion stated above, in the case of A' = 1 and constant solvent-density p. We established 
the existence of global-in-timc weak solutions to a mollification of the model for a general class 
of spring-force-potentials including in particular the FENE potential. We justified also, through 
a rigorous limiting process, certain classical reductions of this model appearing in the literature 
that exclude the centre-of-mass diffusion term from the Fokker-Planck equation on the grounds 
that the diffusion coefficient is small relative to other coefficients featuring in the equation. In the 
case of a corotational drag term we performed a rigorous passage to the limit as the mollifiers in 
the Kramers expression and the drag term converge to identity operators. 

In Barrett & Siili [8] we showed the existence of global-in-time weak solutions to the general 
class of noncorotational FENE type dumbbell models (including the standard FENE dumbbell 
model) with centre-of-mass diffusion, in the case oi K ^ 1 and constant solvent-density p, with 
microsropic cut-off (cf. ( 1.21[ ) and (1.22) below) in the drag term 

(1-20) y,-(g(u)g^) = y,- a(y)gA/C(p) 

Subsequently, in [3] [TO], we removed the presence of the cut-off by passing to the limit i — ^ oo, 
with K > I, and the solvent density p, the viscosity /i and the drag coefficient kept constant. 

In this paper we prove the existence of global-in-time weak solutions to FENE-type models 
without cut-off or mollification, in the general case of K > 1 and with variable solvent-density /o, 
variable viscosity p,{p) and variable drag C(p)- This is achieved by replacing the use of Dubinskii's 
compactness theorem in ^ with the application of the Div-Curl lemma in our proof of relative 
compactness of the sequence of approximating solutions to the Fokker-Planck equation in the 
Maxwellian- weighted space L^{0, T; L\f{n x D)). Since the argument is long and technical, we 
give a brief overview of the main steps of the proof. 

Step 1. Following the approach in Barrett & Siili [8l|9l[T0] and motivated by recent papers of 
Jourdain, Lelievre, Le Bris & Otto 30] and Lin, Liu & Zhang [38| (see also Arnold, Markowich, 
Toscani & Unterreiter [3j, and Desvillettes & Villani jl9j) concerning the convergence of the 
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probability density function ip to its equilibrium value ■i/'ooCs, <z) (corresponding to the 

equilibrium value yoo(s) '■= Q of the velocity field in the case of constant density) in the absence 
of body forces /, we observe that if ip/{({p) M) is bounded above then, for L S ]R>o sufficiently 
large, the drag term (1.201 is equal to 

y,- e(«)(?AfC(p)/3'' 



(1.21) 



where <E C(M) is a cut-off function defined as 



dp) M 



(1.22) 



i{s,L). 



More generally, in the case of -ftr > 1, in analogy with ( 1.21[ ), the drag term with cut-off is defined 

by 



K 



\Cip)M 



It then follows that, for L ^ 1, any solution ip of (1.10), such that 'i/'/(C(p) is bounded above 
by L, also satisfies 

[u ■ y,)^ + g(y) q^ m c(p) 



'di 



i=l 



C(p)M 



(1.23) 



in X £» X (0, T] 



Let 5 A := i'l X 
initial conditions: 



X -Di-i X X Di+i X 



C(P) 

X Dk- We impose the following boundary and 



K 



i=i 



U(P)M 



a(w)g,MC(p) 



U(P)M 



= 



(1.24a) 
(1.24b) 



on 17 X X (0,T], for i = 1, . . . 
on dflx D X (0,T], 



(1.24c) V(-, -,0) = A/(-) C(po(-))/3^(V'o(-, •)/(C(Po(-)) M(-))) > on^xD, 

where qi is normal to dDi, as -D^ is a bounded ball centred at the origin, and n is normal to d^l. 
The initial datum ipQ for the Fokker-Planck equation is nonnegative, defined on x _D, with 



'^o{x,q)dqdx = 1, 



and assumed to have finite relative entropy with respect to the Maxwellian M; i.e. 

, , , . Mx,q) 
V'o(s, q) log ... . dq dx < oo. 
nx_D ~ M{q) ~ 

As we shall suppose throughout that the range of the function C is a compact subinterval [Cmin, Cmax] 
of (0, oo), the finiteness of the relative entropy with respect to the Maxwellian M is equivalent to 
demanding that 

Mx,q) , tpo{x,q) 

log — — ^ ^ dg dx < cx) . 



InxD C(po(x)) ° CiPoix)) M{q) 

Clearly, if there exists L>0 such that < < ^((Po) M, then M ((po) /?^(-0o/(C(Po) *^)) = ^o- 
Henceforth L > 1 is assumed. 



Definition 1.2. T/ie coMpZed problem ( |l.la| -f), ( |l.8[ ), ( |l.9[ ), ( |1.23[ ), ( |l.24a| -c) wi^/ &e referred 
to as model (Pl), or as the general nonhomogeneous FENE-type bead-spring chain model with 
centre-of-mass diffusion and microscopic cut-off, with cut-off parameter L > 1. 
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In order to highlight the dependence on L, in subsequent sections the sohition to ( 1.23 1, ( 1.24a 
c) win be labelled ipL- Because of the coupling of ( 1.23 ) to ( 1.1c I through ( 1.8 1, the density, velocity 



and the pressure will also depend L and we shall therefore denote them in subsequent sections by 
PL, ul and pl- 

The cut-off has a convenient property: when the solvent density p is constant, the couple 
(moojV'oo), defined by Uoo(s) •= Q ^-^d ipooix,q) := M{q), which is an equilibrium solution of the 



hen ( 1.10 1 is rep 



system ( 1.1c -d), (|l.8|), (|l.9|), (|1.10|), with / = 0, is still an equilibrium solution of the cut-off version 
of the problem. 



aced by (1.23) (with boundary and initial conditions ( 1.24a -c)), 
for all L > 0. Thus, unlike the truncation of the (unbounded) spring-potential proposed in the work 
of El-Kareh & Leal [53] in the case of c onstan t-density flows, the introduction of the cut-off function 
into the Fokker-Planck equation (1.10 1 does not alter the equilibrium solution (uocV'oo) of 



the original Navier-Stokes-Fokker-Planck system. In addition, the boundary conditions for on 
dftx Dx (0, T] and n x dD x (0, T] ensure that 



flxD 



4'ix,q,0) dq = 1 



for a.e. t € K>o, in agreement with the requirement that "0 is a probability density function. 

Step 2. Ideally, one would like to pass to the limit L — >• oo in problem (P^) to deduce the 
existence of solutions to (P). Unfortunately, such a direct attack at the problem is (except in the 
special case of d = 2, or in the absence of convection terms from the model,) fraught with technical 
difhculties. Instead, we shall first (semi)discretize problem (Pl) by an implicit Euler scheme with 
respect to t, with step size At; this then results in a time-discrete version (P^*) of (Pl)- By using 
Schauder's fixed point theorem, we will show in Section [s] the existence of solutions to (Pl*)- In 
the course of the proof, for technical reasons, a further cut-off, now from below, is required, with 
a cut-off parameter S £ (0, 1), which we shall let pass to to complete the proof of existence of 
solutions to (Pl*) in the limit of 5 — )■ 0+ (cf. Section [s]). Ultimately, of course, our aim is to show 
existence of weak solutions to the general nonhomogeneous FENE-type bead-spring chain model 
with centre-of-mass diffusion, (P), and that demands passing to the limits At 0+ and L oo; 
this then brings us to the next step in our argument. 

Step 3. We shall link the time step At to the cut-off parameter L > 1 by demanding that 
At = o(L^^), as L — > oo, so that the only parameter in the problem (P^*) is the cut-off parameter 
(the centre-of-mass diffusion parameter e being fixed). We shall show that p^* can be bounded, 
independent of the cut-off parameter L, in L°° {0,T; L°° (il)). By using special energy estimates, 
based on testing the Fokker-Planck equation in (P^*) with the derivative of the relative entropy 
with respect to the Maxwellian of the general nonhomogeneous FENE-type bead-spring chain 
model, we show that m^* can also be bounded, independent of L. Specifically, u^* is bounded 
in the norm of the classical Leray space, independent of L; also, the L°° norm in time of the 
relative entropy of '0l*/C(Pl') ^^'^ the norm in time of the Fisher information of ■j/'l* ■= 
V'l*/(C(Pl*) ^) bounded, independent of L. We then use these L-independent bounds on the 
relative entropy and the Fisher information to derive an L-independent bound on a fractional-order 
in time Nikol'skii norm of u^*. 

Step 4- The collection of L-independent bounds from Step 3, then enables us to extract a 
weakly convergent subsequence of solutions to problem (Pl*) as L oo; and then further strongly 
convergent subsequences {w^*}l>i and {pl*}l>i- The extraction of a strongly convergent sub- 
sequence from the weakly convergent sequence {iPl*}l>i is considerably more complicated: after 
some technical preparation, we apply the Div-Curl lemma, to finally obtain a strongly conver- 
gent subsequence of solutions {p^l'' 



iff*'') to (P^*) with At = o(L-i) as L 



oo, m 

L°°{0,T;LP{n)) X L^{0,T;L'-{n)) x LP{0,T]L\j{n x D)) for any p g [l,oo); any r e [l,oo) 
when d ~ 2 and any r G [1,6) when d ~ 3; enabling us to pass to the limit with the microscopic 
cut-off parameter L in the model (Pl*), with At = o(L^^), as L — > oo, to finally deduce the 
existence of a weak solution to model (P), the general nonhomogeneous FENE-type bead-spring 
chain models with centre-of-mass diffusion. 
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The paper is structured as follows. We begin, in Section [2j by stating (Pl), the coupled non- 
homogeneous Navier-Stokes-Fokker-Planck system with centre-of-mass diffusion and microscopic 
cut-off for a general class of FENE-type spring potentials. In Section [3] we establish the existence 
of solutions to the time-discrete problem (P^*). In Section [4] we derive an L-independent bound 
on the solvent density p^* in L°° {0,T; L°° (V,)); we also derive a set of L-independent bounds on 
M^* in the classical Leray space, together with L-independent bounds on the relative entropy of 
ipL* /C{Pl*) 'with respect to the Maxwellian M, and the L^ norm in time of the Fisher information 
of i/'l* := '0lV(C(Pl*) -^^)- We then use these L-independent bounds on spatial norms to show 
that the Nikol'skii norm 7V^(0, T; L^(f^)) of * is bounded, independent of L and At = o(L-i), 
for a suitable value of 7 e (0, 1). This allows us to prove, via Simon's extension of the Aubin-Lions 
compactness theorem (cf. [SI]), strong-convergence of the sequence {u^'^}l>i in L'^{0,T;L'^{n)) 
for r G [1,00) when d = 2 and r G [1,6) when d = 3. We then use this strong convergence result 
together with the DiPerna-Lions theory of renormalized solutions to linear transport equations 
with nonsmooth transport velocities to deduce the strong convergence of the sequence of approx- 
imate densities {p^*}l>i, and pass to the limit in our approximation to the continuity equation, 
as L — >■ 00, with At = o(L^^). Weak convergence of the sequence {V-'L*}i>i the Maxwellian- 
weighted space L^(0, T; L\,f{il x D)) is an immediate consequence of our entropy estimate, via 
de la Vallee-Poussin's theorem and the Dunford-Pettis theorem. The proof of the strong conver- 
gence of the sequence is however considerably more complicated; it is established in Section [4.5[ by 
first developing interior estimates in standard (unweighted) Lebesgue and Sobolev norms, exploit- 
ing the fact that on nonempty open relatively compact subsets of D the Maxwellian is bounded 
above and below by positive constants. We then use these interior estimates in conjunction with 
the Div-Curl lemma to deduce weak convergence of the sequence (1 -I- V'l *)"^^ on nonempty open 
relatively compact subsets of (0,r) x il x D, where a G (0, 1). Thus we can make use of the fact 
that the continuous functions s e [0,oo) i-> (1 -I- s)^"*"" and s G [0,oo) 1— )■ s" are, respectively, 
strictly convex and strictly concave, and therefore weakly lower (respectively, upper) semicon- 
tinuous, to deduce that {V'L*}i>i converges to a limiting function G L^{0,T; L\f{{l x D)), 
almost everywhere on compact subsets of (0, T) x ft x D; hence, by using a nested sequence of 
nonempty open relatively compact sets, we show that {iPl'^}l>i converges to tp almost everywhere 
on (0, T) X X D. Thanks to the fact that M{q) dq is a probability measure on D, and therefore 
M{q) dq dx dt is a finite measure of (0, T) x f2 x L>, Egoroff 's theorem then implies almost uniform 
convergence of the sequence, and therefore also convergence in measure; thus we can appeal to 
Vitali's theorem to finally deduce strong convergence in L^{0,T; L\j{il x D)) of (a subsequence 
of) the sequence {V'L*}i>i- This strong convergence result then allows us in Sectionjsjto pass to 
the limit with the cut-off parameter L in problem (P^*), with At ~ o{L^^), as L — >■ 00, to deduce 
the existence of a weak solution {p,u,ip M({p)ip) to problem (P), the general nonhomoge- 
neous FENE-type bead-spring chain models with centre-of-mass diffusion. We shall operate within 
Maxwellian-weighted Sobolev spaces, which provide the natural functional-analytic framework for 
the problem. Our proofs require special density and embedding results in these spaces that are 
proved, respectively, in Appendix C and Appendix D of which is an extended version of our 
paper [9J for FENE-type models. 

Our techniques can be easily modified to prove large-data global existence of weak solutions to 
kinetic models with configuration-dependent drag, where instead of being a nonlinear function of 
the unknown density p, as in (1.10), the drag coefficient ^ is a given function of q, bounded 
above and below by positive constants (cf. Hinch [37], Larson [33], Schroder et al. [SD] and 
references therein). The idea behind these models, which have been developed to explain physical 
mechanism by which large stresses rapidly build up in dilute polymer solutions, is that of a 
bead friction coefficient that depends strongly on the inter-bead distance through a nonlinear 
friction law. This principle of conformation-dependent hydrodynamic drag assumes that as a 
chain becomes extended by the flow, the strength of the hydrodynamic friction on the chain will 
also increase; see, [15] for a detailed survey. 
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2. The polymer model (P^) 

Let C M'' be a bounded open set with a Lipschitz-continuous boundary dH,, and supp ose th at 
the set D := Di X ■ ■ ■ X Dk of admissible conformation vectors q :— {qj , . . . ,qJc)"^ in (1-10) is 
such that Di, i — 1, . . . , K , is a.n open ball in M'*, = 2 or 3, centred at the origin with boundary 
dDi and radius ^/bi, bi > 2; let 

K 

(2.1) dD:=\JdD^, where := Z^i x • • • x A-i x 9A: x A+i x • • • x D^. 



l.la-f), dl.SD, ( 


1.9 


), ( 


1.231 and (1.24a- 



boundary- value problem, dependent on the parameter L > 1. As has been already emphasized 
in the Introduction, the centre-of-mass diffusion coefficient e > is a physical parameter and is 
regarded as being fixed. 

(Pi) Find PL ■■ {x,t) enx [0, T] ^ pl{x, t) G M, ul : [x, t) eUx [0, T] ^ ul{x, t) e and 
PL ■■ e O X (0,T] ^pL{x,t) e M such that 



(2.2a) 
(2.2b) 
d{pL ul) 

dt 
(2.2c) 

(2.2d) 
(2.2e) 
(2.2f) 



dpL 
dt 



{ul Pl) 







Pl{x,Q) = po(a;) 

• {pL UL ® Ul) - ■ (Kpl) D{ul)) + VxPL ^ PL f + ^x - t[iIjl) 

Vx ■ Ul ^ 
{plUl){x,0) = {poUo){x) 



infix (0,T], 



in n X (0,T], 
in SI X (0,T], 

on dn X (o,r] 



where V'l : {x,q,t) E n x D x [Q,T] ■iljL{x,q,t) e M, and t('0l) : {x,t) E n x {Q,T] 
i(V'l)(SiO ^ IR''^'' is the symmetric extra-stress tensor defined as 



(2.3) t{iPl) k 

Here /c e M>o, / is the unit dx d tensor, 

(2.4a) C^{^L){x,t) := / V'l(x, g, t) C/;(i Ig.p) ^q, 

(2.4b) f?(^L)(s,t) / VL(s,g,t)dg. 



and 



JD 

The Fokker-Planck equation with microscopic cut-off satisfied by -0^ is: 

K 



dipL 
dt 



(ML-y^)V'L + ^V5, 



a{uL)qrMC{pL)P 



(2.5) 



( i'L \ 

KCiPL) 



K K 



^ 1=1 7 = 1 



VC(pl)m 



in f7 X D X (0,T]. 



Here, for a given L > 1, e C(M) is defined by (1.22), a{v) = u, and 
(2.6) A G M^^^ is symmetric positive definite with smallest eigenvalue oq G 
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We impose the following boundary and initial conditions 

K 



1 ^ 



4A 



IpL 



(2.7a) 
(2.7b) 



n = 



on X X (o,r], i = i,.. 

ondflx D X (0,T], 



(c(>l)) 

(2.7c) ^i(., 0) = M(.) C(po(-)) /3'^(V'o(-, ^/(((Pol-)) ^'^(0)) > onnxD, 

where n is the unit outward normal to dfl. The boundary conditions for V'l on dfl x D x {0,T] 
and O X X (0, T] have been chosen so as to ensure that 



(2.8) 

inxD 

Henceforth, we shall write 



/ ipL{x,q,t)dqdx ^ / 2pL{x,q,0) dq dx Vte(0,T]. 
JnxD ~ ~ ~ ~ JnxD ~ ~ ~ ~ 



C{pl)m' 



C(po)M- 



Thus, for example, (2.7c) in terms of this compact notation becomes: tpL{', 0) — f3^{ipo(-, •)) on 

nx D. 

The notation | • | will be used to signify one of the following. When applied to a real number x, 
\x\ will denote the absolute value of the number x; when applied to a vector v, \v\ will stand for 
the Euclidean norm of the vector v; and, when applied to a square matrix A, \A\ will signify the 
Frobenius norm, [tr(A"'" A)] 5 , of the matrix A, where, for a square matrix B, tx{B) denotes the 
trace of B. 



3. Existence of a solution to the discrete-in-time problem 



Let 



(3.1) H -.^ {w e L^{n) :Yx -w ^0} and V := {w e Hl{n) : W x ■ w ^ 0}, 

where the divergence operator Vj;- is to be understood in the sense of distributions on fi. Let V' 
be the dual of V. 

For later purposes, we recall the following well-known Gagliardo-Nirenberg inequality. Let 
r E [2,oo) if d = 2, and re [2,6] if d = 3 and 9 = ^ (5 ~ Then, there is a constant 
C = C{n, r, d), such that, for aU e H^{n): 

(3.2) II^IIl^(O) <C||77|| 

Let T G C(M>o) be defined by J"(s) := s(logs - 1) -I- 1, s > 0. As lims_^o+ -^(s) = 1, the 
function J- can be considered to be defined and continuous on [0,cx)), where it is a nonnegative, 
strictly convex function with J-{1) — 0. We assume the following: 

mill ; Pmax 

], with pinin > 0; uq € H; ipo > a.e. on x D 



with J'(V'o) e L\j{n X D) and / il^n{x,q) dq € L°°(0); 

Jd ~ ~ ~ 

the Rouse matrix A G M^^-^ satisfies K^; 

], [Mmin,/^max]), C ^ {[P: 



nxD 



ijjo{x,q) dqdx = 1; 



-'mill 1 Pmax. 



•) [Ciiiini CmaxDi with fimini Cmin O7 



(3.3) f e L^iO,T:L^{n)) and A-S(0,5,f), 7, > 1, 



1,...,K in ( 1.5a b); 



where >r > 1 if d 2 and >r = | if d 3. For this range of x, we note from (3.2) that there exists 
a constant E M>o such that 



(3.4) 



wi ■ W2 dx 
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In (3.31, L^m(^ ^ P ^ [liOo): denotes the Maxwellian-weighted LP space over Q x D 

with norm 

Similarly, wc introduce L^,j{D), the Maxwellian-weighted space over D. Letting 



(3.5) 11^11^1,(^2x15) := (/ M 



M 

2 



dq dx 



we then set 

(3.6) X = Hl{nxD) := e Ll,{n x D) : Mhi^,,,d) < ^} 

It is shown in Appendix C of [TT] (with the set X denoted by X there) that 



(3.7) C°°(17 X D) is dense in X. 
We have from Sobolev embedding that 

(3.8) H\n,LliD))^L^in;Ll{D)), 



where s € [1, oo) if d = 2 or s e [1, 6] if d = 3. Similarly to (3.2) we have, with r and 9 as there, 
that there is a constant C, depending only on il, r and d, such that 

(3.9) MLrin:MlAD))<CMl^^^^^^^^^ ^ H\n-Ll,{D)). 

In addition, we note that the embeddings 

(3.10a) H\,{D) ^ LIj{D), 

(3.10b) Hlf{n xD) = L^{Q;Hli{D))r\H\n;Lli{D)) ^ Llj{n x D) = {VL- LIj{D)) 



are compact if 7^ > 1, i = 1, . . . , K , in ( 1.5a b); see Appendix D of ,llj . 

We recall the Aubin-Lions-Simon compactness theorem, see, e.g., Simon [5TJ Theorem 5]. Let 
60, B and Bi be Banach spaces, i3i, i = 0, 1, reflexive, with a compact embedding Bq ^ B and a 
continuous embedding B ^ Bi. Then, any bounded closed subset E of L'^{0,T;Bo), such that 

(3.11) / \\r]{t) - r]{t ~ e)\\l^ dt as 61 ^ 0, uniformly for r/ € 

Jb 

is compact in L^(0,T;B). 

We shall also require the following result (cf. Lemma 1.1, Chapter III, ^1, in |53)). 

Lemma 3.1. Suppose that B is a Banach space with dual space B' , {a,b) is a bounded open 
subinterval ofM., and u and g are two functions in L^{a, b; B). Then, the following three conditions 
are equivalent: 

(i) u is almost everywhere equal to the primitive of g; i.e., 

'"(^) ^ ^ + / di^) ^ B, for a.e. t e [a, b] ; 

J a 

(ii) For each test function Lp G C^{a, b), 

u{t)^{t)dt = - g{t)v{t)dt; 

(iii) For each rj Cz B' , 

-^(m,??) = {g,v) 

in the sense of scalar-valued distributions on (a, b), where (•, •) is the duality pairing between 
B and B'. 

If any of the conditions (i)-(iii) holds, then u is almost everywhere equal to a certain continuous 
function from [a, b] into B. 
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Corollary 3.1. Suppose that (a, 6) is a bounded open subinterval of R and let B be a Banach 
space. Suppose further that u € W^'^{a,b;B); then, 



u{t) — u{a) 



ds 



(s) ds for all t G [a, h] . 



Proof. Thanks to our assumption that u e W^'^{a,b]B), and definmg g = du/At, we have that 
u,g G L^{a,b]B)\ thus by the implication (ii) => (i) in Lemma 3.1 we deduce that 



(3.12) 



40 + 



(s) ds, S. ^ B, for a.e. t £ [a,b]. 



Since, again by Lemma 3.1 u is a continuous function from [a, b] to B, it follows that (3.12 1 holds 
for all t G [a, b]. Hence, by taking t — a in (3.12), we get that ^ = u{a). □ 

The next result is stated without proof as Corollary 1, Ch. XVIII, p. 470 in the book of Dautray 
& Lions [13. 

Corollary 3.2. Suppose that g G L^{a,b;B) and rj G B' . Then, 

5(t)dt,77^ =^ {g{t),v)dt, 

where (•, •) is the duality pairing between B and B' . 
Proof. We define the S-valued distribution u by 



u{t) 



g{s)ds, for a.e. t G [a,b]. 



It then follows that u,g £ L^{o,, b', B) and condition (i) of Lemma 3.1 holds with ^ = 0. Hence, by 
Lemma 3.1 condition (iii) of Lemma |3 . 1 1 also holds, which, after integration over [a, 6], then yields 

fb 



(3.13) 



{u{b),r]) - (u(a),?7) 



{g{s),r]) ds. 



The definition of the duality pairing (•, •) implies that \{g{s),r])\ < ||5(s)||e||?7||e' for all s G [a, 6] 
and all r/ G B'; and, since g G L^{a,b;B), the definition of the Bochner integral implies that 
la \\9i^)\\B ds < oo. It thus follows that right-hand side of (3.13) is finite for all rj G B'. By noting 
that u{b) — gis) ds and u{a) — 0, the identity (3.131 becomes 



5(s)ds,77^ = y (3(s),f7>ds, 
which, after renaming the dummy integration variable s into t, yields the required result. 



Throughout we will assume that (3.3) hold, so that (1.6) and (3.10a b) hold. We note for future 



reference that (2.4a) and (1.6) yield that, for Lp G ^ 



\C,{MV)\^ dx= f f M^U^q.qJdq 
A/(C/;)2|g,|4 dq 



dx 



(3.14) 



< 



< C 



D 



M\Lp\'^dq dx 



nxD 



M\ip\^dq dx 



nxD 



I = 1, 



.,K, 



where C is a positive constant. 

We state a simple integration-by-parts formula. 

Lemma 3.2. Let (p G H\j{D) and suppose that B G W^^"^ 
then, 

K ^ K 

(3.15) 



is a square matrix such that ix{B) = 0; 
/ M V(i?g,)-y,.¥'dg- / M ^Y^^M^^?) I^^l ■ B Aq. 
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Proof. See the proof of Lemma 3.1 in [S]. 



□ 



We now formulate our discrete-in-time approximation of problem (Pl) for fixed parameters 
e e (0, 1] and L > I. For any T > and N > 1, let N At ^ T and t^ = n At, n ^ 0, . . . , N . To 
prove existence of a solution under minimal smoothness requirements on the initial datum uq Cz H 
(recall (3.3)), we assign to it the function m*^ — tf{At) E V, defined as the unique solution of 



(3.16) 

Hence, 

(3.17) 



pou 



AtWxU° : V^v 



Ax 



PqUq ■ V dx 



e V. 



f[po + At |y, ]dx< f po lyoP dx < C. 
Jn Jn 



In addition, we have that po (u" — uq) 'U dx converges to for allv E H in the limit of At — > 0+. 

Analogously to defining ifi E V for a given initial velocity field uq E H, we shall assign a certain 
'smoothed' initial datum, 

= V'°(i,At) E H\j{n X £1), 
to the given initial datum -00 — '0o/(C(po) M) such that 



M 



C(/9o) V'V + At (V^V° • + VgiP" ■ Yiv)] dq dx 

MC(po)/3^(V^o)¥'dgdx Wip E x D) 



nxD 



nxD 

(3.18) 

For p E [1, oo), let 

(3.19) Zp := ^^(p E Lliin X D) : If > a.c. on n X D and J M{q) ip{x, q) dq E L°° (n) 

It is proved in the special case of C(Po(0) = 1 in the Appendix of [T^] that there exists a unique 
V'° e Hli{n X D) satisfying ( |3.18[ ); furthermore, E Z2, 

MC(po)J-(^°)dgdx 



nxD 



-I- 4 At 



M 



nxD 



dg da; 



< 



M C(po) -^(^0) dg dx. 



nxD 



ess.supa^gt^ / M ip^ dq < esssup^en Mip^dq 



D 



(3.20a) 

(3.20b) 
and 

(3.20c) = /3^(^") ^- ^"0 weakly in L\[{VL x D) as L ^ 00, At -^■ 0+. 

In the case of variable C(Po(')) the same properties hold under the assumptions on po and C stated 
in (3.3|. For example, the claim in (3.20c) that ip'^ — /3^(?/'*') a.e. on O x D follows from (3.18) on 
replacing ijp by ■0° — L on the left-hand side of (3.18) and /3^('0o) by (i^{4'o) — L on the right-hand 
side, which preserves the equality. We then take = [^^ — L\+ and note that l3^{tpo) — L < 
to deduce, thanks to the positivity of C(Po) on fi, that — L]+ = a.e. on x D, which then 
implies that < L a.e. on f2 x D. An analogous argument shows that -ijp >Q a.e. onVlx D. In 
particular, -0" = /3^(0°) a.e. ox\Q, x D and ij)'^ E Z2, as was claimed in the line above (3.20a I. In 
fact, V'" E L°^{nx D)f^ H\;{n x D). 

The proof of (3.20b I is based on a similar cut-off argument. By defining, for ijp = 0°(L, At), 
the function A^^ by 

A°(x) := / M{q)i^\x,q)dq, x E 
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applying (3.18) with ip{x, q) — ip{x) <Si !((?) and recalling Fubini's theorem, we have that 



[C(po) xl^ + At v,xl ■ y,(^] dx = / c(po) ^ 

and therefore, for each cj e M and all ip £ _ff^(f2), 

/ [C(po) (A^ - ^+ AtV^iXl - w) • dx 



Af(g)/3^(V'o)dg 



dx y^£H\n), 



Af(g)/3^(^o)dg 



dx. 



(3.21) 

Now thanks to ( |3.3| we have that 

0<C(po)^M(g)/3'^(V^o)dg<C(po) / M(g)V^odq= / Ms,<l)dq e L°°{n), 
and therefore by selecting 



(3.22) uj := ess.sup^go / M iq)ipoix, q) dq = ess.sup^g^ 

we get that 



D 



1 



CiPoix)) Jd 



^o{x,q)dq 



C(Po 



A/(q)/3^(^o)dq 



D 



< 0. 



Thus, by choosing ip — [A^^ — w]+ in ( 3.21| ), we deduce that 



dx < 0. 



Hence, [A^ — = a.e. on fi. In other words, < X^{x) < w a.e. on $7, which then implies 
(3.20b). We shall denote the mapping i/jq i-)- by Saix', i-e., "0° = SAt,L'4'o- 
Let us define 

(3.23) T ■.= {rjeL^irt):Tj£ [p^^^, p^,^] a.e. on fl}. 
It follows for all V, w £ H'^i^) that 

(3.24) w (g) -y : w = [(w • \i'x)w] ■ v — —[{v ■ \7x)v] ■ w + (v ■ \i'x){v ■ w) a.e. in i}. 

Noting the above, our discrete-in-time approximation of (Pl) is then defined as follows. 

(Pf-*) Let N e N>i and define At := T/N; let, further, p^ ~ po e T, ul ~ m° £ If and 
^0 — ^° £ Z2. For n = 1, . . . , N, and given u^"^ G ^nd 

(3.25) eL°°(t„_i,i„;L°°(17))nC([i„_i,t„];i2(r!))n^^i^°°(t„_i,i„;M^^ 



where q € (2, oo) when d = 2 and g G [3,6] when d = 3, such that p^*^|[t„„i,t„](-,tn-i) = p2 
and 



n-l 



dt 



dt- 



Pl 



Vajjydxdi = 



(3.26a) 

where the symbol (•, 



V77GLi(t„_i,t„;W^i'^(f7)), 
denotes the duality pairing between W^'^ (fl)' and W^''^ (V,), 



w • 1-1 (n) 

with q £ (2,oo) when d = 2 and g G [3,6] when d = 3. We then define PL^*'l[t„-i,t„]('i^n) ^ 

T. Given (p2"^ V-^"^) G T x F x Zj, find 
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such that 



/ 



At 



1 Pl-Pl 

2 At 



(3.26b) 



wdx+ fi{pl)D{ul) : D{w)dx 
Jn ~ ~ ~ Jsi ~ ~ ~ ~ ^ ^ 



da; 



At 



M — 



E 



1 ^ ^ 

— A, M ijl-[a{ul)q,] M Cipl) P'^ m 



"S/ q.Lpdq dx 



(3.26c) +/ eMV,^i)l- V,^Lpdqdx = Q S X; 



where, for t G [t„_i, i„) and n — \, . . . ,N , 
(3.27) =/"(•) 



1 

At 



f{;t)dt€L''{n). 



It follows froni (3.31 and (3.27) that 

(3.28) f^^^+ / strongly in L^{Q, T; i^(r2)) as At ^ 0+, 



where k > 1 \i d = 2 and >f = | if d = 3. Note that as the test function w in (|3.26b|) is chosen to 



be divergence-free, the term containing the density g of the polymer chains in the definition of r 



(cf. (2.3|) is eliminated from (3.26b I 



For n € {!,... ,-/V}, and for the functions & V and p^ G T fixed, the existence of 



unique solution 
(3.29) 



pf\^_„t„]&L°°{t^_utn■,L°°{n))r^c{[t^-M■,L'm) 



to (3.26a) satisfying the initial condition P^*' |[t„_i,t„] ('j ^n-i) = follows from Corollaries II. 1 
and II. 2 and the discussion on p. 546 in DiPerna & Lions [20] (with our G V here extended 
from r2 to M.'^ by 0). We refer to Appendix A for a justification that the notion of solution used 



in (3.25), (3.26a) is equivalent to the notion of distributional solution, used by DiPerna & Lions 



in ^20j. The statement 



pr'l[t„-.,t„] G w'-°"it„^^,t„;w'^^{ny) 



in (3.26a), with q € (2, oo) when d = 2 and q e [3, 6] when d = 3, follows from the bound 



r f ^r'pf'^ ■YxV'dxdt 



and the fact that ^p^^'' + V^; • (y^ ^Pl^*') = in the sense of distributions on x (t„-i,tn); 



hence (3.25) 
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A further relevant remark in connection with (3.26b I is that on noting that — 



and p\ ^ — Pl^*' ('j ^ri-i), the second term on its left-hand side can be rewritten as 



(3.30) 



Jn 



At 



ul^ ■ wdxdt = — — — 
~L ~ ~ 2At 



dt 



dt 



2At 



t„_i Jn 



At 



pf*^ul~^ ■Yx{ul-w)dxdt 



pf'^ dt ) [ui"^ ■ Yx{ul ■ w)] dx yw G H\n), 



where q <E (2,oo) w hen d — 2 and q S [3,6] when d — 3. The first equality in (3.30) is a conse- 
quence of Corollary 3.1 and Corollary 3.2 on noting that ^p'^ G L°° {tn~i,tn;W {i})') C 
L^(t„_i,t„; with q g (2,cx)) when d = 2 and g G [3,6] when d = 3. 



The identities (3.24) and (3.30) now motivate the form of the expression in the second line of 



(3.26b). We note here that the requirements that q > 2 when d = 2 and q> 3 when d = 3 are the 



consequence of our demand that the scalar product -w oi the functions u^jW G (fl) belongs 
to W^-^{rt), which is required in ( |3.30 ). 



As p^' e T 



e V (extended from 17 C M"* to R'^ by 0), ' G C{[tn^i,t„]; L^{n)), 
C G C'^([/r?min, Pmax], [Cmin , Cmax] ) , it foUows from CoroUary II. 2 in the paper of DiPerna & Lions 
[20] that C(Pl^*') is a renormalized solution in the sense that 



'5c(pr') 



^ di 

t 



(3.31) 



where q G (2, oo) when d = 2 and g G [3, 6] when d = 3. Hence, on observing that C(Pl) 
apT\-^tn)) and C(pr') = C(pr'(-,<n-i)), we have that 



(3.32) 



/ 

Jn 



api) - cipr') 



ipdxdt 



At 



1 

'At 



dt 



ip) dt 



At 



At 



t„-i \ 

c(pr')y2"'-yx^dxdi 

' ap^t^)dt] [ui'^ ■ dx yip G w''^{n), 



where q G (2, oo) when d ~ 2 an d qr G [3, 6] when d = 3. The first equality in ( 3.32 ) is a consequence 



of Corollary 



3.1 



and Corollary 



3.2 



q G (2, oo) when d = 2 and g G [3,6 



on noting that ^C(pL^*') e Li(t„_i, i„; Vt^^^^ (rj)'), with 



when d = 3. Since G X, it follows from (3.32) that 



1 / c(p£)-c(pr%~„^,,,,. 



yip G X, a.e. q e D, 
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and therefore we can rewrite (3.26c) in the fohowing equivalent form 

'cipi) n - c(p2"') v^r' 1 c(p'i) - c(pr') 



M 



nxD 



At 



At 



If dq dx 



M 



flxD 



At 



K 

,E 



dq dx 

\/ q.Lpdq dx 



(3.33) 



eMV^il^l- V^tpdqdx = Q ^ip e X. 



The following elementary result will play a crucial role in our proofs. 

Lemma 3.3. Suppose that S CZ R is an open interval and let F G W\oci^)- further G denote 
the primitive function of s £ S* i— >■ sF"{s) G M; i.e., G'(s) = sF"{s), s Cz S . Then, the following 
statements hold. 

a) s G I— 7- sF'{s) — F{s) — G{s) G M is a constant function on S; i.e., there exists Cq G M 
such that sF'(s) — F(s) — G{s) = cq for all s € S. 

b) The following identity holds for any a,b G S and any A, B Cz M; 

{Aa-Bb)F'{a)-{A-B)G{a) ^ A{F{a)+ c^)- B{F{b) + cq) + B{h- af [ F" {9a+ {l-0)b) 9 d9. 



c) // in addition B > and there exists a do (zR such that ess.infgg[o,i]-P"'(^'^+ (1 ^^)^) do, 
then 

{Aa - Bb)F'[a) -[A- B)G{a) > A{F{a) + cq) - B{F{b) + co) + ^doB{b ~ a)^ . 

Proof. Before we embark on the proof of the lemma, we note that G{s), the primitive function of 
s G 5 sF"{s) G M, is only defined up to an additive constant. Altering the value of G{s) by 
an additive constant 7, say, does not affect the validity of the above statements: on replacing the 
constant cq by another constant, cq — 7, the statements a), b), c) continue to hold. 

a) By hypothesis the function s G S* H> sF'{s)-F{s)-G{s) G M belongs to W^io^ (5*), and it is 
therefore absolutely continuous on any compact subinterval of S; hence it is differentiable 
almost everywhere on 5. Upon differentiation and using that G'{s) — sF"{s), we deduce 
that the first derivative of s G 5 1— sF' {s)—F{s) — G{s) G M is equal to almost everywhere 
on S. This implies the existence of a constant cq such that sF'{s) — F{s) — G{s) = cq for 
a.e. s G S". Since the function s e S 1-^ sF'{s) — F{s) — G{s) G M is absolutely continuous 
on each compact subinterval of 5, it follows that sF'{s) — F{s) — G{s) = cq for all s G S. 

b) Suppose that a,b € S and A, _B G M. By Taylor series expansion with integral remainder, 



AF{a) - BF{b) = AF{a) - B 
Hence, 



F{a) + {b-a)F'{a) + {b-af / F" {Oa + [l - 9)b) 9 d9 



AF{a) - BF{h) + B{b - af [ F' 

Jo 



{9a + {I - 9)b)9 d9 



= {A- B){F{a) - aF'{a)) + {Aa - Bb)F'{a) 

= -{A - B){aF'{a) - F{a) - G{a)) -{A~ B)G{a) + {Aa - Bb)F'{a) 

= ~{A - B)co -{A- B)G{a) + {Aa - Bb)F'{a). 

By transferring the first term on the right-hand side of the last equality to the left-hand 
side, the stated identity follows, 
c) The inequality is a direct consequence of the identity stated in part b) of the lemma. 
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□ 



Remark 3.1. Lemma 3.3 should he compared with the discussion between eqs. (2.27) and (2.28) 
in reference |5j , which can be seen as a special case of Lemma 

In order to prove the existence of a solution to (P^*), we require the following convex regular- 
of F defined, for any 5 G (0, 1) and L > 1, by 



ization e C^'^ 
(3.34) 
Hence, 
(3.35a) 

(3.35b) 

We note that 
(3.36) 

7L1" 



His) 



^+s(log(5-l) + l 
F(s) =s(logs-l) + l 



2L 



s(logL- 1) + 1 



His) > 




for s < S, 
for S <s <L, 
for L < s. 

for s < 5, 
ioT 5 <s <L, 
for L < s. 



for s < S, 
for 6 < s < L, 
for L < s. 



C{L) 



for s < 0, 
for s > 0; 



and that [F^]"{s) is bounded below by 1/L for all s e M. Finally, we set 

(3.37) f3^(s) {[HTr^s) = max{/3^(s), <5}, 

and observe that Ps{s) is bounded above by L and bounded below by 5 for all s G 
that both (3^ and (3-^ are Lipschitz continuous on M, with Lipschitz constants equal to 1. 



Note also 



3.1. Existence of a solution to (P^*). On recalling the discussion following (3.29), for n. e 

{1, . . . ,iV} we define the function pi'l := PL^*'|[t„_i,t„]('7^n) G L°°(f7); it will be shown below that 

p1 £ T, in fact. With p^*' |[t„_i,t„] thus fixed (and with its values p^"^ ^nd p£ at i = i„_i and 
t = tn, respectively, also fixed,) we rewrite (3.26b) as 

(3-38) b{ul,w) = ^bi^Diw) e V- 

where, for all Wi G H\{Vl), i — 1,2, 




dx 



and, for all (p e i^/(f^ x D) and w e _ffj(f7). 



(3.39b) ib{ip){w) := 
It follows from Korn's inequality 



n—l n—1 

Pl 



Atplf^ 



Atk 



K 

E 

i=l 



Ax 



dx. 



(3.40) 



\D{w)\^dx>cn\\w\\]j^ 



where cq > 0, that, for u^^^ G V and p^^^ , p'l G T fixed, 6(-,-) is a continuous nonsymmetric 
coercive bilinear functional on H\{n) x H\{VL). In addition, for G and p'l^^ , p^l G T 

fixed, thanks to (3.4) and (3.14), tb{Lp){-) is a continuous linear functional on H}^{Vl) for any 
^&Ll,{nxD). 



24 



JOHN W. BARRETT AND ENDRE SULI 



It is also convenient to rewrite (3.26c) (or, equivalently, (3.33)) as 
(3.41) a(^£,^) = 4 (<,/3^(V^E ))(¥') V^eX, 

where, for all (pi ^ X , i — 1, 2, 



a(<y5l,</52) 




4A 



dq dx 



M 



flxD 



i=l j = l 



(3.42a) , 

i=l 3 = 1 

and, for all v G i/H^^), ^ L°^{n x L>) and e ^, 



dq dx, 



(3.42b) 4(i;,??)(^) := 



M 



C(p2"')V'rV + At^[a(z;)g.]C(p2)ry • V,.^ 

4=1 ~ ~ ~ ~ . 

JAt] 



dg dx. 



Hence, on noting (2.6), for G and fixed (and therefore p£ and alsofixed), 

a(-, •) is a coercive bilinear functional on X x In order to show that a(-, •) is a continuous bilinear 
functional on X x X, we shall focus our attention on the case of d = 3; in the case of d = 2 the 
argument is completely analogous; we begin by noting that, by Holder's inequality and the Sobolev 
embedding theorem, 



^^IKl ^\Wi\\'^xV2\dqdx < I MWyJl ^\\L6(n)\Wi\\L3(n)\\Y xP2\\L'^(n)<^q 

~ J D 

< c(fi) ||M2"^||i6(o) / M\\tpi\\Hi{n)\\'i32\\m{n}(iq 



D 

< C(f^) ||w2~^IU'5(0) ll<y5l||L|,(D;Hl(f^)) ll'^2||L|j(D;Hl(n)) 

< \\ul'^\\LB{n) \\y^l\\Hlj(nxD) \\V2\\Hlj(nxD)- 

This, obvious applications of the Cauchy-Schwarz inequality, and the fact that the range of the 
function is the compact subinterval [Cmin, Cmax] of (0, oo), then imply that a(-, •) is a continuous 
bilinear functional on X x X. 

In addition, for all v £ H^{il), -q e L°°{Q, x D) and e X, we have that 



\ia{v,rj){^)\<\\Q{pr')rL-' 



(3.43) 



M\q\^dq 



Llj(nxD) \\f\\Llj{nxD) 

C(Pl)'7||l-(OxD) l|VJ;^;||L2(o) W^^qVllLl^inxD)- 



Therefore, by noting that C{pI~^) i'T'^ e Z2, C(p2) V G L°°{n x D) and recalling (|L3), it follows 
that ^a(u, rj){-) is a continuous linear functional on X for any v € |^^(^^) and rj G L°°(il x £)). 

Before we prove existence of a solution to the problem (3.26b c), i.e., (3.38) and ( |3.41[ ), let 
us first show by induction that the function P^*' |[t„_i.t„], whose existence and uniqueness in the 
function space 

L°°(i„_i,i„;L°°(r!)) n C([t„_i,t„];i2(r!)) n l¥i^°°(t„_i,i„; W^i'i^(f})') 
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has already been established, with q E (2,oo) when d ~ 2 and q G [3,6] when d — 3, satisfies the 
two-sided bound pmin < P'L^\x,t) < Pmax for a.e. x € and every t e [i„_i,f„], i.e., that e T 
for all n G {0, . . . , N}. To this end, for a € (0, 1), we consider the regularized problem 

^ ( ^ ' - Jo ■ 
(3.44a) +af" / V,pf • V,7/da;dt = Vr; e i2(t„_i, i„; 

subject to the initial condition 

(3.44b) pt%t„-^.a;tn~i)=Pr'eT, 

where p^^^ G T was assumed for the purposes of our inductive argument; clearly, := po G 
so the basis of t he indu ction is satisfied. We begin by showing the existence and uniqueness of a 
solution p^^' to (3.44a b) and that p,nin < Pl^q (s, t) < Pmax for a.e. x G and for all t G [tn-i,tn]; 
we shall then pass to the limit a — ^ 0+ to deduce that the limiting function, which we shall show 
to coincide with p^*', satisfies the two-sided bound p„iin < p^*'(.T,<) < pmax for a.e. x E and 

for all t e hence in particular we shall deduce that = Pi^'k'i^n) G Tf. 

The existence of a unique weak solution 

to ( |3.44a b) is immediate; see, for example, Wloka [SS], Thm. 26.1. Further, on selecting, for 
s e {tn-i,tn], the test function rj = X[i„_i,s] p'l^o equation (3.44a|, where for a set 5* C M, xs 
denotes the characteristic function of S, and noting that u^~^ S V, we obtain the energy identity 

\\p^L'M\hin) + '2^ r I iy.pL^*i(^)pdxdt = iip2-^iii.(a), 5e(t„-i,t„], 

which then implies that {p^^^ }ae(o,i) is a bounded set in the function space L°°(i„_i, i:„; L^(ri)) 
and that {-v/ay^PL^Q }Qe(o,i) is a bounded set in L^(t„_i, t„; L'^(J7)). It then follows that 

< — > is a bounded set in L^(t„_i, i„; i7^(il)'). 

I J Qe(o,i) 

Hence there exists an element p^g e L°°(i„_i, t„; L^(f2)) n i/^(t„_i, i„; iJ^(ri)') and a subse- 
quence of {Pi^a }ae(o,i) (^of indicated) such that, as a — > 0+, 
(3.45a) pl^"*! ^ p^""*! wea^ in L°°(t„_i, i„; ^^(r!)), 

(3.45b) «V,p[^'^*l ^ strongly in L^{t^^^,tr,; L^ifl)), 

(3.45c) |p^^*i -> |p^^*l weakly in L^{t^^i,t^; H\nY). 

By a weak parabolic maximum principle based on a cut-off argument, we also have that 

and therefore. 



Pmin < P^a < Pmax, VL > 1, Va G (0, 1), 



Pmin < P^o ^ /'max, VL > 1, Vck G (0, 1). 

Thus, on passing to the limit in (3.44a I we deduce that 
(3.46) V77eL2(t„_i,t„;i?i (!])). 
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As 



Jt„_i Jn 



is a continuous linear functional for all q € (2, oo) when d — 2 and all q G [3, 6] when d — S, the 
application of a density argument to (3.46) yields that 



A dt - 



dt 



W '9-1 (O) 



where q € (2,oo) when d — 2 and q G [3,6] when d ~ 3, and (-jin-i) = ^. As p^'' is 
already known to be the unique weak solution to this problem by the argument from the beginning 
of this section, it follows that = Pl^*') ^nd therefore 

(3.47) pmi„ < P^*'|[t„_i,t„](s,0 < Pmax for a.c. X eil, for all t e [i„-i,t„] and n = 1, . . . , A^. 

In particular, for t = t„, p^in < p'2(x) P^*' |[t„_i,t„] (j, tn) < Pmax for a.e. x € f^; hence, 
p2 e T, as was claimed in the first sentence of this section. 

Thus, for any given p2~^ £ T and w^"^ ^ F' ^ ^ {1, . . . , A^}, we have shown the existence of a 
unique function 

such that p^*'|[t^_j t^](-,i„_i) = p2 with p^ := po e T when n = 1, and 



dt 



dt 



w • 9-1 (O) 



p^"^*! u^"^ • V:rr?dxdt = 



(3.48a) 



V77eLi(i„„i,t„;W^i'^(ri)), 



[At], 



e T. 



where q G (2, oo) when d — 2 and g G [3, 6] when d = 3, and p^ ' l[t„_i,t„] 

We now fix p£(-) := Pl^*' l[t„_i,t„] ('j ^n) S T, and we turn our attention to the proof of existence 
of solu tions to (3.26b c). To this end we consider the following regularized version of the system 
( |3.26b[ c): for a given 6 G (0, 1), find (u^.a^^E,*) e F x ^ such that 

(3.48b) b{ul,, w) = Wls){w) Vu; G V, 

(3.48c) aii^ls, = Uuts.Psii^lMv) ^ ^- 



We emphasize at this point that (3.48a) decouples from (3.48b c); indeed, given p^ G T and 



,n-l 



G V , one can solve (3.48a) uniquely for 

where q G (2, oo) when d = 2 and q G [3,6] when d — Z, and p^^ *^ (•; ^k- i) = P^^{')] by 
defining p2(-) '■— PL^*'l[t„_i,t„]('7*n): can then consider the system (3.48b c) for {^£^,^'£5} 
independently of (3. 



The existence of a solution to (3.48b c) will be proved by using a fixed-point argument. Given 
■0 G L\[{n X D), let (u*, ■0*) G y X A be such that 



(3.49a) 
(3.49b) 



b{u* , w) — £},{ip){w) \/w G V, 



The Lax-Milgram theorem yields the existence of a unique solution u* G to ( 3.49a ) for a given 
■0 G X, and the existence of a unique solution tp* e X to { 3.49b ) for a given pair (u*, ^) x X. 
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Therefore the overall procedure (3.49a b) that maps a function ip £ L'j^{il x D) into tp* E X is 
well defined. 

Lemma 3.4. Let T : L\j(^ x D) ^ X C L \f{il x D) denote the nonlinear map that takes the 
function to tp* = T{ip) via the procedure ( 3.49a b). Then , the mapping T has a fixed point. 
Hence, there exists a solution {u1 g, V'E s) ^ Y ^ ^ (3.48b c). 

Proof. The proof is a simple adaption of the proof of Lemma 3.2 in [S], where — p^^^ = 1 
and C(') is identically equal to a positive constant. Clearly, a fixed point of T yields a solution 
of (3.48b c). In order to show that T has a fixed point, we apply Schauder's fixed-point theorem; 
that is, we need to show that: (i) T : Llj{rt x I?) -> L\.i{il x D) is continuous; (ii) T is compact; 
and (iii) there exists a C* G M>o such that 

(3.50) \mLl,in>cD) < a 

for every ip e L\i{n x D) and k e (0, 1] satisfying V' — kT{iP). 
(i) Let {^/'^^•'}p>o be such that 

^p^P^ 2p strongly in L^/(ri x D) 



(3.51) 



as p — ^ oo. 



It follows immediately from (3.37) and (3.141 that 

(3.52a) M5 /3|'(V'^^^) ^ /3|'(V^) strongly in L''{n x D) as p ^ oo, 

for all r E [1, oo) and, for i = 1, . . . , i^T, 

(3.52b) C^{MC(pl)i^^P'>) ^ C,{MC{pl)i^) strongly in as p -> oo. 

In order to prove that T : i^/(^ x D) — > -^M(i^ x D) is continuous, we need to show that 
(3.53) ^(P) :=r(?A^P)) ^r(V') strongly in X £>) asp^oo. 

We have from the definition of T (see ( 3.49a[ b)) that, for all p > 0, 
(3.54a) ai7i^P\^) = 4(«(f), /3,^(^(f)))(^) G X, 

and w*^^^ G V satisfies 

(3.54b) b{v 



Choosing ip = ry^P^ in (3.54a) yields, on noting the simple identity 
(3.55) 2 (si - S2) si = Si + (si - §2)^ - S2 Vsi,S2e 



and (|3.3|), (|3.37j) that, for aU p > 
M 



L I 

oq At 
2A 



2eAi |Va;ry 



<P)|2 



d(? dx 



(3.56) 



< 



Mc{pi-')m 



nxD 



'dqdx + C{L,X,Cn...,ao')At [ |V,«(ff da;. 



Choosing w = v^p'> in ( |3.54b| , and noting ( |3.55| , ( |3.3| , ( |3.4| , ( |3.14| , ( |3.40[ ) and ( |3.5l| yields, for 
all p > 0, that 



(3.57) 



Pl I" 



< 



(P)|2 



dx + pnun At 



n-1 L,n-l|2 



Pl 



dx + cAt\\r\\i 



CAt 



D{v^P^)\^dx 

M|7/^(p'pd(7 dx < C{L). 



Combining ( |3.56[ ) and (3.57), noting that by (3.3) the range of ( is the compact subinterval 
[Cmin , Cmax] of (0,oo), and recalling (3.40), we have for all p > that 

(3.58) WTj^P^Wx + ||«('')|ki(o) < cm (At)-i) . 
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It follows from (3.581, (3.8) and the compactness of the embedding (3.10b) that there exists a 
subsequence {(TytP*)^ U^^''^)}pk>o and functions ry e X and w G such that, as pk — oo. 



(3.59a) 



weakly in {n; LIj {D)) , 



(3.59b) 


M5 V^Ti'-P"^ 


M3 V^r] 


weakly in L'^{fl x D), 


(3.59c) 






weakly in L^{fl x D), 


(3.59d) 




—7- ry 


strongly in L\.j{n x D), 



(3.59e) 



,ipk) 



V 



weakly in H^{n); 



where s e [1, cx)) if d = 2 or s e [l,6]Jf d = 3. We deduce from ( [slib] ), (|339a)b), ^l9e^ and 
(3.52b) that the functions v^V and if) ^ X satisfy 



(3.60) 



6(u,u;) =4(V')(w) Vu;e\f. 



It follows from (3.54a), (3.42a b), (3.59a-e) and (3.52a) that rj, ip G X and v G V satisfy 
(3.61) ai?j,^) = ia{v,p^{i^)){ip) e c°°{n^). 

Then, noting that a(-, •) is a continuous bilinear functional on X x X, that daiu, f3g{ip)){-) ^ 
continuous linear functional o n X, and r ecalli ng (3.7), we deduce that (3.61 ) holds for all cp G X. 
Combining this X version of (3.61) and (3.60), we have that rj — T{ip) £ X. Therefore the whole 
sequence 

^(P) = r(^A^P)) ^ T{^) strongly in Llj{n x D), 

as p ^ oo, and so (i) holds. 

(ii) Since the embedding X ^ x D) is compact, we directly deduce that the mapping 
T : Llf{il X D) — >■ Lj^,j{il X D) is compact. It therefore remains to show that (iii) holds. 

(iii) Let us suppose that -0 — kT{iI'); then, the pair {v^-ip) gV x X satisfies 



(3.62a) 
(3.62b) 



b{v,'w) ~ £i,{^p){w) \fwGV, 



Choosing w ^ v in (3.62a) yields, similarly to (3.571, that 



n I |2 I n — 1 I n—l\2 n — 1 \ n—l\2 



(3.63) 



At 



dx + At / fi{pl)\D{v)\^dx 
f Pir-vdx-ky2 [ C,{MC{pl)i^):W^vdx . 



Selecting if = [J'i']'(V') in (|3.62b|), defining G W^^^iR) by 

^(s2+52)_l 



(3.64) 



2(5 " 
S - 1 



2L 



a S<S, 

ifsG[5,L], 
if s > L; 
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using that, thanks to ( |3.37| , [g^]'{s) = s/l3f{s) = s[T^]"{s); and that, by virtue of ( |3.32[ ) with 
[^I'KV'), we have 



/ X, / C(pL''*' ) V- • V, [T^n^P) dq dx 



M 



i{pf'^)dt\ul-^ ■ V,[gf](^)dgda; 



(3.65) 



inxD 



The convexity of and Lemma 3.3 with cq on noting that s[Tg]'{s) — J-g{s) — Gsi^) = 0, 
then imply that 



M 



JnxD 



At 
4A 



K K 



-sAtf MV^ij-V^{[F^]'{ij))dqdx 



SlxD 



(3.66) 



r 

<KAtV/ MC(p2)o-(^')9^ • Vg.V^dgda 
JnxD a ~ ^ ~ ^ - 

= KAtJ2 C,iMapl)^)--CTiv)dx, 
. -, Jo ~ ~ ~ 



where in the transition to the final inequality we applied (3.15) with B := g;(jj) (on account of 
it being independent of the variable q), together with the fact that tr(CT(w)) = V^; ■ v = 0, and 

recalled ( |2.4ap . 

Combining (3.63) and (3.66), and noting (3.35b), (3.3), (3.4) and (3.40) yields that 

'Pl\v\'+Pr'\^~<''\'] dx + nAt [ M)\Div)\'dx 



n 



k / Mapl)H{i^)dqdx 
JnxD ~ ~ 



kL-^ At 



M 



nxD 



dq dx 



<KAt [ pir.vdx + ^ f Pi 

Jn ~ ~ ~ Jn 



•"-I'^-^dx 



< 



+k I Mapr')Hi^^^r')dqdx 

K At fl 



nxD 



[ \D{v)fdx + 
Jo ~ ~ 



,2 fi2 

max II rn m 2 



\\r 



iL^in) 



(3.67) 



^ / p2~'K"Tdx + fc / MC{pl-^)T,^{Kijr^)dqdx. 
^ Jo ~ ~ JnxD ~ ~ 



It is easy to show that J^gis) is nonnegative for all s G M, with J^g{l) = 0. Furthermore, for 
any k e (0, 1], T^ins) < T^i-s) if s < or 1 < k s, and also Fgins) < J"i'(0) <lifO<Ks<l. 
Thus we deduce that 



(3.68) 



J'^ins) <J-^is) + l VseM, VKe(0,l]. 
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Hence, the bounds (3.67) and (3.68), on noting (3.36), give rise to the desired bound (3.50) with 

Cmin, /, Wl"^ and Therefore (iii) holds, and so 

□ 



C* dependent only on 6, L, k, 



mm ; pmax i 



T has a fixed point, proving existence of a solution to (3.48b c). 



Choosing w = yJl gin (3.48b) and = s) (3.48c|, and combining and noting (3.3), 

then yields, as in (3.67), with C{L) a positive constant, independent of 5 and At, 



I , n 1 2 I „n — 1 K , n 



Pl\Ul,s\ +Pl 
+ At 



'■L,S 



,n-l\2 



dx + k 



\ / Mp2)I^K,5)l'dx + fcL-ie / A/ iV.VE.^Pdgdx 
1 ~ ~ ~ JnxD ~ ~ ~ 

^^''"^ f M\w,rL.s\'didx 

nxD ~ ~ ~. 



+ 



4A 



< 



n-l I ,n-l|2 



Pl 



dx 



At rP' C"^ 

rmax^^f ||y'^|| 

2/^min Co 



^ II rn 1 1 2 



+fc / MC{pr^)T^wr^)dqdx 



nxD 



(3.69) 



< CiL). 



Equatio n (3.48a) being independent of S, we are now ready to pass to the limit (5 — )• 0+ 
in (3.48b c), to deduce the existence of a solution {(Pi^*'|[t„_i,t„],yi;,'0£)}n=i to (P^*), with 
Pl =pf*'(.,^n) e T, ul eV andi^le Xr]Z2,n = l,...,N. 

Lemma 3.5. There exists a subsequence (not indicated) of {{ul s,iJjl <5)}i5>0; and functions u£ G 
V and ^2 e X n Z2, n e {I, . . . , N}, such that, as 5 0+, 

(3.70a) u2,(5 ~^ ""2 weakly in V, 

(3.70b) < 5 -> < strongly in L''{n), 

where r G [1, 00) if d — 2 and re [1,6) if d = 3; and 



(3.71a) 






weakly in L^{il x D), 


(3.71b) 




^ Af ^ v,V^£ 


weakly in L^(il x D), 


(3.71c) 




^ Af ^ 


weakly in L^{n x D), 


(3.71d) 






strongly in L^(ri x D), 


(3.71e) 






strongly in L*(f2 x D), 


for all s G 


[1, 00); and, /or i — I, . . 


,K, 




(3.71f) 


c.{Mapi)rL,s) 


C^{M dpi) V-Z) strongly m L^{n) 



Furthermore, (p^^*' |[t„_i.t„], Ul, V'E) solves {3.26a~c) for n — 1,...,N. Hence, there exists a 



solution{[pf\^_^^t^-^,ul.rL)}n=i to (Pf ), withpl= p^'i-^tn) & T , ul & V and^l € XnZ2 
for all n — 1, ... , N . 



[At], 



Proof. The weak convergence results (3.70a), (3.71a) and that ip^ > a.e. on x _D follow 



), on noting ( 


3.40 


) and ( 


3.36 



values of r. 
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The results (3.71b c) follow from (3.69); see the proof of Lemma 3.3 in [HI for details. The 
strong convergence result (3.71d) for -!/)£ g follows directly from (3.71a-c) and (3.10b). In addition, 
( [3jTel f) follow fro m (|3.71d[ ) (|3 37| , ( [2^ an d p34|. 

It follows from ( |3.70alb) , ( |3.71b| -f)~ p!39a| b)^ l3^2a| b), ( [Os] ) and (|3J| that we may pass to 
the limit (5 0+ in (3.48b c) to obtain that (m^j-^E) ^ x X with V'E > a.e. on ft x D solves 
( |3.38[ ), ( |3.41| ); that is, ( |3.26b[ c). 

Next we shall show that 

(3.72) / M{q)rL{x,q)dqdxeL^{n), 

Jd ~ ~ ~ 

uniformly with respect to i > 1 for all n G {1, . . . ,N}. Hence we shall deduce in particular that 
V'E S ■^2- We begin by selecting (p{x,q) — (p{x) (X) l{q) in (3.26c) with [p e H^{n), which then 
yields that 



M 



At 



M 



nxD 



C(p[^*l)dt ul-\x)-(W.^)rLdqdx 



el Afy:rV'E ■ Yx^dqd^ = y^eH^iVL). 
JnxD ~ ~ ~ ~ 



Xlix) / M{q)rL{x,q)dq, n = 0, . . . , iV, 



We define 



with tpl :— — and note that G H^{Vl). By Fubini's theorem we then have that 



At 



C(p['^*i)dtU2-^-(y,(?)A2dx 



tn-l 



Jn 



This in particular implies, using the identity (3.32) with ip replaced by a;(/7, w G M, that, for each 
w G M, 



At 



(p{x) dx 



At 



C(p[^*i)dtUri.(y,^) {xi~Lo)dx 



Jn~ 

On selecting ip — [AJ ~ in this identity we then deduce that 

[A2 — uj]+dx 



At 



At 



C(p[^*i)dtUr^.y.([A2-u;]+f dx 



+ e / |y.([A2-a] + )pdx = 0, 
Jn 
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and hence, by omitting the (nonnegative) last term from the left-hand side, we have that 



At 



(3.73) 



At 



C(p[^*i) di ■ {[XI - uj]+f dx < 0. 



t„-l 



As, once again using the identity (3.32 1, with (p ([A£ — + ) this time, we have for each uj g 
that 



At 



At 



C(p[^*i) dt ui-' ■ y. {[XI ~ u]+f dx - 0, 



we can rewrite the second term in (3.73) to deduce that 

api){xi-u:)-c{pr'){xr' 



-UJ 



At 



[A2 — u!]+dx 



(3.74) 



1 MM)_^(£n,|,j_„,j.,,<„^ 



The inequality (3.74 1 can be restated in the following equivalent form: 

[Xl-u:]+dx 



^(^n-i)(A£-^)~(Ar 



At 

c{Pi)-c{pr') 

At 



1 f c{Pi)-c{pr' 



{XI -u) [Xl^L^Udx 

{[Xl-u;] + fdx<0, 



At 

and hence, after simplifying the sum of the second and the third term on the left-hand side, 



{Xl-cu)-{X", 



CO) 



api- r^'^--'-^-^ — ^ [XI - dx 



(3.75) 



1 f c{pi)-c{pr') 

At 



{[X'l-LuUfdx<0. 



Since s e M | ([s — uj]+)'^ G M>o is a convex function, we have that 



{Xl-u:)-{Xr' 
At 



-^[A2--].>^(([A2 



^]+) -([A' 



Using this inequality in the first term of (3.75), on noting that C(Pl ) ^ Cmin > and multiplying 
the resulting inequality by 2At, we get that 



api-'){i[K-^hr-m-'-^]+f}dx+ / iapi)~apr')) {[xi-^]+rdx<o, 

and therefore, upon simplification, 

C(Pl) UK + f dx< I C(pr') ([Ar' - ^] + f dx. n = 1, . . . , ^. 
Jn 

Thus, on recalling that by hypothesis p^, — Poi have that 

(3.76) < / C(p£) {[XI - dx< I C{po) {[Xl - u;]+y dx, n^l,...,N. 
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Now we choose lu as in (3.22), which yields [A° — (^]+ a.e. on il, and therefore, by (3.761, also 



a.e. on 51 for all n e {1, . . . , N}; in other words. 



(3.77) 



0< Kix) <u} := ess.sup^gi^ 



1 



C{po{x)) Jd 



V'ofe,g)dg 



for a.e. a; £ £7 and all n = 0, . . . , A^. Since uj here is independent of L, by recalling the definition 
of A2 we thus deduce that ^ ^2, uniformly with respect to L, as was claimed in the line below 
(|3J2l). 



t. 



Finally, as (p'liIi'li^'l) G T x y x Z2, performing the above existence proof at each time level 
, n = 1,...,7V, yields a solution {(pL""*' - ^2' ^ (Pf ) with pi 
= I, . . . , N, by noting that p^*' thus constructed is an element of C([0, T]; L^(ri)) 



□ 



4. Entropy estimates 



Next , we derive bounds on the solution of (P^*), independent of L. Our starting point is 
Lemma 3.5 concerning the existence of a solution to the problem (P^*). The model (P^*) includes 
'microscopic cut-off' in the drag term of the Fokker-Planck equation, where L > 1 is a (fixed, but 
otherwise arbitrary,) cut-off parameter. Our ultimate objective is to pass to the limits L — > cxd 
and At 0+ in the model (P^*), with L and At linked by the condition At — o{L^^), as 
L -> 00. To that end, we need to develop various bounds on sequences of weak solutions of (P^*) 
that are uniform in the cut-off parameter L and thus permit the extraction of weakly convergent 
subsequences, as L — >■ 00, through the use of a weak compactness argument. The derivation of 
such bounds, based on the use of the relative entropy associated with the Maxwellian M, is our 
main task in this section. 



Let us introduce the following definitions, in line with (3.27): 



(4.1a) 
(4.1b) 



At 



in ^ 



a1 ^(•), t e [tn-l, l^n 



tn], n=l,...,N, 



«f'+(•,^):=^."(•), ut''-i;t):- 



(•), t e (t„_i,i„], n = 1, 



We shall adopt y^*^'*^ as a collective symbol for u^*, m^*'^. The corresponding notations 

and p^**-'*'', and V'lS V'l*'^ ^nd ip^^'"'^'^ are defined analogously. In addition, we define 
the products (pl^^l)^*, (plWl)^*'± and (pl yL)^*^'^); and (C(pl)^l)^*, (C(pl) ^^l)^*'^ and 
(C(Pi) V'l)^*''^'' similarly. The notation p^* signifying the piecewise linear interpolant of pL with 
respect to the variable t is not to be confused with the notation p^*' , which denotes the function 
defined piecewise, over the union of time slabs £7 x n = 1, . . . , iV, as the unique solution 



of the equation (3.26a) subject to the initial condition p^*'(-,i„_i) = p2 ^ , n = \, . . . ,N , with 



p'^ := pq. Finally, we define the functions p\^^^ and C 



>{At} 



by 



(4.2) pr^|(t„_,,„) 



1 

Xt 



1 

At 



apf'^)dt, n = l,...,7V. 



Using the above notation, ( 3.26a -c) summed for rt = 1, . . . , can be restated in the form: find 



{pf*\t),u^\t),i;^\t)) e T X V X {X n Z2), t e [0,T], such that 



7^ 

\ dt 



W (fl) 



T 

Jn 



dt- I I p^^'^u^''- ■V^7]dxdt = 



(4.3a) 
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(4.3b) 



(Plul) 



\ / / Pt'' 

Jo Jn 



2 dt 

I At,~ 



; da; - 



V7 \ At,+ 



w 

K 



Jn 

At,- 



At. A 



dx dt 



f I pf'+/^*^+-z.dx-fcX^ / C,(A^C(pf'+)^f'''):V.u;dx 
Jo Jn ~ ~ ~ « ~ ^ ^ 

Vw e Li(0,T;F), 



dt 



lo JnxD 



M |(C(pl) ^l)^* ^ + a E E • V 



K K 



~At,+ At,- *{At} 7At,+ 



JnxD 

T 



i=l 3 = 1 



dq dx dt 



"S/ xf dq dx dt 



Jo JnxD ,_i 



/ At.+ N 

a[u^ )q. 



(4.3c) 



C(pf'+)/3^(V^f '+) • V,^^dgdxdt = 
V(/7e Li(0,T;X); 



with g e (2,oo) when d — 2 and g S [3,6] when = 3, subject to the initi al con dition s /9f^* (0) — 
po e T, *(0) = u° e F and V'l*(0) = G ^ H ^2, where we recaU and ( |3.18[ ). We 

emphasize that (4.3a-c) is an equivalent restatement of problem (P^*), for which existence of a 
solution has been established (cf. Lemma 3.5). 



In conjunction with defined by (1.22), wc consider the following cut-off version of the 
entropy function F : s ^ ]R>o H> F{s) — s(log s — 1) + 1 £ M>o: 



(4.4) 

Note that 

(4.5) 

and 

(4.6) 

Hence, 
(4.7) 



s(logs-l) + l. 



< s < L, 



2L 



+ s(logL-l) + f, L<s. 

log s, < s < L, 

f+logL-1, L<s, 



(F^fis) = 



< s < L, 
L < s. 



P'^is) = min(s,L) = [(^^)"(s)]-\ .s G ]R>o, 



with the convention l/oo := when s = 0, and 



(4.8) (-F^)"(s) > -F"(5) = se]R>o. 
We shall also require the following inequality, relating to T: 

(4.9) J"^(s) > J"(s), s e 



For < s < 1, ( |4.9[ ) trivially holds, with equality. For s > 1, it follows from ( |4.8[ ), with s replaced 
by a dummy variable cr, after integrating twice over a e [1, s], and noting that (J^^)'(l) = J^'(l) 
and (J"^)(l) = J"(l). 
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4.1. L-independent bounds on the spatial derivatives. We are now ready to embark on the 



derivation of the required bounds, uniform in the cut-off parameter L, on norms of G T, 



(t) e V and '4>f^^'^{t) € XnZ2, t G (0,T]. As far as /o^*^'^ is concerned, it follows by tracing 



At,+ 



the constants in the argument leading to inequality (17) in DiPerna 
from (|3.3|) that po G T, that, for each p € [1, oo]. 



Lions [50] and recalling 



(4.10a) 



\\pT\t)\\L.in)<\\po\\LHn), ie(0,r]. 



and therefore, for each p G [1, oo] 
(4.10b) 



pt''^{t)\\LHn) < IbollLp(a), te (0,T]. 



At. 



Concerning u'^'^'^ , we select w = X[o,t] 

n e {1, . . . ,N}. We then deduce using the identity (3.631 with v 



noting f3^, (3.4), (3.401 and (3.17) that, with t = t„. 



as test function in (4.3b I, with t chosen as 
lit,- 



and 2p = tp' 



(4.11) 



At,+ / 



Pmin / II At,+ At,- 1|2 j 



Jfl 



^,ip^'^+)\Diut''^)\^dxds 



< Po |uopdx- 



„2 (-^2 

/^min Co Jo 
K 



2k f f M(g)X:C/;(ik.P)C(p^+)V'^+'7.g^V.^.^+dgdxds, 
Jo JnxD ~ j^]^ ~ ^ ^ ~ ~ ^ ~ 



where || • || denotes the norm over fl. We intentionally did not bound the final term on the 



bounds on V'l below will furnish an identical term with the opposite sign, so then by combining 



right-hand side of (4.11). As we shall see in what follows, this simple trick will prove helpful: our 



the bounds on w^*'^ and iPl'''^ ^^^^ P^i^ otherwise dangerous, terms will be removed. This 
fortuitous cancellation reflects the balance of total energy in the system. 

Having dealt with u^*'^, we now embark on the less straightforward task of deriving bounds 
on norms of V'l*'"'' that are uniform in the cut-off parameter L. The appropriate choice of test 

'''^*'''~) with t = tn, 



TAt. 



function in (4.3c) for this purpose is (p ~ X[Q.t] {-F^)' [i^i^^''^) with t = tn, n G {1, . . . , iV}; this can 
be seen by noting that with such a iy9, at least formally, the final term on the left-hand side of 



(4.3c) can be manipulated to become identical to the final term in (4.11), but with the opposite 
sign. While Lemma 3.5 guarantees that ip^*'^{t) belongs to Z2 for all t G [0,r], and is therefore 



nonnegative a.e. on fl x D x [0,T], there is unfortunately no reason why -0^ ' should be strictly 
positive onflxDx [0, T], and therefore the expression (T^)' {ip'^^'~^) may in general be undefined; 
the same is true of (J"^)"(V'i*'"^), which also appears in the algebraic manipulations. We shall 



circumvent this problem by working with (J^^)' {tp'^*'^ + a) instead of (J^-^)'(-0^*' 



since V-'l*'"'" is known to be nonnegative from Lemma 

are well-defined. After deriving the relevant bounds, which will involve J-^{il!^*''^ 
shall pass to the limit a — )• 0+, noting that, unhke {F^)' and {F^)"{%pf^ 



3.5 



(J-^)'(^f ^++a) and {F^)" {^'^'-^ + a) 



), where a > 0; 

At,-|- 



a) only, we 
), the function 

(J"^)(iAf*'"^) is well-defined for any nonnegative V'l*'^- Thus, the core of the idea is to take 
any a € (0, 1), whereby < a < 1 < L, and choose — X[o.i\ (-^^)'('/'l + «); with t = 
n e {1, . . . , A^}, as test function in (4.3c), and then pass to the limit a — >■ 0+. An equivalent but 
slightly more transparent approach is to start from ( |3.26c ) with the indices n and n — 1 in (3.26c) 
replaced by k and /c — 1, respectively, choose ip — (J^^)'(^/;^ -j-a) as test function, sum the resulting 
expressions through k = 1, . . . ,n, with n S {1, . . . , N}, and then pass to the limit a — > 0+. For 
reasons of clarity, we shall adopt the latter approach. 
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Thus, for k = 1, . . . ,n and n G {1, . . . , N}, we arrive at the following identity 



At 



M 



nxD 



At 



nxD 



K 

E 



-i.(v,[(^^)'(V'i + a)]) ^idg 



dx 



4A 



/Tfe V7 \ I TrL \i I 



V,J(^^)'(^!+«)]dgda; 



+ / eMV,^2-V,[(J-^)'(V£ + a)]dgdx = 0. 



(4.12) 



We shall manipulate each of the terms on the left-hand side of (4.12 1. We begin by considering 

M' 



At 



[(J-^)'(^i + a)] dgdx 



and 



M 



flxD 



y.[(-F^)'(v^2 + ")] rLdqdx 



in tandem. By noting that, thanks to (|4.7[), 



(V'i+a)-a [/3^(^^+a)]-iV,(^i + a) 



f;^(^^+a)-a[(J-^)'(^i + a)] 



where 



s-1, 

^ + ^ - 1 

2L ^ 2 ^' 



s < L, 
L<s, 








"g^(^i + a)-a[(J-^)'(Vi + a)]' 




dj 




J D 









By applying ( [332| with <yJ = M[g^{%j}l + a) ~ a (T^Yi^'l + a)] dq G where g G 

(2, oo) when d = 2 and q G [3, 6] when d = 3, we then have that 



To = - 



M 



flxD 



dpi) ^ apt') 

At 



We note in passing that the statement M[Q^{'il)\ + a)-a{F^)'{il^\ + a)] dq G W^'"-^ (f2) above 
follows, for all a G (0, 1), from the following considerations. Since V'i G -'^^j also € X, 

and hence j M {J^)' {iil + a)Aq G i7i(r2) C for ah q G (2,cx)) when d = 2 and 

q G [3,6] when = 3. Furthermore, since ^/^l G X, also T := JjjMQ^{if'l + a)dq G i^(r2) and, 
since X C H'^iQ; LIj{D)) and, by the Sobolev embedding Q, H\n;Llj{D)) CI {n; LIj (D)) 
for the range of q under consideration here, the definition of G^, a straightforward application 
of the Cauchy-Schwarz inequality to the integral over D involved in the definition of T and the 
application of Holder's inequality to the integral over ft involved in the definition of the (ft) 
norm, imply that Va^F G L'^{U). Finally, by a Gagliardo-Nirenberg inequality applied to the 
function F - /j.jFdx (cf. (2.9) and (2.10) on p.45 in [33]; or inequality (2.19) on p.73 of [31] in 
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conjunction with Poincare's inequality in the L^-i (17) norm; or Theorem 2.2 and Remark 2.1 in 
[52] , where the proof of the GagUardo-Nirenberg inequahty can also be found) : 



rda; 



<c{q,d)\\r\\lti,^^\\Y.T\r^ 



Li-i (fj) 



hence T e L^{i}), which, together with G L^{il), implies that T e W^'^{n). Thus we 
deduce that MlG^i^^ + a) - a (J"^)'(V'i + «)] dq G W^'^{Q), where q € (2, oo) when d = 2 
and q e [3, 6] when c? = 3, as was claimed above. 
By rewriting Ti as 



Ti = 



M aA-') + + a)] d, d. 



M 



At 



and adding this to the expression for T2 yields that 



To = 



Afc(p: 



M 



At 



C{pim+a)-apl-'){A-'+a) 



At 



[{T'^y{^j1+a)]dqdx 
{^l + «)[(-F^)'(^^^ + a)] ~ + a)] dg dx 

[(J-^)'(V^i+a)]d(zdx 



^i^i^^^t;^(V^i + a)dgdx. 



By applying part c) of Lemma |3.3| with F(s) = J"^(s), G(s) = ^^(s), A = C{pl), B = C(Pl~ )' 
a = V'i + a, 6 = -0^"^ + a, noting that s (J"^)'(s) - -F^(s) - Q^{s) = := cq for aU s e (0, cx)), 
and ess.infs>o(J-"'^)"(s) = 1/L := do, it follows that 

1 



T2 > 



At 



MC(pi)^^(^i + a)dgda; 



flxD 



nxD 



Mapt')^\A-' 



a) dq dx 



(4.13) 



1 



2AtL 



MC{pl''){i^l-tP'r')'dqdx. 



nxD 



We now move on to the next term in (4.12): thanks to (2.6), we have that 

K K 



T3:=A / EE^»jMV,^V^i-V,J(J-^)'(^i+a)]dgdx 

- / M [(-F^)"(^i + ")] E E Y^.^i • Y^^^' 



(4.14) 



> 



4A 

Oo 

4A 



SlxD 



1=1 j = l 

Af [(J-^)"(^i+a)]|V,V^ipd(zdx. 



It is tempting to bound [{J^^)" {'tp'l + a)] below further by (i/)| + a) ^ using (4.8 1. We have 
refrained from doing so as the precise form of (4.14) will be required to absorb an extraneous term 
that the process of shifting ip'l by the addition of a > generates in term T5 below (cf. the last 
line in (4.16)). Once the extraneous term has been absorbed into the right-hand side of (4.14), 
we shall apply inequality (4.8) to the resulting expression to bound it below further. 
The next term that has to be dealt with, this time by a direct use of (4.8), is 

(4.15) Ti-.^e f A/y,^^y,[(J■^)'(^^ + a)] >£ / Af dg dx. 



nxD 



QxD -ipf + a 
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It remains to consider the critical final term 



Tr := - 



(4.16) 



K 



J2 M^D <l^] M C(pi) • V,. [(J-^)'(V'! + a)] dq dx 



K 



(IxD 



J2 9.] M dpi) + a) ■ V,, [(-F^)'(V^i + a)] dq dx 

i=l ~ ~ ~ ~ ~ ~ 

K 



nxD ; 



apt: 



dx 



MCipt 



nxD 



K 



^[(V, ui)g,].V,,Vidg da 



i=l 



Thus, by applying the integration-by-parts formula (3.15 1 to the expression in the square brackets 



in the penultimate line of (4.161, we deduce that 



- / MY^U'{\\q^'')aP^)i^t[{q^ql)■■'^.ul]dqdx 



(4.17) 



MCipt 



flxD 



/?^(Vi+a) 



K 



^[(V,ui)g,] • V.^Vidqdx. 



By summing (4.131, (4.141, (4.15) and (4.17) we obtain 



(4.18) 



1 

At 



nxD 
1 



AfC(pi)J-^(V^^+a)dgda;- / Af ((Pl"') -F^(^l-' + «) dg dx 



nxD 



2AtL 



MapT')i^j'l-i^l-')'dqdx 



nxD 



+ 



4A 



M [(J-^)"(7/^i + a)] |V,Vir dqdx + e I M 



k 1 2 



SlxD 
K 



nxD tp'l + a 



■ dq dx 



< I U'{^\q^\^) C(pi) V't [{q^ qj) : V, ul] dq dx 

inxD ^ ~ 



1=1 



Mcipi) 



(IxD 



K 



^[(V^ M|)g,] • V,,^^ dq dx. 



i=l 



As each term in (4.18) can be seen as the value of a piecewise constant function on the interval 



(tfc_i,ife), multiplication of (4.18) by At and summation over the indices k — l,...,n, where 
n e {1, . . . , A^}, yields on noting that ■0f'*(O) =4)'^ = /3^(V'°), for t = tr,, that 
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inxD 



MC(pf '+ W) + a) dq dx 



1 



2AtL Jo 
4Aio 



nxD 



M ■+ + a)] ly.^r'^l' dq dxd. 



At,+ |2 



f2x_D 



< 



M^I^^^dqdxds 
JnxD Vl ' + 

MC(po)^^(/3^(V^°) + a)d(? dx 



ilxD 



Jo JnxD .-_n 

t 



(4.19) 



JnxD 



1 



TAt,- 



1=1 



At.+ N 



y^i^'i dqdxds. 



We refer to [HI (4.14)-(4.18)] for the details of a similar, but somewhat simpler, argument in the 
case of C = 1. The denominator in the prefactor of the second integral motivates us to link /St to 
L so that AtL = o(l) as At— >0+ (or, equivalently, At = o{L~^) as L — oo), in order to drive the 
integral multiplied by the prefactor to in the limit of i — > oo, once the product of the two has 
been bounded above by a constant, independent of L. 

Comparing (4.19 ) with (4.11 ) we see that after multiplying (4.19 ) by 2k and adding the resulting 
inequality to ( 4.11[ ) the final term in (4.11 ) is cancelled by 2k times the second term on the right- 
hand side of (4.19 1. Hence, for any t = tn, with n G {1, . . . , N}, we deduce that 

ft 



pt''\t)\ur^^{t)\-^dx+ 



Pmin 

'Af 



At,+ 



At,- ||2 J 

Ul II ds- 



At, + \ 



\Piu' 



At,+ \|2 



dx ds 



Jn 



-2k 



M(:{pf'+{t))F\i^f'+{t)+a)dqdx 

Cm in k 



AtL 



JnxD 
t 

2k e 



M(V'i*'+-V'L*'~)^dg dxds 



Af ^^^^T-f -dq dxds 



< Po lyoP dx + 



2k 



JnxD ij}']^ 

an k 
H — — 
2X 

rmax ^>ir 
Mmin Cq Jo 

At 



Af,+ 



JnxD 
t 



M (J-^)"(^^*^+ + a) ly^^Z-f '+P d£ dxds 

MC(po)^^(/3^(^°) + a) dq dx 



nxD 



MCipl 



JnxD 



At, + \ 



/3i(V'f'+ + a) 



y^.'0^*'~'' dq dx ds. 



(4.20) 



Let b :— {bi, . . . , bx), recall (3.3), and 6 := |6|i := bi 
of the last term on the right-hand side of (4.20) by 



-bji, then we can bound the magnitude 



ao k 
~4X 



JnxD 



M(J-^)"(V'f'+ + a)|V,Vf'''l' dqdxds 



(4.21) 



4Afc6C, 



2 

max 



ao 



Mf*'+Pdxds 



Jn 
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see [ini (4.20)] for details in the case of C = 1- Noting ( |4.21| ), p.40[ ), and using (|48]) to bound 
+ from below by + = + and ^ to bound 



by F{'4>^*''^ + a) from below yields, for all t = t„ , n <E {!,..., N}, that 



n 



Pmin 







1 At,+ At.-||2 



ds+ f I /.(pf-+)|Z)(yf^+)pdxds 
Jo Jo 



+ 2A I AfC(pf'+(t))-F(^f'+(<)+a)dgdx+^ 



+ 2fce 



/ ;At.+ I 

JnxD Vl ' + 



A \ I I TAt.+ , 

Jo JoxD tpi^ + a 



M - tjj^'-^f dq dxds 

JnxD 

M'-^;i4^ L dgdigds 



Jn A^minCo Jo ~ 



(4.22) 



2A: 



/^min Co JQ 

M C(po) J"^(/3^(V'^) + a)dqdx + a 



ll|?(yr^)lPd.. 



Next we note that an analogous argument to the one that was used to derive W, (4.25)] yields that 



(4.23) 



nxD 



AfC(/9o)-^^(/3^(^") + a)dgda;< ^Cmax|f^|+ / Af C(po) ^(^" + «) dg dx. 



The only restriction we have imposed on a so far is that it belongs to the open interval (0, 1); let 
us now restrict the range of a further by demanding that, in fact, 



(4.24) 



< a < min 1, 



Mmin ao Co 



where cq is the constant appearing in the Korn inequality (3.40). Then, the last term on the 
right-hand side of (4.22 1 can be absorbed into the third term on the left-hand side, giving, on 
noting (|3.40p and (|4.23|, for t ^ <„ and n e {1, . . . , N}, 



(4.25) 



pt'-+{t)\ut'-^{trdx + 
•t 



Pmh: 



JSl 



KPl ) - a 



Ai Jo 



At,+ At,-||2j 

-"l II ds 



2 

max 



Co Co 



|I?(uf ^+)pdxds 



2k 



MC{pt''+{t))Hi^t''^{t)+a)dqdx 



Cmin ^ 

AtL 



M(V^^*'+ - VL*'")^dg dxds 



JnxD 



2k e 



M 



, At , dq dx ds 

JnxD ^pj^ ' + a ~ ~ 



oo k 



M 



4A Jo JnxD V'l*'^ + 



dq dx ds 



< Pq \uofdx + 



„2 ^2 ft 

rinax ^ >t: 



Mmin Co J q 



11/^*^+11^(0) d. 



+ 3afcCmax|f^| + 2fc / MCipo)J^ii;° + a)dqdx. 

JnxD ~ ~ 



The key observation at this point is that the right-hand side of (4.25) is completely independent 
of the cut-off parameter L. 
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On noting |H1 pp. 1243-44], we can pass to the limit a — > 0+ in (4.251 to obtain, for all t = t„, 
n e {1, . . .,N}, that 



At,+ /,x|2 . 



Pmin 



1 At,+ At,-||2 



d5 + r / /^(pf'+)|I)(yf'+)pd3;ds 
Jo Jn 



'2k 



Cm in ^ 



M - V'l *'~)^ dg ds 



JnxD 



+ 8ks 



JnxD 
t 



JnxD 

2 



Af |y,A/tAf'*'+P dgdjds 



(4.26a) < / pol^ol'dx+ / 

Jn Mmin Co Jo 

< / polyoI'dx+^-'i^^^ , 

Jn Mmin Co Jo 

(4.26b) 



/^*'+lli^(o)d.s + 2fc 



MC(po)-^(V'°)dgd5 



SlxD 



l/lli><(n)ds + 2fc 



AfC(po)-^(^o)dg d^ 



nxD 



--: [B(yo,/,V'o)]', 



where, in the last line, we used (3.20a) to bound the third term in (4.26a|, and that t G [0,r] 
together with the definition (3.271 of /'^*'+ to bound the second term. 



4.2. L-independent bound on the time-derivative of p^*'. In this section we shall derive 
two L-independent bounds on dp^^^^ /dt. We begin by observing that 



(4.27) 



[At] 



dt 



dt = - 



W '9-1 (O) 



Jn 



Pl 



T^^^dxdt V77eCo°°(0,T;W^i'i^(f])), 



dt 



which implies that coincides with the derivative of 

^[At] ^ L°°{Q,T]L°^{Q.)) C P'(0,r;VF^'^(17)') 

with respect to t in the sense of M^^'^(f7)'-valued distributions on (0,T). 
We deduce from (|43a]) that, for any 77 e 2.^(0, T; W^'^{^)), 



7^ 

\ dt 



dt 



(4.28) 



< l|p!^^*'l|L-(0,T;L~(f2)) ll^L*' IU^(0,T;Li(n)) llVj;??]] 



L2(o,T;L9-i (n)) 



where g G (2, 00) when d = 2 and g G [3, 6] when d = 3. 



Next we shall show that the expression on the right-hand side of (4.28 1 is bounded, independent 



of L and At. Indeed, we have by (3.47) that p^*' G T and therefore Hp'^"^' ||l~(o,T;L°°(o)) 5; Pn 
Thus, by (|3.40|), (lO and (|4.26bp, we then have that 



JAt] 



II V:, 

(4.29) 



At,+ ||2 



lL2(o,T;L2(o)) 



< 



Co Mmin Jo Jn 



p{pt''^)\D{u'^''+)\'dxds< 



Mniin 



[B(uo,/>o)]' 
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and, by (3.171, (4.26b) and (4.291, we have 



IV7 At.-||2 



L2(o,T;L2(o)) 



0||2 



= At||V.,u"||2 



T 

At 
T-At 



Ai, 



•ds 



A^+||2 



ds 



(4.30) 



Bfeo,/>o)]' 



At.+ ||2 



^0 Mmin 



Thus we deduce from (4.291, (4.30) and Poincare's inequahty that \\u^*''"'^\\L'^(Q^T:H^(n)) ^ C**: 
where C* is a positive constant that depends solely on e, pmim Pmax, Mmin, Cmin, Cmax, 
ao, Co, Cjt, fc, A, K and 6; in particular, C* is independent of L and At. Hence, by Sobolev's 
embedding theorem. 



(4.31) ||uf*^''^^||L2(o,T;L-(n)) < C'*, for ah 

and hence also for all s = q, with q as above, where C* is independent of L and At. Using these 
bounds, we deduce from (4.281 that 



se[l,oo) if<i==2, 
se[l,6] ifd = 3. 



dp 



[At] 



at 



dt 



W '9-1 (O) 



V77 e i2(0,T;M^i'i^(J7)), 



(4.32) 



which then implies that 9p^*'/9t is bounded, independent of L and At in the function space 
L2(0, T; M^^'i^ (rj)'), where q G (2, cx)) when d = 2 and q G [3, 6] when d = 3. 

Note also that, for any 77 e ^^(G, T; iJi(17)) C £^(0, T; H/^'?^ (17)), equation (|4.3a[) yields 



at 



,77) ^ dt 



< l|Pjf*'llL~(0,T;L~(a)) IIMl*' |1l~(0,T;L2(o)) ||y:r?7lUi(0,T:L2(n)), 



(4.33) 



where q G (2, 00) when d — 2 and g G [3, 6] when d = 3. Also, p^*'^ > /Omin(> 0) a.e. on x (0, T], 
which, together with (4.26b), implies that ||w^*'^||L°°(o,T;L2(n)) < C**, where C* denotes a positive 
constant that is independent of L and At. Since, by (3.16), < ^^^||mo|Pj it follows that 

l!i°°(o r;L2(si)) < C*, where C* is a positive constant, independent of L and At. Thus, for 
any 77 G Li(0,T; iJi(f7)) C Li(0, T; M^^'^ (!!)), 



(4.34) 



5t 



dt 



W '9-1 (o) 



where 9 G (2,cx)) when c? = 2 and q G [3,6] when ct = 3; i.e., 9p^*'/9t, when considered as an 
element of [L\0,T; H^{fl))y = L°°(0, T; 7Ji(f])'), is bounded, independent of L and At. 

In summary then, we have shown that 9/9^*'/^^ € V'{0,T;W^^^ (n)') is bounded, indepen- 
dent of L and At: 

(i) when considered as an element of the function space L^{0,T;W^'^{ny), where q G 
(2, 00) when d = 2 and q G [3, 6] when d = 3; and 

(ii) when considered as an element of the function space L°°{0,T;H^{iiy). 



EXISTENCE OF GLOBAL WEAK SOLUTIONS FOR DILUTE POLYMERS 



43 



4.3. L-independent bound on the time-derivative of u^*. We begin by deriving some pre- 
liminary L-independent bounds on tlie time-derivatives of u^* and [plul)^^- On noting from 
and \A.\h\ that 



At 



t - t. 



n-1 At,+ f-n 



At,. 



an elementary calculation yields that 

r 

6 

(4.35) 



', t G (t„_i,t„], n = 1, ... ,7V, 



At,- ||2 



ds < 



< 



V7 At,+ ||2 I IIV7 Af,-||2 
^xU^ II + lly^rtt^ II 



Note that, by (3.471, pmin < Pl*^'^"^ ^ Pmax a.e. on J7 x [0,T] for all At and L. Hence we have 

that ||pf *^''''^||L°°(0,T;L°°(n)) < As 



At --^^ ^ At --^ ^ 



for all t G [tn-i, tn] and n 1, . . . , A'^, which in turn implies that 
\\{pLUL)''\;t)\\LHn) < l|pf'+(-,i)l|L~(n)||yf (•,0IIl=(O) 

+ {\\pt'''i;t)\\L^in) + \\pt''^i-ML^m) Wut'^-i-ML^m 

< Pmax||uL*(-,i)||L-(n) + 2pmax || U^*'" (• , t) || L= (fj) 

for all t G (tn-i, tn] and n = 1, . . . , N. By squaring both sides, using the algebraic-geometric mean 
inequality and integrating the resulting inequality over t G [0, T], we deduce on noting (4.31 1 that 



(4.36) 



ll(PLML)'^*l!L2(n,T;L=(n)) <C*, for 



s G [l,oo) if d = 2, 
sg[1,6] ifd = 3, 



where C* is a positive constant, independent of L and At; (4.31) and (4.361 then imply that 



and are bounded in H-\0,T; L'^ifl)), where 

independent of L and At. 



s G [l,oo) if d = 2, 
sG[1,6] ifd = 3, 



We shall use (4.3b) to improve the bound (4.36). To this end we first note that using (4.3a I in 



(|4^ yields, for ah w G L\0,T; V), that 

fl 



^ ' --At „.. 1 [At] At,- V7 r.At,+ 



dt 

■ T 



21 

„At,+ ^ r^f_ At,+ N 



u. 



dx dt 



+ / / pipt'nD{urn--DMdxdt 

Jo Jn ~ ~ ~ ~ ^ 

i-T , 

{At} 

Pl 



Jn 



I At,- v7 N At,- 



At,+ 



da; dt 



(4.37) = r / pf'+/'^*'+-ii'dx-fc^ / C,(AfC(pf'+)V^f'+):V,u;dx 

Jo Jn ~ ~ ~ „_i in ~ ~ ~ ~ 



dt. 



It follows that, for any t',t" G [0,r] such that < t' < t" < T and choosing w = Xlt'.t"]^' with 
t; G If , where, as in the discussion following (3.44b), X[t'.t"] denotes the characteristic function of 
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the interval [t',t"], we have 



t' Jn 



^' -At „. 1 [At] At- w /_.At, 



dx dt 



t' Jn 



fi{pt''+)D{u^*-+):D{v)dxdt 



It' Jn 

rt" 



{At} 

Pl 



I At,- v7 \ 



At,+ 
L 



At,+ 



dx dt 



- / / pT'^J^'^^-vdx-kY^ f C,(MC(pf'+)^^+):V.z;dx 
Jt' Jn ~ ~ ~ -^^ Jo ~ K ~ ~ 

and hence, equivalently (cf. ( |4.2[ ) for the dcfinitfon of 



dt Vw G V, 



• w dx = — 



£)(i;)da; dt 



(P 



{At} ^[At]^ At - 



At,+ 



t' JO 

f 



rt" 



t' Jn 

V dx dt 

At.H 



Jf Jn ~ ~ ~ ^^-^ Jt' Jn ~ 



(4.38) 



Ui + U2 + U3 + U4 + Us 



We shaU suppose that t" = t' + (S, where 5 G (0,T — t'], and bound each of the terms Ui, 
I = 1, . . . , 5, in turn. We note in particular that, by the definition (4.2 ) of p^^*'^, the term U2 = 
when t',t" e {{) = tn,ti,---,tN^i,tN =T). 

For Ui, by the Cauchy-Schwarz inequality, we have that 



|Ui| < 



t'+& 



\D{u't''^{t))\^dx\ dt 



n 

At, + \ 



\\D{v)\\ 



< PmnxS^ ^ )l|L2(t',t'+5;L2(0)) IIP(^)II G F. 

Further, for any q G (2, 00) when d — 2 and any q G [3, 6] when c? = 3, we have that 



IU2I <Pn 



t'+S 



\\ut*'-m^^,^jY.nT-^mdt 



Wnhiin) G V. 



The Gagliardo-Nirenberg inequality (3.2) and Korn's inequality (3.40) imply that 

L 9 ^ (S2) 



for all g € (2, oo) when d = 2 and g € [3,6] when d = 3. Therefore, by Korn's inequality (3.40) 
again, and by Sobolev embedding and Holder's inequality, we have that 



\V2\<C{d,q,n) 

pmax 



t'+S 



\\D{ut'-m\U\D{ut'^^m\dt 



\\D{v) 



< C{d,q,n) Pniax\\U]^ ' llLoo'(o,T;L2(n)) 



X \\p{ut''-)\\l.^^,^,,+S;LHn)) \\p{uT^^)\\L^t'.t'+S;mn)) \\p{v)\\ G V. 



At,+ 
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Hence, 

PK - )ll 

X s"^ \\piuf''+)\\L2^,,^t'+5;mn)) Wpivjw yv e V, 

for all q e (2, oo) when d = 2 and all q e [3, 6] when d — 3. An identical argument yields that 



IU3I < C{d,q,n)p. 



1 

maxlltti IIl°^(0,T: 



X ^^^i^ ||p(%*'+)||L^(t',*'+i;L2(f,)) ||Z?(^;)|| Vv e V, 

for all q G (2, 00) when d = 2 and all g S [3, 6] when d = 3. 

For U4, on noting (3.3), (3.4), Korn's inequality (3.40) and Holder's inequality (with respect to 
the variable t') yield that 

IU4I < C{co,>c,n)p^,^6-^ ll/^*^+llL^(r,tW(o)) wpivjw Vv e V, 
with >f > 1 if d = 2 and >f = | if d = 3. 



Finally, for term U5, from (|2.4a|), (|3.15|), the Cauchy-Schwarz inequality, we have that 
IU5I 



i=l 
rt'+S 



MC(pf W) W C/;(||g.|') qj -.V.v dq dxdt 



TAt. 



flxD 



-2k 



nxD 



K 



dq dx dt 



t'+5 



< 2fcC„ 

< 2fcCn 



t'+S 



t'+6 



|V,w| 



|V:,w| 



D 



K 



A/|q|2^f-+(0dg 



1 

2 




v,\/^f'+(t) 


2 




JD 







dg da; dt 



dx dt 



M 



dg 



dx dt 



X ess. 



suP(K,t)eax(o,T) 



Migpvf'+dg 



< 2fcCmax52 ||Vq^?AL*''^llL2(t',t'+5;L|,(axD)) IIYxW|| 

X eSS.SUP(2,_f)gj2x(0,T) 



D 



M|gpV.f'+dg 



Since |g|2 = |qi| 



Igxl < fci + • • • + =: ^, the inequality (3.77) immediately implies that 



the final factor is bounded by VaTfo. Korn's inequality (3.401 then implies that 



IU5I < C(fc,^,6,C,„ax,Co)54 ||y,Y'V^L*'^llL2(t^t-+5;L|^(OxD)) WP{v)\\ G V. 

By collecting the upper bounds on the terms Ui, i = 1, . . . , 5, and noting the upper bounds on 
the first and the third term on the left-hand side of (4.26b) in terms of [B(uo, /, "^o)]^: together 



with the uniform lower bounds p^*"^ {x,t) > p^^^^ and p{p'^'^''^ {x,t)) > Prain for {x,t) £ fix [0,T] 
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we thus have from (4.38) that 



<C\\D{v)\\ {S^ +S^)\\D{i 



V dx 



\L-^it',t'+S:L^{n)) + (5^ |l/^*'^||L2(t',t'+5;L''(n)) 



(4.39) 



where C is a positive constant, independent of At, L and S; and 6 = t" — t' E (0, T — t']. 

In what follows, we shall suppose that t',t" e {0 = ^Oj^i, • ■ • ,^Af-i,ijv = T} with t' < t" , so 
that 6 := t" — t' is an integer multiple of At; then 

(4.40) = p^tHn [unn - unt')] + [p^tHn ^ p^tHf)] ut\t'). 



By selecting -q = X[t' ,t"] {t) (yL*(i') 'il) (4.3a), with t" = t' + S, we have using Korn's inequality 
(|3.40|) that, for any veV, 



< 



PmaxtJ^ ||L2(t',t'+d-;L<!(n)) II ("i * (*') •^)ll^^ 



<Pmax'5'|kf^ \\LHt',t'+S;L-.m) f 1 1 V^wf (i') 1 1 1 1 « 1 1 ^ + 1 1 1 1 1 1 (t') 1 1 ^ ^^^^ 



< C5^ 



(4.41) 



X ||p(uf (t'))|| ||^;||^ ||^;||2;(^^) + ||I?(z;)|| (i')ll^ (i')ll2;(n) 



where q G [4, cxi) when d — 2 and g G [4,6] when d — 3, and C C(co,/9max), where cq is the 
constant in Korn's inequality (3.40). 

By dotting ( 4.40[ ) with v £ V integrating over Q and substituting (4.39) and (4.41) into the 
resulting identity, we deduce that 



P^t'\t' + S) [ut'{t' + S)~ut'{t')]-vdx 



<C(52||u^*' \\L2{t',t'+5;Li{n)} 



X ^\\D{ut'im\ WvJ—^ Ml^n) + \\D{v)\\ ll^f (t')ll^ lkL*(t')lll;(n) 
+ C \\D{v)\\ (^(<55 + J^) \\D{ut'^+)U.(,,^,,+s-,LHn)) + \\f^'^^\\mt',t'+5;L'<m 



(4.42) 



+ ||Vg^l/'L*'^llL2(t',t'+5;L|,(f2xD)) 



for all V € V, where q G [4,oo) when d — 2 and q e [4,6] when d — 3, and C is a positive 
constant, independent of At, L and S. Here 6 — iAt, where i = 1, . . . , N — m, and t' — mAt for 
TO = 0, . . . , — 1. The symbol S will be understood to have the same meaning throughout the 
rest of this section, unless otherwise stated. 

We now select v = Uj {f + 5) - uf\t') in ( |442| ), sum the resulting collection of inequalities 
over t' G {0,^1, ... ,r — (5} (denoting the sum over all t' contained in this set by X]t'=o)' ^'^'^ "li-Ote 
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the following obvious inequalities: 



At. 



L2(t',t'+(5:L9(0)) 



<c\ 



At.— II q-2 II At,- 

fL \\L^{t',t'+S:Li{^)y\'^L 

At 2 

fi ' lli,2(0,T;L9(n))ll;P(^L 



<;-4 

I 9-2 

L2(t',t'+5:L'!(0)) 

5-4 
I ?-2 



\\ut\t' + 5)- (t')ll^ < (2 llyf IU^(o.T;L^(o))) , 



|yf (i')ll^ < llyf ll2-4o,T;L^(a)) and (t' + <5)||f^ < \\um^^,.T.,L^^n)V 



the first of which follows by Sobolev embedding and Korn's inequality (3.40), to deduce from 



(4.42) that 



< C 5^ \\u^ ' ll22(^o,T;L9(n)) ll3fi*llL~(o,T;L2(n)) 



At|| q-2 



T-S 



At E II^C 



g-4 

.At, -Ml 



,At /.A 



lL2 



t'=0 



T-5 



T-(5 



HQ)) 



+ <5HI/llL2(t',t'+5;L«(a)) + l|Vg^VL*^^llL=(t\t'+5;L|j(axn)) 



(4.43) 



T-5 



=: vi (V2 + vg) + a At E V4(t') Mt') + Mt') + Mt') , 



where q € [4, oo) when d = 2 and q € [4, 6] when d = 3, and C is a positive constant, independent 
of At, L and S; Vi denotes the expression in the first line on the right-hand side of ( |4.43 ); V2 and 
V3 denote the two terms in the bracketed expression multiplied by Vi; V4(/:') is the factor in front of 
the bracket in the fourth line on the right-hand side of (4.43); and Vc^{t'), ^eit') and ^rit') are the 



three terms in the bracketed expression multiplied by V4(i'). We shall consider each of the terms 



Vi, V2, V3, V4(t'), . . . T^rit') separately. We begin by noting that by (4.31) and by the definition of 



M^* in conjunction with the uniform bounds \\u'^*'^\\]^^(^Q X;L^{n)) < C* (cf. the discussion in the 



paragraph between (4.33) and (4.34|), 



(4.44) 



2 g-4 

Af,-||5-2 II At||g-2 <r r'As 



where q is as above, and C is a positive constant, independent of L, At and S. In the rest of this 
section we shall assume that q € (4, 00) when d = 2 and q G (4, 6] when d = 3. 
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Next, for the term V2, Holder's inequality with respective exponents 2{q — 2)/{q — 4), 2 and 
- 2 for the three factors under the summation sign in this term yields that 



T-<5 



2(9-2) 



T-S 



L2(o)) 



At J2 WDiut'im 



T-5 



t'=0 

1 

q-2 



AtY,\\piut\t' + S))-Diut'im\ 



t'=0 

t'=0"^*' ~ ~ / \ ~ ~ t'=At 



9-4 
2(9-2) 



T-S 



X 



AtJ2\\piut\t' + 5))-piut\m\ 



\ t'=0 / 

=■ V21 V22 V23. 

We shall consider the three factors on the right-hand side of this inequality separately. For the 
first factor, we shall use the following elementary result, whose proof is omitted. 

Lemma 4.1. Suppose that g G L^{0,T), g > 0, is a piecewise constant, left- continuous function 
on the partition {0 — io? ^ij • ■ • j ^iv-ij = T} of the interval [0,T] with step size At = T/N, 
where N e N>i; and let S = £At, where I ^ {1, ... ,7V} (i.e., g{tk) = g{tk-) for k = 1, . . . ,N). 
Then, 



AtY, g{s)As = {Atf J2 J2 9{tk)<i{AtrJ29itk)^S g{s)ds 

t'=0 *' k = l s=k k=l •'^ 

where the < sign can be replaced by an equality sign when £ — I and £ = N . 



N-e+i k+e-i 



On applying Lemma 4.1 with 5 : s G (0, T] i-> (s))!!^ in conjunction with (3.171 and 



(4.26b) we deduce that 



V2i<l^/ \\Diut'''is))rds 



< s^t^ I At\\D{u'>)f + I \\D{u^''+{s))W' ds 



where C is a positive constant, independent of L, At and 5. Analogously, since m^* and u^*'^ 
coincide at all points t' = £At, £ = 1, . . . , N, we have, again by (3.171 and (4.26b), that 

v,, = (At\\Diu^)r+At j2 wpiut'^^mA 



= (^At\\D{u°)f + At 



t'=At 

T-S 



\\p{ut''+{s))rdsj <c, 

where C is a positive constant, independent of L, At and d. For the term V23, we have by the 
triangle inequality, shifting indices in the summation, and noting, once again, (3.171 and (4.26b), 
that 

V23<2i^ (^Atj2j\Diut\t'))ry 

= 2^ [At Wpiu'W + At \\p{ut''+it'mA < C, 

\ ~ ~ t'=At ~ ~ / 
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where C is a positive constant, independent of L, At and S. Thus we deduce that 

V2 = V21V22V23 < CS^^, 
where C is a positive constant, independent of L, At and S. An identical argument yields that 

V3 < C^^iMi), 

and therefore 
(4.45) 



Vi (V2 + V3) <CS^ S-^ii-^-) , 



where C is a positive constant, independent of L, At and S. 

Finally, by the Cauchy-Schwarz inequality applied to the sum starting in the fifth line of (4.43 ), 
the bound on the term V23 above, using Lemma (4.1) with 



and the bound (4.26b), we deduce that 

(4.46) C At^Viit') {V^it') +Ve{t') +Vj{t')) < C (s'^ + s"^) + CS^ + ^ + 06^ + ^ , 



where C is a positive constant, independent of L, At and S, with S ^ £ At, 1= 1, . . . ,N. 
On substituting ( |4.45[ ) and ( |4.46[ ) into ( |4.43[ ), we thus have that 

T-S 



E / pfHt' + S)\ut\t' + 6) 



ut'{t')\'dx 



< C(55+3(^ +C((55 +(5^)(55 < C(5i"i^ with 



qe(4,oo) ifd = 2 
q e (4, 6] if d = 3, 



where C is a positive constant, independent of L, At and 5, with 5 — £ At, I — 1, . . . ,N. 
On recalling that p^^*\t' + 6) > pmin for all t' G [0,T — S], we finally have that 

T-S 

E ||uL*(i' + '5)-UL*(i')ll' <CS^-^ with 



t'=o 



(7e(4,oo) ifd = 2 
q e (4, 6] if d = 3, 



where C is a positive constant, independent of At, L and 6; with 5 = lAt, £ = 1, . . . , N. As 
uf\£At) = u^^-+{£At) and u^\{£-l) At) = uf^'^{£At), £ ^ 1, . . . ,N, and uf*'^ are piecewise 
constant functions on the partition {0 = tQ,ti, . . . , tjv-ii tN — T}, we thus deduce that 



T-S 



(4.47) Ai E ||M^*'=^(t' + <5) -u^*'=^(i')|P < C<5^-^ with 



t'=At 



g e (4, 00) if d = 2 
q e (4, 6] if d = 3, 



where C is a positive constant, independent of At, L and 6, with S — £ At, £ — 1, . . . , N — 1. By 
selecting q as large as possible and taking the square root of the previous inequality, we have that 



(4.48) 



|u^*'=^(- + 5) - -"l*'"^!- )IIl2(0,T-5;L2(j2)) <C5'^ 



for all 6 — £ At, £ = 1, . . . , N — 1, where C is a positive constant independent of step size At, L 
and (5; < 7 < 1/2 when d = 2 and < 7 < 3/8 when d = 3. 

We shall now extend the validity of ( 4.48| ) to values of ^ G (0, T] that are not necessarily integer 
multiples of At. We shall therefore at this point alter our original notational convention for S, and 
will consider 6 = v At, with v E (0, N]. Let us define to this end £ — [v] max{fc E N : k < v}, 

:= v-[v] e [0, 1), and for i e (0,r] let m e {0, . . .,N~£~2} be such that t e (m At, (rn+l) At]. 
Hence, 



(4.49) 



At,± 



{t + v At) = 



,At,± 



{t + £At) 



if t G (to At, TO At + (1 - ?9) At] 



(t+(£+l)At) if t e (TOAt + (l-?9)At, (TO+l)At], 
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which then impHes on noting that 



s e M>o ^ e 



is a concave function, that 



At,±, 



> ''{■ + 6)~u 
< (l-i?) At 



At.±f 
L \ 



\L^(0.T-S-L^(Q.)) 



T~iAt 

E I 

t'=At 



At,± 



T-vAt 



At,±, 



vAt) 



At,±, 



■dt 



(t' +£At)-u2-{t') 



^At 



T-{l+l)At 

E 

t'=At 



\\ut''^it' + {i + l)At)-ut*^^{t')\\ 



< (1 - C {£ Aty + d C {{£ + 1) Aty 
<C[{l-'d)eAt + {>{£ + l) At]'' 

= c[{e + ^) At] = c{vAty' ^ cs'', 

where S = vAt, v e (0,A^]; < 7 < 1/2 when d ^ 2 and < 7 < 3/8 when d = 3; and C is a 
positive constant, independent of At, L and S. The second inequahty in the chain of inequahties 
above follows by applying (4.471 first with S — £ At and then with 5 = {£ + 1) At. 
Consequently, 



At,± 



At,± 



( ■ )IU2(0,T-(5;L2(n)) < C S"' 



for all S e (0,r], where C is a positive constant independent of At and L; < 7 < 1/2 when 
d 2 and < 7 < 3/8 when d = 3. 

Thus we have established the following Nikol'skii norm estimate: 



(4.50) 



|ArT-2(0,T;L2(o)) 



sup S 

0<5<T 



At,± , 



(• + (5) - • )\\mo,T-5;L^n)) < C, 



where C is a positive constant, independent of L and At; < 7 < 1/2 when d — 2 and < 7 < 3/8 
when d = 3. 

Remark 4.1. We note in passing that in the special case of an incompressible Newtonian fluid 
with variable density, when the extra stress tensor appearing on the right-hand side of (1.1c I is 
identically zero, ( |4.50[ ) continues to hold and improves the Nikol'skii index 7 = 1/4 obtained in 
the work of Simon [52, p. 1100, Proposition 8 (ii)] and [521 p. 1103, Theorem 9 (ii)] (under the 
hypothesis f € L^(0, T; L^(f2)) compared with f S L^(0, T; L*'(J7)) assumed here with k> 1 when 

d = 2 and >f > § when d = 3, and under the same assumptions on the initial data uq and po as 
in (3.3 I here). 

4.4. Strong convergence of the sequences {Pi*^'*^}L>i) {Hl*^^ ^''}l>1i and weak conver- 
gence of {iPl'^^'^^}l>i- We begin by collecting a number of relevant bounds on the sequences 

r At(,±)i r At(.±)i , r ,At(.±)i 

{Pl }l>i, {Hl }l>i, and {V'i '}l>i- 



First, we recaU that pf\t) e T for all te [0, T] and from (4.10a) that 



(4.51a) 



llP^'WIUnf^) < IIPollLP(n), te{0,T], 



and therefore, for each p e [1, 00], 
(4.51b) 



t& (0,T]. 



Next, noting (4.1a b), a simple calculation yields that [see (6.32)-(6.34) in |TT] for details]: 
MlV^V^Pdgdxdi < 2 / / M \\V^\fj^\^ + \S/^\f^^\^] dqdxdt, 



(4.52) 



InxD 
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and an analogous result with Vx replaced by Vg. Then the bound (4.26b), on noting (3.3), (4.50), 
(3.40), (4.1a b), (3.17), (3.20a), (4.52) and the convexity of J-, imply the existence of a constant 

Cmax, T, \A\, flo, Co, C^, 



Ci, > 0, depending only on B(uo, /, V'o) and on e, p 
k, X, K and b, but not on L or At, such that: 



mini Pmaxi piniii; Smini 



Mil 



eSS.SUPtg[o,T]llwL 



1 

At 



At,+ 



At,±||2 



At.- 







ess.suptgfo-T] 

T 



\\Ur 
\^dt 

rAt(,±) 



AfT-2(0,T;L2(n)) 



nxD 



1 

AtL 



lo JnxD 



At, 



' dq dx dt 



(4.53) 



JflxD 



At(,±)|2 



dg da; dt 



JnxD 



rAf(,±)|2 



dg dx dt < C* 



where < 7 < 1/2 when d = 2 and < 7 < 3/8 when d = 3. 
Henceforth, we shall assume that 



(4.54) 



At = o(L ^) as L -> 



Requiring, for example, that < At < Co/{L \ogL), L > 1, with an arbitrary (but fixed) constant 
Co will suffice to ensure that (4.54| holds. The sequences {/3^*'}l>i, {Pl*''^''}l>1j {pi^*^}L>i 
{?iL*^'^''}i>i {V'l**'''''^}l>i as well as all sequences of spatial and temporal derivatives of the 
entries of these sequences will thus be, indirectly, indexed by L alone, although for reasons of 
consistency with our previous notation we shall not introduce new, compressed, notation with At 
omitted from the superscripts. Instead, whenever L — !■ 00 in the rest of this section, it will be 
understood that At tends to according to (4.54). We are now ready to embark on the passage 
to limit with L 00. 



Theorem 4.1. Suppose that the assumptions (3.3) and the condition (4.54), relating At to L, 
hold. Then, there exists a subsequence 0/ {(p^*, u^*, '!/'^*)}l>i (not indicated) with At = o{L~^), 
and functions {p,u,ili), with ■0^0 a.e. on D, x D x [0,T], such that 

p e L°°(0, T; T) n C([0, T];LP{nj), u e L°°(0, T; n L^{Q, T; V), 

where p Cz [l,oo), and 

^fieL\0,T;Hlj{nxD)), 
such that, as L — > CXI ( and thereby At — > 0+ ), 

(4.55a) 1 ^ p 

[At] 

Pl -^P 



(4.55b) 
(4.55c) 
(4.55d) 
(4.55e) 

where p G [1, 00); 
(4.56a) 
(4.56b) 
(4.56c) 



At(,±) {At} ^ 

Pl ^ Pl P 
p{pt'^-^^)^p{p) 

apT\ apT'-^\ ci'^'' ^ ap) 



weak* in L°" {0,T; L°° {fl)), 
strongly m L°° {0,T; LP{fl)), 
strongly m L°° {0,T; LP{n)), 
strongly m L°° {0,T; LP{fl)), 
strongly m L°° {0,T; LP{n)), 



At(,±) 



Af(,±) , 

— ?> M 



At(,±) 



weak* in L°°(0, T; ^^(17)), 
weakly in L^{0, T; V), 
strongly in L^{Q,T;U {Vl)), 
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where r G [1, oo) if d : 
(4.57a) 



(4.57b) 
(4.57c) 



V, 
M5 V„ 



2 andr e [1,6) if d 
M5 V 



TAt(,±) 



3; and 

weakly m L^{0,T; Lli{n x D)), 
weakly m L^{Q,T;L^{Vl x D)), 
weakly m L^{Q,T;L^{n x D)). 



Proof. The weak convergence results (4.56a b) follow directly from the first and fourth bounds 



in ( |4.53| ). We deduce the strong convergence result ( |4.56c| ) in the case of u'^*'~^ on noting the 
second and fourth bounds in (4.53 1, (3.11 ), and the compact embedding of V into L'^{^)nH, with 
the values of r as in the statement of the theorem. In particular, with r = 2, a subsequence of 
{u'^*'~^}l>i converges to u, strongly in L^{0, T; L^{il,)) as L oo, with At = o(i~^). Then, by the 
third bound in (4.53 ), we deduce that the three corresponding subsequences { 



}l>i converge 

to u, strongly in ^"^(0, T; L'^{Q)) as L — > oo (and thereby At — ?> 0+). Since these subsequences are 
bounded in L'^{0,T; H^{fl)) (cf. the bound on the fourth term in ( |4.53 1) and strongly convergent 
in L^(0, T; I/^(f2)), it follows from (3.2) that (4.56c I holds, with the values of r as in the statement 
of the theorem. Thus we have proved ( |4.56a -c). 

The convergence result (4.55a) and the fact that p e L°°(0,T;T) follow immediately from 
(|4.51a[) and as *(t) G T for all t e [0, T]. Further p.26a[) imphes that 



Pl 



[At] dr] 



dt 



dx dt — 



t„_i Jn 



Pl 



[At] , At, 



V xTj dx di = — 



with q S (2, oo) when d = 2 and q G [3, 6] when d = 3. 
and noting that p\ = po, we then deduce that 



Pl 



[At] At,- 



V7/e Ci([t„_i,t„];W^i'^(r!)), 
Upon summation through n = I, . . . , N, 

■VjjTydxdt^ / porjdx 



yj] £ C\[0,T];W^-^ (n)) s.t. ry(-,r) = 0, 



with q G (2,oo) when d = 2 and q G [3,6] when d — i. Hence, on letting L — ^ oo, with 
At = o(L^^), and noting (4.55a) and (4.56c), we deduce that 



p -TT- dx dt 
^ dt ~ 



pu ■ V xV dx dt 



Poridx 



(4.58) 

with q G (2, oo) when d 



V77G Ci([0,T];M/i^^(r!)) s.t. 7j{-,T) = 0, 



2 and g G [3, 6] when d = 3. Thus we have shown that p is a weak 



solution to (l.lab). One can now apply the theory of DiPerna & Lions [20] to (l.lab). As 
u G L'^{0,T;V), it follows from Corollaries II. 1 and II. 2, and p. 546, in [5D], that there exists a 
unique solution to (l.lab) for this given u, which must therefore coincide with p. In addition, 
p G C{[0,T],L'P{i})) for p G [1, 00), and the following equality holds: 

(4.59) llp(i)IUp(o) = IIpoIUp(o), te(o,r]. 



Thanks to (4.51a) and (4.55a), by the weak* lower semicontinuity of the norm function, and (4.59) 



with p = 2, we have that, for a.e. t G [0, T], 

||p(t)|p<liminf||pfl(t)|p<||po|p 

i^— >-oo 

This then implies for a.e. t G [0,r], 



(4.60) 



Mm 



lim \\pl 

1/— )-oo 



[At] 



\Po\ 
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Thus we have proved (4.55b I in the case of p = 2, which, on extracting a further subsequence, 
impUes that 



(4.61) 



r [At] 
hm p'l ' = p 



a.e. on r2 X (0,r). 



By recalhng (3.471, we then deduce (4.55b| for slip e [1, oo) by Lebesgue's dominated convergence 
theorem. 

It follows from ( |3.48a[ ), with ij = X[i„_i,t) V', t e (t„_i,i„], and r] = X(t,t„] t & [i„-i,tn), and 
ip g W^''^ (fl), where q € (2, oo) when d = 2 and q e [3, 6] when d = 3, that 



(4.62) 



I [At] 

\Pl 



At.±| 

Pl I 



L^{0,T;W 



< C max 

(f2)') n=l,...,Af 



At. 



dt < C{At)^, 



where we have noted (4.51a), (3.2 1 and (4.53). Similarly, on recalling (4.2), we obtain 
(4.63) 



[At] 

\Pl 



oAf II 

Pr 1 g 

'^^ "L°=(0,T;W ' 9-1 (fj)') 



and on noting (3.31 1 
(4.64) 



iia^D-cr^ii 



1^ <c(Aty- 

JAf] 



Applying (4.55b) with p ~ q, we have that converges to p strongly in L°° {0,T; L'i{Q)) — 

L°°{0,T;L^{ny) C L°°(0,r;VF^'^(f7)'), where 9 e (2, 00) when d = 2 and g G [3,6] when 
d = 3; thus, p^*' converges to p strongly in L°°{0,T;W^'^ (fl)'), where q G (2, 00) when d = 2 
and q G [3, 6] when d = 3. By means of a triangle inequality in the norm of L°°(0, T; W^'^ i^Y) 
and noting (4.62), we thus deduce that p^*'^ converges to p strongly in L°° {0,T;W'^''^ (il)'), 
where q G (2, cx)) when d — 2 and q G [3,6] when d = 3; analogously, using (4.63) this time, 
p^* and p\f^^\ converge to p strongly in L°°{0,T;W^'^ (fl)'), where q G (2, 00) when d = 2 
and q G [3,6] when d = 3. Since, thanks to ( |4.51b| , {p^*''''}l>i, {Pl*}l>i and {p|^'^*^}l>i are 
weak* compact in i°°(0, T; L°°(f7)) C L°° {0,T;W^'^ (n)'), where q G (2, 00) when d = 2 and 
q G [3, 6] when d = 3, it follows that the weak* limits of the weak* convergent subsequences 

extracted from {p^*'^}l>i, {Pl*}l>i and {Pi'^*^}L>i have the same limit as p^*': the element 
pG L°° {0,T; L°° (il)). The weak* lower semicontinuity of the norm function and inequality (4.51bl 
together imply that 

|lp(i)|l2<liminf||pf('±)(t)|l2<|lpo|l2 = |lp(t)|l2. 
Hence, proceeding as above in the case of p'^*', we deduce (4.55c). 



Concerning (4.55d) and (4.55e), these follow from (4.55c) and (4.64) via Lebesgue's dominated 



convergence theorem by possibly extracting a further subsequence, thanks to our assumptions in 
(3.3) on p and C. 

We complete the proof by establishing ( 4.57a| -c). According to (4.53) and (3.3), 



2fcCn 



nxD 



M J^{i,^*^'^\t)) dq dx < [B(uo, /, i'o)? 



for all t G [0, T]; hence, on noting that s log(s -j- 1) < 2 [J-{s) + 1] for all s G M>o, we have that 

[B(uo,/>o)]'+2|r!|. 



(4.65) max 



M (i) logi^t""'^^ it) + l)dqdx< ^ 



As s G M>o !—>■ slog(s + 1) G M>o is nonnegative, strictly monotonic increasing and con- 
vex, it follows from de la Vallee-Poussin's theorem that the sequence of nonnegative functions 

rAt(,±) 



{'^l^'''^'}l>i, with At = o{L ^) is uniformly integrable on 17 x D x (0,r) with respect to 
the measure dv := M{q) dq dxdt. Hence, by the Dunford-Pettis theorem, {■(/'l **''*^}i>ii with 
At = o(L-i), is weakly relatively compact in L^{n x D x {0,T);v) = L^{G,T]L\j{n x D)); i.e.. 
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there exists a nonnegative function V G L^i^jT; L]^{n x D)) and a subsequence (not indicated) 
such that 



(4.66) 



TAt(,±) 



tj) weakly in Li(0,T;L}^(r2 X £>)), 



as L — T' CXI, where At = o{L The fact that the hmits of the subsequences of {V'l *^'*''}i>i 



the same follows from the sixth bound in (4.53). Thus we have shown that (4.57a) holds, and that 
the limiting function "0 is nonnegative. 



The extraction of the convergence results (4.57b c) from (4.53) can be found in Step 2 in the 
proof of Theorem 6.1 in |^. □ 



In the next section we shall strengthen (4.57a) by showing that (a subsequence of) the sequence 
{V'l *^'^''}l>i is strongly convergent to ■0 in -^^(0, T; L\.i{Q. x D)) as L — oo, with At = o(L~^). 



4.5. 



Strong convergence of ^A^*'''*' in L^[Q^T;L\j{p, x D)). In Section 



4.3 



we derived an L- 



independent bound on the Nikol'skii norm, based on time-shifts, of the sequence {; 



At(,±) 



L>1 



of approximate velocities. The Nikol'skii norm bound ( |4.50 ) was used in the previous section in 
conjunction with the bounds on spatial derivatives of {uif'^'^^} l>i established in Section U.ll to 
deduce, via Simon's extension of the Aubin-Lions theorem |51j . strong convergence of {y^*''^ /i>i 
in L?'{Q,T;U{9?)) a.s. L ^ oo, with Ai = o(L"^), for 1 < r < oo when d = 2 and 1 < r < 6 when 



2, which we shall then use to pass to the limit in nonlinear terms in (4.3b) in conjunction with 



weak convergence results for the sequence, which suffice for passage to the limit in those terms in 
(4.3b) that depend hnearly on {u^*''^'^^}i 



In [9] we used a similar argument for the sequence of approximations to the solution of the 
Fokker-Planck equation, except that due to the form of the KuUback-Leibler relative entropy and 
the associated Fisher information in the bounds on spatial norms of the sequence resulting from 
our entropy-based testing, (which, in turn, was motivated by the natural energy balance between 
the Navier-Stokes and Fokker-Planck equations in the coupled system, that manifests itself in 
a fortuitous cancellation of the extra-stress tensor in the Navier-Stokes equation with the drag 
term in the Fokker-Planck equation in the course of the entropy-testing), we had to appeal to 
Dubinskii's extension to seminormed cones in Banach spaces of the original Aubin-Lions theorem 
to deduce strong convergence of the approximating sequence of probability density functions. 

Unfortunately, in the present setting, the appearance of the nonlinear drag C,{p) in the Fokker- 
Planck equation obstructs the application of Dubinskii's compactness theorem, and the approach 
based on Nikol'skii norm estimates, that was used in Section [473] in the density-dependent Navier- 
Stokes equation, also fails, because — in order to compensate for the rather weak spatial control 
in (4.26b) of the KuUback-Leibler relative entropy and the Fisher information — its application 
ultimately requires a uniform L°°(0, T; L°°(r2 x D)) bound on the sequence of approximations 
to the probability density function, which is not available. We shall therefore adopt a different 
approach here. Since the argument below that finally delivers the desired compactness of the 
sequence {V'l*''^''}l>i (with At = o{L~^) as L — >■ oo), in L^{Q,T]L\j{Vt x D)) is long and rather 
technical, we begin with a brief overview of the key steps. 

First, by ( 4.56a[ ) (a subsequence of) the sequence {'0i**''^''}-L>i is weakly convergent in the 
space Ll(0,^;L]^^(O x D)) to G L^{Q,T]L\,;{VL x D)) as L oo, with At = o{L-^). We 
shall then make use of the property that if $ is a strictly convex weakly lower-semicontinuous 
function defined on a convex open set U of M, and the weak limit of $(^^*'-'^'') is equal to ^{4>)^ 
then the sequence {0'i*^'*^}L>i converges almost everywhere on (0, T) x fi x L) as L oo, with 
At = o(L-i) (cf. Theorem 10.20 on p.339 of flM). According to ( |4.53| ), 



max 

*G[0,T] 



At(,±) 



dq dx 



is bounded, uniformly in L and At; in addition is strictly convex. Thus T may appear as a 
logical first candidate for the choice of the function $. Unfortunately, we do not know at this 
point if the weak limit of the sequence {-7^(V'l **^''''^)}l>i in -^^(0, T; L\j{fl x D)) is equal to T{ip), 
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and therefore the argument outlined in the previous paragraph is not directly applicable with 
the choice $(s) = -7^(s). We shall therefore make a different choice: we select the strictly convex 
function <i>(s) = (1 + s)^^", s > 0, where a € (0, 1) is a suitable (small) positive real number. We 
note in passing, as this will be important in the argument that will follow, that s (1 + s)", s > 0, 
is a strictly concave function on M>o for a € (0, 1). Although we do not know at this point if, with 
the latter choice of <&, the weak limit of the sequence {^iipL*^'^^)}L>i in L^iO,T; L\,j{n x D)) is 
equal to tf>('0), and therefore this particular $ may seem no better than the original suggestion of 
^(s) = J'{s), we note that by using estimates interior to $7 x 13 on subdomains fig x Dq fl x D 
on which the Maxwellian weight is bounded above and below by positive constants, and therefore 
the uniform bounds in Maxwellian- weighted norms that result from ( 4.26b[ ) become bounds in 
standard, unweighted, Lebesgue and Sobolev norms, one can use function space interpolation 
between these unweighted norms to deduce a uniform bound on the L^^^ {0,T; L^~^^ {flo x Dq)) 
norm of fo'" ^ suitable (small) value of 6, which, with a G (0,5) and an application of 

the Div-Curl lemma, then implies that the weak limit in L^{0,T; L^{flo x Do)) of ^'(V'^*'"'^'') is 
equal to $(^) as L — >■ oo, with At — o(L^^). Hence, by the argument, outlined in the previous 
paragraph, we deduce almost everywhere convergence of a subsequence on (0,T) x ilo x Dq, and 
finally, using an increasing sequence of nested Lipschitz subdomains (0, T) x fife x Dk, fc = 1, 2, . . . , 
and extracting a diagonal sequence from ?/'^**''^'', we arrive at a subsequence of {"0^ **^'^''}i>i that 
converges almost everywhere on (0, T) x x Z3 to as L — > oo, with At = o{L^^). Since the 
set J) := (0, T) X X 13 has finite measure, according to Egoroff's theorem (cf. Theorem 2.22 on 
p. 149 of |25j ) almost everywhere convergence implies almost uniform convergence, and in particular 
convergence in measure. Thus, by Vitali's convergence theorem (cf. Theorem 2.24 on p. 150 of |25j). 
and thanks to the uniform integrability of the sequence {4'l*^^^^}l>i in -^^(0, T; L\j{fl x D)), we 
finally deduce the desired strong convergence of the sequence {^/'^ }l>i in L^(0, T; L\j{fl x D)) 

below to strong 



4.2 



as L oo, with At = o(L ^). We will further strengthen this by using Lemma 
convergence in the LP{0, T; L\j{il x D)) norm, for any p €E [1, oo). 

We now embark on the programme outlined above by observing that, since on each compact 
subset of D the Maxwellian M is bounded above and below by positive constants (depending on 
the choice of the compact subset), it follows from (4.57a) that {V'l**''^''}-L>1j with At = o(L^^), 
is weakly relatively compact in Ll^^{n x D x (0,T)). Hence, by uniqueness of the weak limit, 

(4.67) V-f'^'^^^^V' weakly in ii^^(0,T;Li(f7 X D)). 
We shall show that in fact 

(4.68) ip^*'-'^'^ ^ a.e. on (0,r) X X D. 

Let O := ft X D , and suppose that Oq is a Lipschitz subdomain of O such that Oq (s O. As 
s \og{s + 1) + 1 > s for all s e M>o we have from (4.65 1, the bounds on the seventh and the eighth 
term on the left-hand side of (4.53), and noting once again that M is bounded below on Oq by a 



positive constant (which may depend on Oq), that 
(4.69) max 



WIII/i(Oo)di<t^(Oo), 



where C{Oq) is a positive constant, which may depend on Oq but is independent of L and At. 
It then follows from the bound on the second term on the left-hand side of (4.69) by Sobolev 
embedding applied on the bounded Lipschitz domain Oq <e O C that 



(4.70) 



LiK + l)d-2 



dt< (7(00). 



Interpolation between the first inequality in (4.69) and the inequality (4.70) then yields that 

uf^^^\t)\\ dt= iiv/^f('±)(t)ii ^<™,, dt<c(Oo). 

L (K + i)d (Op) Jo L (K+i)d (Oq) 



(4.71) 
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By writing ip^^^'^^t) = [y V'l *^'*'']^ and applying Holder's inequality we then deduce for any 
p£ [1,2) that 

i-T 





2P ^,iK+l)d- 



L^-p (Oo) 



dt 



P 



{K + l)d 



Let po be the largest number in the range [1,2) that satisfies (4.731; we thus deduce from (4.72| 
that 

i-T 



(4.74) 



dt < C{Oo), 



(K + l)d + 2 , , 
with Po := TT^T-TTTT-T e (1'2). 



(is: + l)rf+l 



Thanks to (3.771 we also have that 



IIV'L*^''''^llL°°((0,T)xf1;Li,j(Li)) < C', 



and therefore, 

(4.75) ||V^f^'^^||L~((0,T)xO„;LMA,)) <C 

for any two Lipschitz subdomains[^ilo <£ ^ and d _D; here C is a positive constant, independent 
of L and Ai. Fixing Oq = f^o x and interpolating between (4.75) and (4.71), which states that 



(g+l)d+2 

i Cif+i)d ((o,T)xOoxZ5o) 



we deduce that for any two real numbers q\ and 92, with 

2 

< (71 < cxD and 1 < 92 < 1 



and satisfying the relation 



[K + l)d 



91 1 



1 

92 



(i^ + l)rf' 



(4.77) 

we have that 
(4.78) 

Note further that since p^*'^ > Pmin a-e. on x [0,r] and /i(p^*'^) > Mmin a.e. on i7 x [0,T], 
interpolation between the bounds (cf. (4.53)) 
(4.79) 



lV'L*^''^^l!L'!i((0,T)xao;i'2(i5o)) - C'(Co)- 



max 



[B(yo,/,V'o)]' 



l^?(yL*^'*^)WI'dxdi < 



[B(yo,/>o)] = 



yields, using Korn's inequality (3.40), that 

(•T 

At(,±) 



(4.80) 



\u. 



L=(f2) 



dt < C, where 



s < 



^ = 4 when d = 2 
when d = 3, 



^ 2(d+2) ^ 10 
d 3 



^ Let us suppose that O := f2 X D, where C and D are bounded open Lipschitz domains in and R" 
respectively. Then, O is a bounded open Lipschitz domain in R™+". This follows on noting that the Cartesian 
product of two bounded open sets is itself is open. That O is a Lipschitz domain follows by combining Theorem 
3.1 in the Ph.D. Thesis of Hochmuth I28| . which implies that the Cartesian product of a finite number of bounded 
domains, each satisfying the uniform cone property, is a bounded domain satisfying the uniform cone property; and 
Theorem 1.2.2.2 in the book of Grisvard |26l . which states that a bounded open set in R" has the uniform cone 
property if, and only if, its boundary is Lipschitz. 
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At,-il+5 



1+.5 



dq dx dt 



and C is a positive constant, independent of L and At. Hence, by Holder's inequality, (4.80) and 
( |4.78[ ), we get that 

i-T 

(4.81) < C„ 





TAt.+ 


1+5 


UDo 







dx dt 



At, 



l+<5 



dq I dxdt 



where (5 > is to be chosen, and 1/a + 1/b = 1, 1 < a, b < oo, with 



(l + (5)a< + 



To this end, we define 



and we select any 



r := 



91 > 1 



{K + l)d 
d 



1 



Note that qi > l + r = l + 2/{{K + l)d); and in particular qi > r. We then define q2 ■= I + S, 
where S :— r/{qi — r); hence qi > q2- We let a := qi/{qi — 92), b 91/92- Clearly, with such a 
choice of qi and 52, we have that qi{l — (I/92)) — r and l<q2<l + ^ = l + '^./{{K + < qi. 
We thus deduce from ( |4.81[ ) using ( |4.80| ) and ( [4.78^ that 

(4.82) llyf' 

where 5 > is as defined above; Co = ^0 x ^^o; and C(C'o) is a positive constant, independent of 
L and At. 
Analogously, 

K 



Ci'^*^ V'l*'^IIli+'5((o,t)xOo) - C'(Co), 



(4.83) 



i=l 



At. 



c(pf^-)/3^(v.r'") 



At,+ - 



< C(Oo), 



Li+*((0,T)xOo) 



this time with < S < r/{r + 2), where, as above, r = 2/((K + l)d). This follows on noting the 
second inequality in (4.79), (3.40) and that, thanks to ( 4.78[ ), 



I TAt(,±) 



< \\->jjL*'"'^^\\L''i{{0,T)xno\Li2(Do}) < C{Oa), 



L'!i((0,T)xOo;L92(_Do)) 

with qi = 2(1 + (5)/(l — 5); % = 1 + 6 with < S < r / {r + 2); q2 ~ $2; 9i related to q2 via 
(4.77), and noting that since < r/(r + 2) < r < 1, we have I + r < qi < qi = r(l + S)/S < 00, 
I<q2=q2<l + r. Here, again, (7(00) is a positive constant, independent of L and At. 

With the bounds we have established on ipf^*^^'^\ we are now in a position to apply the Div-Curl 
lemma, which we next state (cf., for example, ^24], p. 343, Theorem 10.21). 



Theorem 4.2. Suppose that 2) C is a bounded open Lipschitz domain and E N>2. Let, for 
any real number s > 1, W~^'''{'D) and W~^''' (Tl;M/^^'^) denote the duals of the Sobolev spaces 

^^^), respectively. Assume that 



where 



weakly in LP(£>;M^), 
weakly in L'^{T)\] 



111. 

- + - = - < 1. 
p q r 
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Suppose also that there exists a real number s > 1 such that 

div ffn = V • ffn is precompact in W^^'^{'D), and 

curl g„ = (y9„ - (ygn)'^) is precompact in W^-i^''(D; M'^^'^). 

Then, 

Hn • Qn ^ H ' Q weakly in U'{1)). 



Remark 4.2. Concerning the extension of Theorem \4.^ to the case of r — 1, we refer to the recent 
paper by Conti, Dolzmann and Muller |16| . For the our purposes here the version of the Div-Curl 
lemma as stated above will suffice. 



Consider the following sequences of *Tt = I + d + Kd component vector functions defined on the 
Lipschitz domain D :^ {0,T) x Qq x Dq C R'^ , 



H 



AtX 



H 



it) 



H 



{x,l) 



At,L ' ^^At,L 



{x,d) . jj-(qi,l) 



(91, d) 



' ^^At,L ' ^^At,L 1 ■ • ■ ' ^^At,L I ■ ■ • I ^^At,L ; • ■ • ; At,L 



H 



(<1K,1) 



{qK,d) 



where, 



H 



it) 

AtX 

{X,l) jj(x,d)s 

At,L' ■ • • ^^Atx) 

iq,A) TTiq,,d)~^ 

At,L ' ■ • ■ ' ^At,L 



At 



M(C(pl)^l) 
Ul MCI V-L 



K 



1 = 1, 



and 



Q 



At.L 



(1 + ^1*)",0,...,0 



with a e (0, i) fixed. 



d + Kd 
times 



Thanks to ( |4.78| , ( |4.72| , ( |4.82| and ( |4.83p , there exists a real number e (1,2) such that the 
sequence {i?At,L}L>i, with = o{L^^), is bounded in i^* (S, M^); hence there exists an element 
H € and a subsequence, not indicated, such that HAt.L ~> -H^, weakly in L^* {D,M.^). 

Also, the sequence {QAt,L}L>i, with At = o(L~-^), is bounded in L"^' {D,R'^) with = 1/a; 
hence there exists a Q € L''*(2),M^) and a subsequence, not indicated, such that QAt,L Q, 
weakly in L'?*(D,M^). 

With our definition of {i?At.L}L>i, we have that 

divt,j;,, HAt,L = 0. 

Therefore the sequence {divt.x,q HAt.L}L>i, with At — o(i~^), is precompact in W~^''^{D) for all 
s > 1. 

Further, since a £ (0, it follows from (4.691 that {QAt,L}L>i satisfies 



I curl 



(0,T)xOo 



t,x,q ' 



]At,Lrdqdxdt < C 



IVx.qii + ^tTrdq dxdt 



(0,T)xOo 



V:r,,\/VL I d? dt 



< c 

'(0,T)xOo 

< C{Oo). 

Therefore, the sequence {curl^ j. ^ (5At,L}L>ii with At = o(L^^), is precompact in the function 
space VF-i'2(j);M^nxm)_ 

We thus deduce from Theorem 14.21 that 



HAt,L ■ Qaix — HAt,L ■ QAtx, 
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where the overhne • signifies the weak hmit in (J)) of the sequence appearing under the overline; 
thus, 



(4.84) 
As 



(C(pl) V^l)^* (1 + )" = (C(pl) ^l)^* (1 + V-f 



(C(pL)^Lr*(-,0-C£(-) 



t-t 



n-1 



tji t 



At ■ At '^'^ 



+ (cr'(-)-c£(-)) ^r\-) 



tji t 

~Ar 



for all t e [tn-i, tn] and n — 1, . . . , N , which in turn implies that 

(4.85) iCiPL) = C(pf '+) V^f + (C(pf ^") - C(pf ^+)) V^f ^A*, 

where is the nonnegative discontinuous piecewise hnear function defined on (0, T] by 



^At(i) 



At 



1, 



the fact that, by (4.55el, 

At,-N At 



C(Pl ' ) ^ C(Pl ' )llL=°(0,T;LP(f2)) 



< IIC(Pl*^ ) - C(p)llL-(0,T;LP(n)) + I1C(Pl'''^) ^ C(p)llL-(O.T;LP(n)) 



At, + \ 



as i — 7- oo (with At = o{L ^)), 1 <p < oo, implies, on noting (4.71), that 
(4.86) ||(C(pf •~)-C(pf '"^^aJlm^)) ^0, 

as i cx) (with At = o(L-i)). Further, by ( |4.55e[ ), C(Pl*'^) ^ C e i°°(0,T; L°°(r2)) strongly in 
L'^{0,T;LP{n)), 1 < p < oo, and, by (|467|, V'l* -> V' weakly in Li,^(0, T; Li(r2 x D)); thus we 
deduce, on noting (4.71), that C(Pl*'^) V'l* converges to C,{p)il>, weakly in L^(£l). Hence we have 
shown that 

(C(pl)^l)^* - C(pf'+)V^f = C(P) ^ e Li(J)). 

Consequently, we have from (4.84 1 that 

(C(PL)^^L)^*(1 + V^f)"^ 

Noting ( |4.85[ ) and ( |4.71| ) then yields that 



C(p)V^ (1 + V^f*) 



weakly in L-'^(S)). 



C(Pl '+) ^L* (1 + V-L*)" ^ C(P) (1 + 'Af )" weakly in LH®)- 



Thus, by the strong convergence C(Pl*'^) ^ C{p) e i°°(0, T; L°°(f7)) in L°°(0, T; iP(rj)), l<p< 
oo, which, thanks to our assumptions on the function C stated in ( |3.3[ ) implies that 1/C(Pl*'^) ^ 
l/C(/5) e L°°{0,T;L°°{n)) in L°^{0,T; LP{n)), 1< p < oo, we finally have, on noting ( |4Jl| ), that 

weakly in L^{Ti). 



(4.87) 



As, by definition, (1 + V^f *)" (1 + V'l*)", weakly in ^^(S), by adding this to ( |4.87| we have 
that 



(1 + = (1 + V'L*) (1 + V'l')" ^ (1 + V') (1 + V-L*)" weakly in L\D). 

Thanks to the weak lower-semicontinuity of the continuous convex function s G [0, oo) i— > s"+^ G 
[0,oo) it follows (cf. Theorem 10.20 on p.339 of ^) that 



Consequently, 
(4.88) 



(1 + ^)1+" < (l + ^)(l + ^At)c 
(1 + ^)" < (l+^At)a, 



On the other hand, the function s S [0, oo) i— >■ s" G [0, oo) is continuous and concave, and 
therefore s G [0, oo) i— >■ — s" G (— oo, 0] is continuous and convex; thus, once again by the weak 
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lower-semicontinuity of continuous convex functions, we immediately have (cf. Theorem 10.20 on 
p.339 of [Ml) that 



(4.89) 



(1 + V')"<-(1 + V'f)". 



We deduce from ( |4.88[ ) and ( |4.89[ ) that 
(4.90) 



(1 + V-)" 



+ )", 



and consequently, since the function s £ [0, oo) —s" e (— oo, 0] is continuous and strictly convex, 
and its domain of definition, [0,cx)), is a convex set. Theorem 10.20 on p.339 of [24^ implies that 
there exists a subsequence (not relabelled) such that 



a.e. in D. 



(4.91) Vi' ^ ^ 

Next, we select an increasing nested sequence 

{^ol^i of bounded open Lipschitz domains 
Dq = (0, T) X X Dq, where {f^oifc^i ^^'^ {^o}T=i ^^'^ increasing nested sequences of bounded 
open Lipschitz domains in fl and D, respectively, such that Ufe^i ®o = (0,T) x O x D. Since for 
each k we have pointwise convergence on 2)q of a subsequence of {■(/'l*}l>i, by using a diagonal 
procedure, we can extract from {V'l*}l>i a subsequence, (which is, once again, not relabelled) 
such that 

-^f"* i) a.e. in (0, T) x x D 
with respect to the Lebesgue measure on (0, T) x x D. Let, for any Borel subset Aoi{0,T)xQ.xD, 



M dq Ax dt. 



Since M E L^{{0, T)xilx D), the measure i' is absolutely continuous with respect to the Lebesgue 
measure, which then implies that ■0^ * ~^ V'j almost everywhere with respect to the measure (or, 
briefly, v almost everywhere) in (0,T) x x D. Since i/((0,T) x ft x D) < oo, according to 
EgorofF's theorem (cf. Theorem 2.22 on p. 149 of |5Fj) ly almost everywhere convergence of -0^* to 
■0 implies v almost uniform convergence of '0^* to ■0, and in particular v convergence in measure 
of -0^* to 0. Finally, by Vitali's convergence theorem (cf. Theorem 2.24 on p. 150 of [IS]), the 
uniform integrability of the sequence {V'l*}l>i in i^((0,T) x fl x D;v) = £^(0, T; L]^^(^^ x D)) 
and I' convergence in measure of "0^ * to "0 together imply that 



(4.92) 



strongly in L^{0,T; Llj{n x D)). 



It follows from the (|4.92|) and the sixth bound in (|4.53|) that 
(4.93) 



7Ai(,±) 



strongly in £^(0, T; L]v^(f2 x D)). 



In fact, (4.93) can be further strengthened: it follows from Lemma 4.2 below and (4.93) that 

strongly in U'{^,T;L\^{^ x D)) Vp G [l,oo). 



(4.94) 



TAt(,±) 



Lemma 4.2. Suppose that a sequence {</?n}5^Li converges in L^{0,T; L\j{fl x D)) to a function 
(fi G L^{0,T; L\,j{fl X D)), and is bounded in L°° {0,T; L\j{fl x D)); i.e., there exists Kq > such 
that ||Vn||L°=(o,r;L],,(OxD)) < Ko for all n>l. Then, ip G LP{Q,T;L]^{Q. x D)) for all p E [l,oo), 
and the sequence {(pn\n>i converges to ^ in LP(0,T; L\j{n x D)) for all p G [l,c»). 

Proof. See the proof of Lemma 5.1 in [TT]. □ 

This then completes our proof of strong convergence of the sequence {'0L*'"'^''}i>i to the func- 
tion G L^{0,T;Llj{n X D)) in the norm of the space LP{Q,T;L\i{n x D)), for ah p G [l,oo). 

5. Passage to the limit L — >• oo: existence of weak solutions to this FENE chain 
model with variable density and viscosity 



We are now ready to pass to the limit with L — >■ oo. 
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5.1. Passage to the limit L -> co. In this section we prove the central result of the paper. 



Theorem 5.1. Suppose that the assumptions (3.3) and the condition (4.54), relating At to L, 
hold. Then, there exists a subsequence o/ {(p^*, m^*, '!/'^*)}l>i (not indicated) with At = o{L^^), 
and functions (p, u, ip) such that 

p e L°°(0, T; T) n C([0, T]; if u G L°°(0, T; n ^^(0, T; F), 

where p ^ [l.oo), and 

with ijj > a.e. on fl x D x [0,T], satisfying 
(5.1) 

M{q) i^(x, q,t) dq < ess.sup^^^i — — -—r -0o(x,g)dg) for a.e. {x,t) £ n x [0,T], 
whereby "ip G L°° (0,T; L\j(fl x D)); and finite relative entropy and Fisher information, with 

(5.2) F(4!)£L°^{Q,T-L\j{nxD)) and e L^{Q,T-H\j{Vt x D)), 

such that, as L oo (and thereby At — > 0+ ), 

(5.3a) P"t^^ P weak* in L°°{0,T;L°°{n)), 

(5.3b) pf'\ pt'^'^\ p[^'^ ^ p strongly in L^{Q,T-L^{n)), 

(5.3c) /i(p^**'=^') ^ Ai(/o) strongly inL°°{Q,T- LP {^)), 

(5.3d) C(pr'),C(pf^'*^Cr*^-^C(p) strongly in L^{{\T- LP {n)), 

where p £ [1, oo); 

(5.4a) weak* in L°°{{),T;L^{Vl)), 

(5.4b) uT^'^^ u weakly m L^{0, T; V), 

(5.4c) Ml*^'"^^ ^ w stronc/Z?/ in L2(0, T; L''(17)), 

where r G [1, oo) if d = 2 and r G [1,6) if d — 3; anc? 

(5.5a) M^V^^/ij^^^ M^V^\/^ weakly m L^{0,T; L^{n x D)), 



(5.5b) M5 Vg^i/'L*^''^^ ^ V,y?A weakly in L^{0,T;L^{nx D)), 

(5.5c) strongly m LP {0,T;Llj{nx D)), 

(5.5d) /3^(V^f'*^'=^') ^ strongly m LP {0,T;Lli{nx D)), 
for all p G [l,oo); anrf, 

K if 
(5.5e) V, • ^ C,(M C(Pl*'+) ^l*'+) "> V, • ^ C,(MC(p) ^) weafcZ?/ m 2.^(0, T; y'). 

i=\ ~ ~ i=i ~ 

The triple {p,u,ip) is a global weak solution to problem (P), in the sense that 

— [ [ p-^dx&t— [ [ pu -V xTjAxAt ^ [ pq{x) r]{x,0) dx 
Jo Jn c)t ~ Jo Jn ~ ~ ~ Jn ~ ~ 

(5.6a) Vt? G W^''^{0,T;W^-^in)) s.t. ri{-,T) = 0, 
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with q e (2, oo) when d = 2 and q € [3, 6] when d — 3, 



Jn 



p u ■ Ax dt 

~ at ~ 



Jn 



fi{p) D{u) : D{w) — p{u® u) -.y 



dx dt 



Pq{x) uo{x) ■ w{x, 0) dx + 



Jn 



K 



pf ■w-k^C,{MC{,p)i^) -.V^^wdx 



1=1 



dx dt 



(5.6b) 

and 



Vw e W^'^{Q,T;V) s.t. w{-,T) = 0, 



JnxD 

T 



MC{p)ip-^dq dx dt 



M 



JnxD 



eVj.^ — u({p)'ip ■ Vx^dqdxdt 

\ I I M V^^ V^ • \/q^ip dq dxdt 

Jo JnxD,, ~ ~ ~ ~ 



K K 



Jo JnxD ,_i 



i=i j=i 

K 

(j{u)q 



nxD 



(5.6c) 



C{p)'ip ■ "^qi^ dqdxdt 

M{q) C{pa{x)) ifo{x, q) ip{x, q, 0) dq dx 

yip e W^'\0,T;H'{n X D)) s.t. Lp{-,-,T) = 0, 



with s > 1 + i [K + 1) d. In addition, the weak solution {p, u, ip) satisfies for all t G [0, T]: 



(5.7a) 



\p{t)\Pdx= f \po\Pdx, 
~ Jn 



for p G [1, oo), and the following energy inequality for a.e. t Cz [0, T] 



p{t)\u{t)\''dx+ p{p)\D{u)\^dxds+ 2k M(:{p{t))F{il:{t))dqdx 



Jn 



8ke 

/o JnxD 
< I po \uo\'^dx + 



iV^VVl^dg da;ds- 



Oo k 



2 I rmax^^c 



/||i.(f,)d.s + 2/c 



nxD ~ ~ 

[ M\Vg^\^ dqdxds 
JnxD ~ ~ ~ 

MC{po)J'i^o)dqdx 
nxD ~ ~ 



(5.7b) < [B(uo,/,^o)]' 



with J^{s) = s(logs — 1) + 1, s > 0, and [B(uo, /, i/'o)]^ as defined in (4.26b). 

Proof. We split the proof into a number of steps. 

Step A. The convergence res ults (|5.3a -d), ( |5.4a| -c) and ( 5.5a[ b) were proved in Theorem 4.1 



The strong convergence result (5.5c) was established in Section 4.5 see (4.94). Next from the 
Lipschitz continuity of /3^, we obtain for any p £ [l,oo) that 

IIV'-/3^(V'f*^''^^)llLP(0,T;Li,(nxD)) 

< U-l^'-mLHO.T-.LUnxD)) + IW'^i^) - /3''{^t'^-^^)\\LH0.T;Ll,inxD)) 



(5.8) 



>,T-L\^(nxD) 



< 11-0-/3 W\\LP{0,T;LHnxD)) + H-i'L ILp(0,T;L1,(Ox_D))- 

The first term on the right-hand side of (5.8) converges to zero as L ^ oo on noting that /3^('0) 
converges to ^p almost everywhere on fix 13 x (0, T) and applying Lebesgue's dominated convergence 
theorem, the second term converges to on noting (5.5c). Hence, we obtain the desired result 

As the sequences {'>Pl'^'"''^^}l>i converge to ip strongly in L^(0, T; L\j{fl x D)), it follows (upon 
extraction of suitable subsequences) that they converge to a.e. on x D x [0,r]. This then, 
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in turn, implies that the sequences {-7-'('0^*'''*'')}l>i converge to J-{ip) a.e. on Q x D x [0,T]; in 
particular, for a.e. t e [0,r], the sequences { J"(V'^'^'^^ (■, ■,i))}L>i converge to J-{tp{-, ■,t)) a.e. on 
ri X £). Since is nonnegative, Fatou's lemma then implies that, for a.e. t e [0,r], 

M{q)T{i^{x,q,t))dx dq 



nxD 



(5.9) 



< lim infi- 



M{q)J'{^^'^'^\x,q,t))dxdq < C*, 



nxD 



where the second inequality in (5.9) stems from the bound on the fifth term on the left-hand side 
of (4.53). Hen ce the first result in ( 5.2 ) hold s, and the second was established in The orem 4.1 
Similarly, (5.1) is established on noting (3.77) and that V'l^'*'' ^ 0- Analogously to (5.9), one can 
establish for a.e. t e [0, T] that 

M{q) C{p{x)) J^{4i{x,q,t)) dxdq 

< lim infi_>.c 



nxD 



(5.10) 



nxD 



M{q) ({pT^^^^ i^)) J'aiijt'^'^^x, q, t)) dx dq. 



Finally, it is shown in Step 3.7 in the proof of Theorem 6.1 in 9J on noting (3.15 ), (5.5b c) and 



(|5.1|), that in the case C = 1 



(5.11) fc/ V / Ci(AfC(Pi*'+)V'L*'^) : V^wdxdi^fc / V / C,{M C{p)^p) -.S/^wdxdt 



as L — ^ oo, for any divergence-free function w £ C^{[0,T];C^{Q)). The proof there is easily 
generalised to the present variable C, on noting (5.3d). This implies (5.5e), thanks to the denseness 
of these smooth divergence functions in the function space L'^{0,T;V), and on showing that the 
right-hand side of Ku\ is well-defined for w € L^{Q,T; V), on noting (|3.15l), K2\ and 



Step B. We are now ready to return to (4.3a-c) and pass to the limit L — >■ oo (and thereby 
also At 0+). We shall discuss them one at a time, starting with equation (4.3a). We have 
already passed to the limit L — )■ oo in (4.3a) using (5.3b) and (5.4c) to obtain (4.58) in the proof 
of Theorem 



4.1 



The desired result (5.6a) then follows from (4.58) on noting the denseness of 
the set of all functions contained in C^{[0,T];W^'^ {fl)) and vanishing at i = T in the set of 
all functions contained in W'^'^{0,T;W^ (ft)) and vanishing at t = T. In addition, the energy 



equality (5.7a) was proved in the proof of Theorem 4.1 see (4.59) 



Step C. Having dealt with (4.3a), we now turn to (4.3b), with the aim to pass to the limit with 
L (and At). We choose as our test function 

(5.12) w e C\[0,T];C^{n)) with w{-,T) = 0, and • w = on Q for all t e [0,T]. 



Clearly, any such w belongs to L^(0,T; V) and is therefore a legitimate choice of test function in 
(4.3b). Integration by parts with respect to t on the first term in (4.3b), and noting (3.30) and 



(4.2) for the second term yields that 
Jo Jn 



(pl ul)^^ ■ -^dxdt 



dw 



Jn 



pI^'^*"^ u^*' • Va;(w^*'^ • w) dx di 



Jn 

T 



ti{p^'^+)D(u'^'+):D[w)dxdt 



Jn 



{At} 

Pl 



[u 



At,. 



v7 ^ At,+' 



I At,- 



At,-|- 



dx dt 
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/o Jn 



= f pou°-wiO)dx+ r f pf''+ f^*-+ -wdx 
Jn ~ ~ ~ Jo Jn ~ ~ ~ 



dt 



(5.13) 



1=1 



Jn 



c,{MapT'^)i'T'^)-^xwdxdt. 



Next we note from (4.1a b) that for t e (t„_i, t„], n = 1, . . . , iV 
(5.14) {pLULt'=pt'u 



.+ (tn - t) 



At,- ( At,+ At 



Pl ( 



Air,— \ 



It foll ows that we can pass to the limit L — > oo, (Ai -> 0), in (5.13), on noting (5.14), (5.3b c), 
(5.4b c), (5.5el, (3.28) and that converges to uq weakly in iJ, to obtain 



"'O 



pu ■ 



dw 
di 



dx dt 



Jn 



Jn 



p{p) D{u) : D{w) dx dt — ^ / / pu ■ \/ x{u ■ w) dx dt 

Jq Jn ~ ~ ~ ~ ~ 

(u ■ V x)u ■ w — {u ■ V.j.)w^ ■ wj dxdt 



1 K T 

(5.15) = / pouo-w{0)dx+ I I pf-wdxdt^kV] [ [ C,(Af C(p) V') : w dx dt. 

Jn ~ ~ ~ Jo Jn ~ ~ ~ -^-^ Jo Jn~ » ~ ~ 

The desir ed res ult (5.6b) then follows from (5.15) on no ting (13.24 1, the denseness of the test 
hmctions ( |5.12[ ) in W^O, T; V) and that all the terms in ( |5.6b ) are well-definec . 

Step D. Similarly to (5.13), we obtain from performing integration by parts with respect to 



time on the first term in (|4.3c|) that 



JnxD 



-Af (C(pL)^i 



dip 1 



K K 



1=1 j=i 



dq dx dt 



IQ J(lxD 



r 7irv7 TAt,+ At.- ,,r/.{At} ~At. 



Vx^dq dx dt 



Jo JnxD - 1 



/ At, + \ 



C(pf'+)/3^(^f' + ) • V,,vdqdxdt 



(5.16) 



AfC(/3o)V'V(0)dgda; \/p C^{[Q,T]-C°° {VL x D)) s.t. ip{-,-,T)^Q. 



We now pass to the limit L — > oo (and At — > 0+) in ( 5.16[) to obtain (5.6c) for the smooth ip 
of (|5.16 ) using the convergence results (5.3dl, ( 5.4b[ c), (5.5a-d) and ( |3.20c ). The desired result 
(5.6c) then follows on noting the denseness of the test functions ip £ C^{[0, T], C°°{fl x D)) with 
(/?(•, •,T) = in W^^^{0,T;H''{n x D)) with pi-,-,T) = 0, for s > 1 + i (JsT + 1) d, and that all the 
terms in (5.6c) are well-defined. 



Step E. The energy inequality (5.7b) is a direct consequence of the convergence results (5.3a~d), 
(5.4a b) and ( 5.5a[ b), on noting (3.28), (5.10) and the (weak) lower-semicontinuity of the terms 
on the left-hand side of ( 4.26b| ). For example, it follows from (4.26b) for a.e. t S [0,T] that 



[p^*'+(t)]5 u^*'+(t) ^ 5(t) weakly in L^(r2) as L -y oo 



(5.17) 



hminf/ pt''^{t)\u2''^{t)Ydx> / \g{t)\'dx. 



It follows from (5.3b) and (5.4a) that g{t) = [p{t)]^- u{t) . 



□ 
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Appendix A. Equivalence of notions of solution to the continuity equation 



The purpose of this Appendix is to show that the notion of solution that was used in (3.25), 



(3.26a) coincides with the notion of distributional solution to the continuity equation used by 
DiPerna & Lions [20], whose work we referred to following (3.29) for a proof of the existence of a 
unique solution to (13. 251), (|3.26a|). 
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As previously, we shall assume throughout this Appendix that is a bounded open Lipschitz 
domain in M'', d € {2, 3}. We shall suppose that h e LF'{Q,T\V) and pq ^ T and we extend 6 
by outside VL x (0,T) to the whole of E'' x (0,T). The resulting function, denoted by 6*, then 
belongs to L'^{Q,T]H'^{W^)) and is divergence-free on the whole of x (0,T). We define pi as 
being equal to po on Vl and to in M*^ \ SI. Clearly, pi € LP(]R'*) for all p e [1, oo]. 

We consider the transport equation 

dp* 
~dt 

subject to the initial condition p*(-,0) — Po(')- This initial- value problem is to be understood in 
the sense of distributions of M'' x (0,r); i.e.. 



+ Yx-{b*p*) = 



l'^ X (0,T), 



(A.l) 



P at 
• at 



k* P* -YxV* dxdt = 



for all test functions 77* E C°°(M'' x [0,r]) with compact support in M'' x [0,T); we wiU denote 
this test space by ViW^ x [0,T)) and note that V{R'^ x [0,r)) = C^([0,T); C^{R'^)). 

According to Proposition II. 1 in the work of DiPerna & Lions [50], there exists a solution 

I''')) corresponding to the initial condition 



p* to the initial-value problem (A.ll in L°°{0,T;. 
pI e LP{R'^) for each p £ [l,oo]. Furthermore, by Corollary II. 1 in [20] the solution of ( A.l[ ) is 
unique in L°°{0, T; LP(M'^)) for all p € [1, 00], and by Corollary II.2 in [20], p* € C([0, T];LP{R'^)) 
for all p e [1, 00). 

Since &* = in (M'^ \ H) x (0,r) and = in M'* \ n, trivially, p* = in (M'' \n)x {0,T). 



Therefore, (A.ll is equivalent to the following: 



(A.2) 







p -P- dxdt 
' dt ~ 



PoVix,0) dx 



n 



pb ■ YxV^xdt = 



for all test functions 77 e C^([0,T); C°°(il)), where p denotes the restriction of p* to fl x (0,r). 
Thus, the solution p to (A.2 1, subject to the initial datum po € T, exists and is unique in 
L°°{0,T;LP{n)) for all p e [l,oo], and it belongs to C([0, T]; for all p e [l,oo). 

By selecting -q G Cq°((0, T); C°°(r2)) in (A.2), the second term on the left-hand side vanishes, 
and by noting that, thanks to the Sobolev embedding theorem, iJ^(ri) is continuously embedded 
in L'^{n) for q S (2, c») when d = 2 and for q E [3, 6] when d = 3, we then have that 



(A.3) 



A A^ 

P ^ d J dt 
ot 



< c{q) ||p||L~(o,T;L-(n)) \\b\\L^o.T;LHn)) WYxVW 



L'^{0,T;Li-i (n)) 



for all 77 E C(^{{0,T);C°°{rt)), where c{q) is a positive constant; and hence, by the density of 
C°°(n) in W^'^in), also for ah t] E C(f'((0, T); 1^^'^ (17)). 

Let us consider, for b E L^{0,T;V) and po E T fixed, and therefore p E L°°{0,T; L°°{n)) also 
fixed. 



By (A.3I, 



m\<ciq,P,b)M 



Pg^dxdt, TjEC^iiO,T);W''^m 



1^ \fr]EC^((0,T);W^'^Un)). 

By the Hahn-Banach theorem the continuous linear functional i can be extended from the linear 
subspace C^((0, T); W^^'?^ (f])) of ^^(o, T; VF^'^ (fi)) to a continuous linear functional on the 
whole of L2(0,T;I1K1'^ (!])), and since Co°°((0, T); IF^'i^ (f])) is dense in L^{0,T-W^^^^ {n)) 
the extension is unique; we denote it byr. Hence, ^(77) = £(7;) for all 77 e C^f ((0, T); M/^^'^ (17)), 
and 

ir (77)1 < c{q, p, b) V77 E L\0, T- W^^^^m, 

and therefore t E L'^{{),T;W^^^ {9))' . Since for the range of q under consideration both 
W^'^(O) and W^'T^^{9)' are reflexive Banach spaces, it follows from a result of Bochner & 
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Taylor (cf. Remark 2.22.5 on p. 125 in Kufner, John & Fucik [3T]; and a result of Gel'fand cited in 
item (b) on the same page in [31]) that L^(0, T; W^''^ {il)) is a reflexive Banach space and there 
exists an element g e L'^{0,T;W^''^ (V,)') such that 



In particular since p ^ is an integrable function on ft, we have that 



drj 



dt 



w 1-'^ (n) 

T 



dt = -l{r]) ^ -t{ri) 



Consequently, by selecting test functions rj of the form r/ — if{t)£^{x), where ip £ C§°{0,T) and 
we deduce that 



and therefore, by the bilinearity of the duality pairing 



also 



dip 



dt = 



W '9-1 (O) 



W '9-1 (f2) 



By applying Corollary |3.2| to both sides of the last equality, noting in particular that p and g both 
belong to L\0,T;W^'^{ny), the latter since g G L^{0,T;W^'^ {n)') C ^^(O, T; VF^'i^ (f^)'), 
and the former since p e L°°(0, T; L°°(f])), and W^'^(n) C ^^(f^) implies L°°{n) C 
we deduce that 







dt 



w '9-1 (n) 



Hence, 



p(t)^(t)dt 



5vdt,0 ^ V^eCo°°(0,r), veeiyi'^(r!). 

/ VK^'?^(f2) 



g(t)(^(t)dt V^eC6"(0,T), 



3.1 



as an equality in the Banach space W^''^ {fl)' . By the implication (ii) ^ (i) in Lemma 
then deduce that p is almost everywhere on (0,r) equal to the primitive of g; i.e., g = dp/ dt. 
Therefore, 



dp 
dt 



dt \f7j e L^{0,T;W^'^^{n)), 



w '9-1 (a) 

and (reverting from the notation dp/ dt to dp/dt), in particular, 



Jn 



p^dxdt = eiTj)^tirj) = 



dp 
dt 



W '9-1 (O) 



dt \/i]eC^{0,T;W^^^^{n)). 



We now substitute this into (A. 2) to deduce that 
(A.4) 



/ (^,v) , dt- / [ pb-V,rjdxdt = yr^eC^iO,T;W'^^m. 
JQ \ /w^-^{n) Ja Jit 



Hence, by density. 



(A.5) 



/o \df 



dt~ [ [ pb-VxVdxdt^O \f7] € L^{0,T;W^'^^{n)). 
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In particular if 6 G L°°{0,T; V) (as is the case when b G V is independent of t; see, for example 
(3.26a) where u"^^ S V was used as b), the choice of the test function in (A.4| can be further 
relaxed thanks to the density of C^(0,T) in L^{0,T); so (A.4) then implies that 



^,r,\ ^ dt- [ / p6 -y^Tydxdt-O V?/e Li(0,T;M^i'^(r!)), 

and in this case dpjdt e £°°(0, T; VF^'^ i.e., p G H/i'°°(0, T; VF^'^ (1))'). This is the 

notion of solution that was used in ( 3.25[ ), (3.26a). 

Conversely, if p e L°°(0, T; L°°(f])) n C([0, T]; n M^i'°°(0, T; W^'i^ (f])') solves ( [XS ) 

subject to the initial condition /9(-,0) = po(') € T, then, by choosing test functions r\ of the form 
T]{x,t) = with ip € C^(0,T) and i € C§°{n), we have from that 



(A.6) 



9p 



dt 



)dt- /" ip f pb-Vx^dxdt^O V(p e Co°°(0,r), e C^f (17). 
Jo Jn ~ ~ 

Since p e VKi'°°(0, T; M^^'i^ (17)') and so both p and g = dp/ dt belong to ^^(O, T; PF^^?^ (17)'), 
it follows from Corollary |3 . 1 1 that 

pit) = P0 + ^{s)ds, poe L°°{n) cW'-^{n)\ fora.c. [0,T]. 

The implication (i) (iii) of Lemma 3.1 applied to the first term in (A.6) then implies that (A.6) 
can be rewritten as 



(A.7) 



£ f^{p,0^i^^^^^^dt^ J\ l^^pb-Y.^dxdt^O e Co°° (0, T) , e Co" (17) . 



After partial integration in the first term in (A.7) we have that 
(A.8) 



dt- f ip f pb-YxC<ixdt = V(/je C^(0,T), V^e C5"(17), 
Jo Jn 



and then, by rewriting the duality pairing in the first term of (A.8) as an integral over 17 using 
that the product p^ is integrable on 17, 



(A.9) 







p^ dx 



n 



dip 
~dt 



dt 



[ p [ pb-Yx(dxdt=^0 V(/jeC5"(0,T), VeeC5"(17). 
Jo Jn 



Upon rearrangement (A.9) becomes 
(A.IO) 
The set 



-/ I P^^^'^^' I [ pb-YxivOdxdt^o WpeC^{0,T),yieC^{n). 

Jo Jo CIt Jq 

l?]eC^{nx{0,T)) : r]{x,t)^ ^ p,{t)^i{x), for some iV G N 

l<i<N 



and some p^ G C^f (0,T), ^.(x) e Cg°(17), i = 1, 



is dense in C^(17 x (0,T )) (cf. Lemma in Ch 1., ^3, Sec. 2 of Vladimirov [M], for example). We 
thus deduce from (A.IO) with p and ^ replaced by pi and respectively, forming finite linear 
combinations of the resulting equations, and then using a density argument, that 

(A.ll) -[ [ P^dt- [ [ pb-YxVdxdt^O V77 e Co°°(17 X (0,r)). 

Jo Jn Jo Jn ~ ~ 

This implies that p G i°°(0, T; L°°(17)) n C([0, T]; LP(17)) n W^^°° {0,T]W^'^ (ny) satisfies 

dp 



dt 



+ Yx ■ (bp) = 



2?'(17x (0,T)), 
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together with the initial condition p{-,0) — Po(')i where po € T, i.e., that p e L°° {0,T; L°° (il)) n 
C{[0,T];LP{n)) n W'^^°°{0,T;W^'^{ny) is a solution to the continuity equation in the sense of 
distributions on x (0,T). 

Alternatively, on extending 5 by to b* and pq and p by to Pq and p*, respectively, as at the 
start of this Appendix, we deduce from (A. 11 1 that p* satisfies 

% + Y.-ik*P*)^0 in2?'(]R'^x (0,T)), 
ot 

together with the initial condition p*{-,0) — Po(-)i where Pq e iP(M'') for all p E [l,oo], and 
p* e L°°{0,T;LP{R'^)) with p e [l,oo]; p* e C{[0,T]; Lp{R'')) n W'^^°°{0,T;W^'^ {R'^Y) for 
p € [l,c») and for the range of q under consideration. Hence, p* is a solution to the continuity 
equation in the sense of distributions on M'' x (0,T), as in [50] , 
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